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PRECEDING page blank not filmed 


ABSTRACT 

Partial-reflection data collected for the eclipse of July 10, 1972 as 
well as for July 9 and 11, 1972, are analyzed to determine eclipse effects 
on £ -region electron densities* The partial -reflection experiment was set 
up to collect data using an on-line PDP-15 computer and DECtape storage* 

Except for a couple of changes, the experiment was the same setup as used 
by Birley and Seahrist [1971] . The electron-density profiles show good agree- 
ment with results from other eclipses* The partial-reflection programs were 
changed after the eclipse data collection to improve the operation of the 
partial -ref 1 ection system. These changes were mainly due fo expanded 
computer hardware and have simplified the operations of the system consider- 
ably. 
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1. INTRODUCTION 

A solar eclipse can be thought of as the obscuration of solar radiation 
by the intervention of the moon between the sun and a point on the earth. This 
obscuring of the sun is a function of time which varies with the location on the 
earth, altitude, and the type of radiation. Depending on the wavelength of 
solar radiation and the ionospheric constituents, solar radiation can cause three 
chemical processes known as dissociation, ionization and excitation [Whitten and 
Poppoffj 1971] . The variation in solar radiation with time during a solar eclipse 
is given as an obscuration function and varies according to the different wave- 
lengths of radiation. The obscuration function for visible light is easily 
calculated, being just that for the visible disk. Figure 1.1 shows this 
obscuration function for the eclipse of July 10, 1972 at 75 km altitude above 
the University of Illinois Aeronomy Field Station located near Urbana. At this 
location the eclipse was partial, with about 60% maximum obscuration. The 
obscuration functions for various other radiations during a total eclipse are 
shown in Figure 1.2 [Sears, 1972], Notice the large difference between the 
obscuration functions for ultraviolet radiation and X-rays. 

O 

Solar radiation with wavelengths less than 2900A causes various chemical 
reactions in the ionosphere [Whitten and Poppoff , 1971] with the most pronounced 
effects occurring in the D- region (50 to 90 km). For example, Turao and Seehrist 
[1970] show two orders of magnitude change in the electron density and more than 
three orders of magnitude change in CO^ - and C0 4 " at 75 km during sunrise. Certain 
solar radiations greatly enhance the concentration of positive ions as well as the 
electron density so that during the daytime, except for during enhanced particle 
precipitation [Lauter and Knuth s 1967], the main ionization source above 70 km 



PERCENT OBSCURATION 



TIME (CST) 




OBSCURATION FUNCTION 



Figure 1.2 Obscuration functions for visible light (V), Lyman alpha (L ), ultraviolet (UV) , 
and X-ray (X) ionizing fluxes for the 1966 solar eclipse from Sears [1972]. 
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is solar radiation as given in Section 2.1. Therefore by correlating the elec- 
tron densities with the obscuration function for the ionizing radiation in a 
solar eclipse, values for the production and loss of positive ions and confirma- 
tion of the ionizing sources can be obtained. 

Data from the D region have been obtained by both rocket measurements and 
ground-based techniques. Although rocket measurements seem to be more accurate 
[■ Meehtly , et al. t 1967], the amount of data is limited by cost. Ground-based 
techniques can be set up anywhere and can gather large amounts of data, although 
the accuracy is not as great, and they are primarily limited to evaluating elec- 
tron densities. One type of ground-based technique which is discussed in this 
paper is called the partial-reflection experiment. Data are collected using 
vertical incident radio waves which are partially reflected from the D region. 
The information obtained can be in one of two forms: differential absorption 

[Pirnat and Bowhill, 1968] and differential phase [ Wiersma and Seahrist 3 1972]. 
Partial -reflection data using the differential absorption method were collected 
from 1200 to 1700 CST for the solar eclipse of July 10, 1972, as well as July 
9 and 11 as control days. The experiment was set up as described by Bivley and 
Sedhvist [1971] with two exceptions as described in Chapter 3. The solar and 
ionospheric conditions for this experiment are given in Chapter 2. 
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2. PRODUCTION AND LOSS OF THE 2?- REG ION IONIZATION 

Recently several papers have summarized the knowledge of the D region of 
the ionosphere. Thomas [1971] presents an overall review of the D region while 
theoretical models of the D region are presented by Sechrist [1972] , Ferguson 
[1971], Donahue [1972], and Radiceiia and Stowe [1970]. The D region is perhaps 
the most complicated part of the ionosphere as well as the most difficult part 
from which to obtain accurate data. The chemical composition is dependent on 
height and solar zenith angle [Thomas , 1971]; although it consists of neutral 
constituents, positive ions, negative ions, and free electrons, this chapter is 
mainly concerned with the processes of formation and loss of free electrons 
during the daytime (solar zenith angles less than 90°) and during a solar eclipse. 
Using results obtained from measurements on other eclipses, the expected results 
from the partial-reflection experiment are given. 

2 . 1 Ionization Sources 

Although there is general agreement on what ionizes the neutral fi-region 
constituents, there is some doubt as to the relative importance of each source. 

The ionization sources for the daytime D region at midlatitudes, as given by 
Mitra and Rowe [1972] and by Aikin [1972] are: 

O 

1) Lyman-a (1216A) ionizing nitric oxide (NO) 

2) 1-8A X-rays ionizing all constituents 

° i 

3) 1027-1118A ultraviolet radiation ionizing metastable 0 (A ) 

z g 

4) Galactic cosmic rays ionizing all constituents. 

Along with these sources precipitating electrons may be considered 
another source of free electrons, but is of prime importance only in the polar 
regions, at night, or after a magnetic storm [Lauter and Knuth, 1967] and will 
not be considered in this paper. 
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The primary ionization source below 70 km is considered to be galactic 
cosmic rays [Seahristj 1972] , although its effect may extend as high as 75 km 
[Keneshea, 1967], The primary ionization source above 70 km is either (1) or (3) 
depending upon the nitric oxide distribution adopted. Few measurements of nitric 
oxide have been made, so most distributions available are from theoretical models. 
Distributions measured by Barth [1966] and Pearce [1969] are at least an order 
of magnitude greater than distributions calculated from theoretical models of 
the ionosphere [Mitra 3 1966] , but distributions measured by Meira [1971] below 
85 km are about the same as those calculated by Shimazaki and Laird [1970] . 

Using distribution by Barth [1966] for NO, the primary ionization source between 
70 and 80 km is Lyman-a ionizing NO, but using nitric oxide distributions given 
by Shimazaki and Laird [1972] and photoionization rates for 8 iven 

Hunten and MoElroy [1968] , the main ionization source between 70 and 80 km is 
° 1 

1027-1118 A UV radiation ionizing 0 9 ( A ) [ Thomas , 1971]. Somoyajulu and 

* 9 

Avadbanulu [1972] pointed out that according to measurements by Huffman 3 

et at. [1971], photoionization of 0 9 (*A ) is important only above 80 km 

A g 

making Lyman-a the main ionization source. Figure 2.1 from Seahrist [1972] 
shows ion-pair production rates for various radiation during solar minimum. 

In any case the distribution of NO is important to the rate of production of 
free electrons between 70 and 80 km, and the distribution by Meira [1971] is 
used in this paper. 

The variation of ionization sources (1) and (3) with respect to solar 


activity is small [Thomas , [1971], but 2-8 A X-ray flux can change by several 

orders of magnitude. Typical X-ray fluxes for different solar activity as 

. . -3 -2 -1 
given by Avkin [1972] are less than 4 x 10 ergs cm sec for a quiet sun, 

-4 -3 -2 -1 

between 4 x 10 and 4 x 10 ergs cm sec for moderate sun, and greater 
-3 -2 -1 

than 4 x 10 ergs cm sec for an active sun. A solar flare on July 11 at 




Figure 2.1 Ion-pair production rates from various Z?-region ionization sources 
as given by Seahrist [1972] . 
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0 -2 -2 -1 
8:10 AM CST produced a 2-8 A X-ray flux of 1.5 x 10 ergs cm sec . With 

O 

an active sun or a solar flare 2-8 A X-ray ionization can become the primary 

O 

source of ionization. The 2-8 A flux for July 9, 10, 11 in Figure 2.2 from 
the Solar Geophysical Data, 1973 (U. S. Department of Commerce) shows the 
solar activity to be quiet to moderate. The X-ray flux is expected to have 
little or no correlation with the electron density of the upper D region 
except for the X-ray burst near 1435 on July 11. 

2 . 2 Formation of Ions in the D Region 

The electron density between 70 and 85 km is dependent on the formation 
of positive ions. The three main ionization reactions for this region are: 

A) 0 2 + hv -+ 0 2 + + e 

B) NO + hv + N0 + + e 

C) N 2 + hv + N 2 + + e 

as seen in Figure 2.3 adapted from Mitra and Rowe [1972] and Donahue [1972], 
which is a block diagram of the positive-ion chemistry at 75 km. The main loss 
process for N 2 + is by the charge-exchange reaction: 

D) N 2 * ♦ 0 2 *N 2 ♦ 0 2 . 

This reaction is very fast (1 x lO”^ cm^ sec”'*') [ Fehsenfeld , et al. 3 1965], 
Therefore concentrations of N 2 + are small and the production of 0 9 + is either 
by photoionization or by charge transfer. Electron production, therefore can 
be determined by the production of NO and 0 2 minus the formation of NO by 
charge exchange reactions shown in Figure 2.3. Since the production of N0 + is 
dependent on NO distributions, the production rate of the free electrons also 
depends on the NO distribution which can differ by at least an order of 
magnitude (Section 2.1). 

The main positive ions between 70 and 80 km are hydrated ions of the form 
H (H 2 0)^, n being some number greater than zero [Naroisi and Bailey , 1965], 



SOLAR ZENITH ANGLE (deg) 





Figure 2.3 Flow diagram of the formation of positive ions including conversion rates [Donahue s 1972]. 

Three-body rate constants are in units of 10-28 C m6 sec-1; two-body rate constants are in 
units of 10 -9 cm3 sec - l. Rate constants not given by Donahue are from Good , et al, [1970]. 
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Two basic reaction schemes for the formation of water cluster ions as presented 

by Fehsenfeld and Ferguson [1969] are from N0 + and beginning with the reaction 

°2 ++ + M * °4 + + M wliere M is a body. Both schemes are given in 

Figure 2.3. Each scheme raised several questions which are dealt with by Donahue 

[1972]. According to Figure 2.3, N0 + creates hydrates with masses of 55 and 

higher, yet 19 + and 37 + are the dominant hydrates detected. Also the first three 

reactions with NO + are too slow relative to the loss rate. Problems with the 0^ + 

scheme are: it seems to ignore the large NO concentration and the ionization 

of 0-(*A ) seems to be an overestimation according to Huffman , et al. [1971], 

*• . 9 

but this may be the main source of water clusters between 77 and 85 km [Donahue, 

9 

1972]. Even with the large number of hydrated ions, the rapid recombination 
rate competes with the formation of hydrated ions [ Thomas , 1971]. This recombi- 
nation represents the main loss process for free electrons between 70 and 80 km. 

The formation of negative ions would constitute a loss of free electrons 
by the attachment reaction; 

E) e + 0 2 + 0 2 -> 0 2 ' + 0 7 . 

Figure 2.4 by Thomas [1971], giving a scheme for the daytime negative electrons 
at 65 km, shows reaction (E) to be fast, but the loss reactions 

F) 0 2 ‘ + 0 •* 0^ + e 

G) 0 2 " + O^A ) * 2 °2 + 6 

are much faster. Although the formation of 0^ is rapid, there is rapid return 
to 0 2 ". The negative ion chemistry is dependent on atomic oxygen and 
concentrations. At night these concentrations decrease so that reaction (E) 
constitutes an important loss process for free electrons. 

At eclipse totality free electron production is reduced to that comparable 
of nighttime electron production, and the production of atomic oxygen and meta- 
stable 0 7 ( 1 A ) are also greatly reduced [Shimazaki and Laird, 1972], By 
^ 9 




Figure 2.4 Block diagram [Thomas 3 1971] showing the negative ion chemistry during the 
day. The lifetimes of electrons and each ion are for a height of 65 km. 


tsj 
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comparison of eclipse data, Mechtly , et at. [1972] shows the possibility of 
attachment reactions as being the main loss process at totality. This would 



[1965] during an eclipse shows less than an order of magnitude change in atomic 
oxygen. More measurements of atomic oxygen are needed during eclipses to deter- 
mine more accurately the loss process for free electrons during totality of a 
solar eclipse. 

2 . 3 Recombination 

Above 70 ion during the daytime, negative-ion chemistry is not important; 
so the main loss process of free electrons above 70 km is by recombination 
with positive ions , The continuity equation for electrons as given by Whitten 
and Poppoff [1971] is: 



-flL- 

1 + X 


■) - (a 


2 '[«]■* 

dt 




( 2 . 1 ) 


where [ e ] is the electron density, X is the ratio of negative ion concentrations 
to electron densities, q is the ionization rate, is the ion-electron 

recombination coefficient, and a- is the ion-ion recombination coefficient. 

v 

With the assumption that variation in X is insignificant, then dX/dt = 0 
and defining an effective recombination coefficient as = a n + 

Equation (2.1) reduced to: 


die]... , 
dt " 


1 + X" 


a eff 


( 2 . 2 ) 


During a solar eclipse at totality, the electron production decreases by 
several orders of magnitude. Using an ionization rate of zero (<j = 0) , a ^ 
can be obtained from Equation (2.3) for short intervals of time. 
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eff 


. AM 

At 


[ey 


(2.3) 


With small changes in the electron density a ef£ can be obtained by the ap- 
proximation [Mitra and Rowe, 1972] 

“eff = c ? / te] 2 (l + A) (2.4) 

Below 70 km the problem is complicated by the presence of negative ions 
[Mitra and Rowe , 1972] for which a time dependent analysis of the negative 

reaction scheme has to be used [Thomas, 1971]. As discussed in Section 2.2, 
there is the possibility of loss by attachment. Many problems about the loss 
process still remain unsolved including the question of the NO distribution. 

2.4 Expected Results 

Figure 1.2 by Sears [1972] gives the obscuration function for different 
Z>-region solar ionization sources from the eclipse of 1966. Lyman-a and visible 
light have the same obscuration function but not so with UV and X-rays. The 
obscuration function for visible light at Urbana, Illinois for July 10, 1972 
(Figure 1.1) is therefore expected to be different from the obscuration func- 
tion for ultraviolet radiation and X-rays. Using the maps of the sun given in 
Solar-Geophysical Data, 1972 (U.S. Department of Commerce) and the moon's movement 
across the sun's disk, an idea of the obscuration function for different solar 
radiations can be obtained. Since the solar activity during the eclipse was 
quiet to moderate, the predominate ionization source between 70 and 80 km is 
expected to be Lyman-a. 

The total obscuration is about 60%, therefore data is used from previous 
eclipses with a similar obscuration and about the same solar zenith angle. The 
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solar zenith angle is shown in Figure 2.5 to be about 37°. Figure 2.6 by 
Deeks [1966] gives various electron densities for an eclipse during March 
equinox noon at sunspot minimum. Figure 2.7 by Smith , et al. [1965] gives 
electron density distributions for various obscurations of the eclipse of 
July 20, 1963. In Figure 2.6 the electron density for 60% obscuration shows 
little change until above 70 km. For Figure 2.7 at 40% obscuration the elec- 
tron density at 75 km has no change while above and below this altitude show 
marked changes. Below 75 km the change is, therefore, expected to be no 
larger than above 75 km and the change is expected to be approximately 36% 

(from equation (2.4]). Due to the changing solar zenith angle, the magnitude 
of the slope of the changing electron densities before the maximum obscuration 
of the sun is expected to be greater than the slope after maximum obscuration. 
2.5 Statement of the Problem 

The purpose of this paper is to present the setting up, collection, and 
analysis of the partial-reflection data taken before, during and after a solar 
eclipse and to present changes made in the partial-reflection computer programs 
in order to simplify the operation and more effectively reject noise. 



-RAY FLUX (ICTergscm 
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TIME (CST) 

Figure 2.5 The variation of the solar zenith angle for July 10, 1972. The 
partial-reflection data collected period is shown as well as the 
time of maximum obscuration for the eclipse. 
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Figure 2.6 Variation of electron density during a solar eclipse 
at March equinox, mid-day, and sunspot minimum at 
middle latitudes [peeks, 1966] . 



ALTITUDE (km) 



ELECTRON DENSITY (cm* 3 ) 

Figure 2.7 Electron-density profiles for the eclipse of July 20, 1963 [Smiths al , 3 1965], 

Profiles 1, 2, 3, and 4 refer to obscurations of 92%, 86%, 40%, and 2%, respectively. 
The solar zenith angle was 55° at totality and 61° at 40% obscuration. 
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3. EXPERIMENTAL TECHNIQUE 

The partial-reflection experiment was first performed by Gardner and 
Pawsey [1953] . Electron densities were deduced for 65 to 82 km from partially 
reflected, circularly polarized radio waves. The transmitter operated at 1 kw 
during each 30 usee pulse with a center frequency of 2.28 MHz, and the partially 
reflected signals were displayed on an A-scan oscilloscope. Several improvements 
have been made in the experiment and are discussed by Pimat and Bowhill [1968] . 

Gregory [1956] used an increase in transmitter power of 4 kw and a decrease 
in the transmitter pulse width to 9 ysec. These changes improved the amplitude 
and resolution of the partial reflections. Fejer and Vi.ee [1959] developed an 
improved receiving and storing method using a dual-beam cathode-ray tube oscil- 
loscope and camera. The system was operated at 1.83 and 2.63 MHz. Belrose and 
Burke [1964] also operated at two different frequencies (2.66 and 6.275 MHz) and 
transmitter power of 1 Mw, were able to obtain electron densities from the D and 
E region. Belrose and Burke [1964] were the first to use the generalized 
Appleton-Hartree formulas by Sen and Wyller [1960] for partial-reflection 
application. 

Using the generalized Appleton-Hartree formulas and several approximations, 

the ratios of partially reflected extraordinary waves [A ) to the partially re- 

cc 

fleeted ordinary wave [A q ) for two heights can be used to calculate electron 
densities [Pimat and Bowhillj 1968 and Reynolds and Seehrist , 1970] . The ratio 
A l A at each height is inversely related to the absorption by the expression 
exp(2/^ from which the name differential absorption originates . At the 

University of Illinois the electron density was calculated directly from these 
ratios, and as seen in Chapter 4, small changes in these ratios can produce large 
variations in the electron densities. 
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Henry [1966] designed and built the hardware for the partial -reflection 
experiment at the University of Illinois. The transmitter that is presently 
being used was built for the purpose of making shipboard measurements. This 
transmitter operates at 40 kw during each 20 ysec pulse and with 5 pulses per 
second. The center frequency is 2.66 MHz with a 50-ohm unbalanced output. 

Fijjure 5,1 shows a block diagram of the transmitter. The reduction of power 
from the initial 50 kw vised is to give longer life to the tubes used, and the 
pulse is shortened from 50 ysec used by Henry [1966] for better height resolution. 

Figure 3.2 shows the two antenna arrays used to transmit and receive cir- 
cularly polarized signals. Each array consists of 30 half-wave dipoles in the 
north-south direction and 30 in the east-west direction [Wiersma and Sechrist y 
1972] . Each direction has matching networks that differ by 90° from the other 

direction of the same array to give a circularly polarized radio wave as shown 
in Figure 3.3. Each array gives approximately 22 dB gain with the main beam 
in the vertical direction. The first sidelobe is down 14 dB. Since both 
arrays are the same, this is a decrease of approximately 30 dB in the sidelobes 
relative to the main signal which has 44 dB gain. Further details on the 
antennas are given by Pimat and Bowhill [1968] and Reynolds arid Seehrist [1970]. 
3.1 Development of Receiving ccnd Storing Data 

The receiver, storage and timing controls have had two main changes in 
the development of the partial-reflection system. The experiment was originally 
set up using photographic film to store the partially reflected signals as dis- 
played on an oscilloscope (see Figure 3.4). The controlling circuitry or pulser 
sent pulses of 30 volts to the transmitter, receiver, and camera. The pulser 
has remained the same with the exception of the addition of extra control cir- 
cuitry depending on the storage method. The amplitudes of the received signals 
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Figure 3.1 Block diagram of the partial-reflection transmitter 



















Figure 3.2 Partial -reflection antenna arrays 
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Figure 3.3 Block diagram of the partial -ref lection system 
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Figure 3.4 Typical frame of data as collected by Henry [1966]. 
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were later measured visually and electron densities were obtained. Pimat and 
Bcwhill [1968] shows that there is good correlation between electron densities 
calculated from the partial- reflection data and from rocket measurements with 
the transmitter operating at 25 kw of power during a 50 usee pulse. This system 
of collection and storage is inexpensive, but the processing of the data to 
obtain electron densities is very slow and preparation and operation are 
complicated. 

Reynolds and Sechrist [1970] set up data storage on paper tape. Ordinary 
and extraordinary samples were punched on paper tape for heights corresponding 
to 75 km and 80 km. Data can be stored at a rate of 30 values of each sample 
in one minute. From the paper tape the data can then be read into a computer 
and processed. This data on paper can be used to obtain an electron density 
for between 75 and 80 km. Reynolds and Sechrist [1970] show the results using 
paper tape compares favorably with results from rocket measurements and with the 
results published by Belrose and Burke [1964]. Although the system has a faster 
operation than the original system, it produces only one electron density and 
the added control circuitry is very complex. 

Birley and Sechrist [1971] set up the partial-reflection experiment using 
a PDP-15 computer. The received signals were transmitted to the computer via an 
analog to digital converter and stored on DECtape to be processed later. The 
data consisted of four noise samples from 45 to 49.5 km and 21 data samples from 
60 to 90 km in 1.5 km increments. The collection rate is 5 sets of 26 samples 
sec This collection is done alternating between ordinary partial reflec- 
tion and extraordinary partial reflections. Electron densities obtained by 
Birley and Sechrist [1971] show good agreement with electron densities obtained 
from rocket measurements between 67.5 and 82.5 km. The other heights suffered 
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from too many rejections due to noise and saturation of the analog to digital 
converter, small signal to noise ratios, or inaccurate A Ik ratios. Computer 
storage offers several advantages : 

1) A fast rate of data collection (presently limited to the 
transmitter speed) 

2) Data can be stored more compactly and in much larger 
quantities 

3) The controlling circuitry is greatly simplified 

4) The data processing is faster 

5) [e] can be obtained for every 1.5 km 

This type of system also poses several disadvantages: 

1) High cost 

2) Development of computer software 

3) Loss of accuracy in digitizing the data 

4) Development of new circuitry and modification of the old 
for adaption to the A/D converter 

5) More complicated operations (operator must know computer 

operation) 

These disadvantages have been reduced with additional equipment and development 
as given in Section 3.3. 

3.2 Partial-Reflection Data Collection for the Solar Eclipse 

The partial-reflection receiver was interfaced into the PDP-15 computer 
to obtain data to be processed as described by Birley and Sechrist [1971] . 
Several changes in the receiver and controlling circuitry and the addition of 
an analog-to-digital converter were required prior to using the computer. A 
block diagram of the original receiver is shown on page 18 of Aeronomy Report 
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No. 13, [Henry, 1966]. The analog-to-digital converter saturates with an input 
of one volt or greater and will be damaged with inputs greater than five volts. 

The maximum output of the receiver was therefore reduced from 10 volts to 1.5 
volts by one of the IF amplifiers, and the full-wave bridge diode detector 
was replaced by a single diode to reduce the nonlinearity of the receiver. A 
second blanking gate was inserted with the mixer in the RF amplifier module to 
more completely remove the initial effects of the transmitter pulse. The 
polarity reversal circuitry was not used but was left intact while the differen- 
tial amplifier and inverter were replaced by two DC amplifiers on integrated chips. 

The block diagram of the modified receiver is shown in Figure 3.5. Figure 
3.6 shows the RF module with the extra blanking gate and Figure 3.7 shows the 
IF amplifier/DC amplifier module with the revisions. Both modules were modifi- 
cationsof the RF-3 module and IF-6 module respectively, given by Henry [1966]. 

The receiver power supply was unchanged as set up by Henry [1966], Encode pulses 
as shown in Figure 3.8 were used to control the operation of the A/D converter after 
Birley and Seohrist [1971] . The encode pulse circuitry consists of a 5-volt 
power supply and 4 monostable multivibrators (Figure 3.9) with a variable timing 
for length of noise and signal pulses and the delay of each. 

Two main modifications were made in the software set up by Birley and 
Seohrist [1971] . For the first change D. E. Ward [private communication] set up 
a computer-controlled synchronization with the external pulser. The timing 
shown in Figure 3.9 is used to determine which radio wave mode has been received. 
The computer programs are set up to store only pairs of sets of 26 numbers read 
from the A/D converter. A set of numbers is read in and assumed to be from a 
radio wave of ordinary mode. The computer's clock is set for 150 ysec and the 
computer waits for another set of numbers. If another set is not read in prior 
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Figure 3.5 Block diagram of the revised receiver used to operate with a PDP-15 computer. 
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Figure 3.6 The RF amplifier module for the receiver. 




Figure 3.7 The IF and DC amplifier module. 
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Figure 3.8 The encode pulses as set up by Birley and Sechrist [1971] used to collect data during 
the eclipse, and the revised encode pulses used by the present programs. 
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Figure 3.9 The encode pulse circuitry used to produce the former and present encode pulses 
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to the 150 usee, the set was from an extradordinary radio wave and is rejected. 
Othewise, both sets are accepted and the computer is synchronized with the pulser. 
This process is done only when the computer has a possibility of being out of 
synchronization with the pulser which are: 

1) Beginning of every file 

2) After the transfer of a block of data to disk 

3) After collection is stopped and restarted by console 
control switch 

4) During a timing error (no longer a terminal error, see 
Section 3.3) 

5) When the computer "forgets to read" (discussed in Section 3.3) 

The second change is to account for the nonlinearity of the receiver as 
seen in Figure 3.10 and was initially set up to adjust the data during processing 
[Wiersma and Sedhrist t 1972] . Due to the time needed for the calibrating opera- 
tion (approximately a half day), the computer is used which increases the speed 
of the process while making it possible to account for inaccuracies in the 
analog to digital converter. This process takes about 40 minutes (including 
30 minutes for the receiver warm up). The adjustment to the data is done by 

using a table look-up method in the collection programs. Since the data stored 
on the disk are linearized data, the table is not needed after the collection is 
done and can be deleted after all the data are stored. The method is to convert 
the A/D converter output to the corresponding normalized receiver input. This 
is done by injecting a CW signal of a known value using an attenuator with one 
dB increments and storing the output in the computer using the set up shown in 
Figure 3.11. Straight line segment approximations to the curve in Figure 3.10 
are obtained as shown in Table 3.1. Using outputs from 0 to 511 the corresponding 
inputs are determined normalized to 511 maximum, stored in a table as shown in 




Figure 3.10 Graphs of the input versus output of the receiver used for eclipse data collection 
(old receiver) and the receiver presently being used. The input and output values 
have been normalized to the maximum of the A/D converter (511) . 
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Table 3.1 

Straight line segment approximation to the relationship of 
receiver input to receiver output. 


Slope 

Input 

Output 

Attenuation 

Used 

S< 0=50.906 

TU( 1) = 

0.000 

TUCK !)= 4.786 

99DB 

S< 2)= 3.912 

TUC 2)= 

5.620 

TUOC 2 ) = 4.896 

45DB 

S( 3)= 2,394 

TU ( 3) = 

6,310 

TUCK 3) = 5.073 

44DB 

S( 4) = 7.762 

TU( 4)= 

7.080 

TUOC 4) = 5.394 

43 DB 

S< 5) = 9.151 

TU ( 5) = 

7.940 

TUOC 5) = 5.505 

42DB 

S( 6)= 2.576 

TU( 6) = 

8,91 0 

TUOC 6) = 5.611 

41DB 

S( 7) = 4.190 

TU( 7) = 

1 0.000 

TUOC 7) = 6.034 

40DB 

S( 8) = 1 .983 

TU( 8) = 

11.220 

TUOC 8) = 6.326 

39DB 

S( 9) = 2.255 

TU ( 9)= 

12.590 

TUOC 9) = 7.016 

38DB 

S(10)= 1.676 

TU ( 1 0) = 

14.130 

TU 0 ( 1 0 ) = 7.699 

37DB 

S( 1 1 >= 2.576 

TU ( 1 1 )= 

15.850 

TU0C1 1 )= 8.726 

36DB 

S(!2) = 1.946 

TU ( 1 2) = 

17.780 

TUO( 12) = 9.475 

35DB 

S(13>= 1.356 

TUC 1 3) = 

19.950 

TUO (1 3) = 10.590 

34D3 

S ( 1 4) = 1.754 

TU ( 1 4) = 

22.390 

TU 0 ( 1 4 ) = 12.389 

33DB 

S(15)= 1.381 

TU C 1 5 ) = 

25.120 

TUOCl 5) = 13.945 

32DB 

S(I6) = 1.012 

TU ( 1 6) = 

28.180 

TU0C1 6)= 16.161 

31DB 

S( 1 7) = 1.161 

TU(17)= 

31 .620 

TUOCl 7) = 19.560 

30DB 

S ( 1 8) = 0.969 

TU ( 1 8) = 

35.480 

TUOCl R)= 22.886 

29DB 

S(19)= 0.910 

TU ( 1 9) = 

39.810 

TUOCl 9) = 27.353 

28DB 

S(20)= 1.113 

TU (20) = 

44,670 

TUO (20) = 32.694 

27DB 

S(2l)= 0.987 

TU ( 2 1 ) = 

50.120 

TUO (2 1 ) = 37.591 

26DB 

S(22)= 0.771 

TU (22) = 

56.240 

TUO (22)= 43.794 

25DB 

S (23 > = 1.172 

TU(23)s 

63.100 

TU0C23) = 52.693 

24DB 

S<24)= 0.708 

TU(24)= 

70.803 

TU0C24) = 59.260 

23DB 

S C25 ) = 0.880 

TU (25) = 

79.430 

TU0C25) = 71.451 

22DB 

SC26) = 0.832 

TU (26) = 

89.130 

TU0(26)= 82.478 

21DB 

S (27) r 0.980 

TU(27)= 

100.000 

TUO (27) = 95.540 

20DB 

S(28)= 0.789 

TU (28) = 

112.200 

TU0C28) = 107.992 

1 9DB 

S (29) = 0.900 

TU(29)= 

125.900 

TUO (29) = 125.364 

I8DB 

S (30) = 0.816 

TU (30) = 

141,250 

TUO(30)= 142.426 

17DB 

S ( 3 1 )= 0.918 

TU(31 ) = 

158.490 

TU0C3I ) = 163.559 

16DB 

S (32) = 0.766 

TU(32)= 

177.830 

TU0C32) = 184.627 

I5DB 

S(33)= 0.944 

TUC 33)= 

199.530 

TU0(33)= 212.947 

14DB 

S (34) = 0.913 

TU (34) = 

223.870 

TUO (34) = 238.730 

I3DB 

S (35) = 0.933 

TU(35)= 

251 .190 

TU0(35)= 268.645 

12DB 

S (36) = 1 .030 

TU (36) = 

281 .840 

TU0C36) = 301 .499 

1 1DB 

S (37) = 1.111 

TU(37)= 

31 6,230 

TU0C37) = 334.895 

10OB 

S (38) = 1.260 

TU (38)= 

354,820 

TUO (38)= 369.618 

9DB 

S (39) = 1.306 

TU(39)= 

398.1 10 

TU0C39) = 403.979 

8DB 

S(40)= 1 .471 

TU(40)= 

446.680 

TUO(40)= 441 .155 

7DB 

S(4 I ) = 1.809 

TU ( 4 1 ) = 

501 .190 

TU0C41 ) = 478.202 

6DB 


TU(42)= 

562.340 

TU0(42)= 512.000 

5DB 
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Table 3.2, and placed on a storage device (normally a disk). The program DLOGF 
(given in the Appendix in MACRO language) reads Table 3.2 into the computer, and 
the table is used during collection of the received partial- reflection signal. 
Using the table, the MACRO subroutine LIN does the linearization of the numbers 
read from the analog to digital converter. The programs responsible for the 
formation of these two tables are TBFORL (FORTRAN IV), LINAP (FORTRAN IV), RADC 
(MACRO) , and TTM (MACRO) . 

The system as it has been described was used to collect and process the 
partial-reflection data for the three-day eclipse period of July 9, 10, and 11, 
1972. The rest of this chapter will describe further changes and developments of 
the system. These changes have been due to an increase of 16 K core memory, the 
addition of 2 disk units capable of storing 262.144 words each, and the changing 
from a single user monitor system to a background/foreground monitor disk system. 
3.3 Real-Time Data Storage and Automatic .Processing 

A computer operates on its own timing system and if this timing system oper- 
ates along with events outside the computer that affect the operation of the 
computer, then the computer is said to be operating in real time. For instance, 
if the computer reads in a set of 26 samples and is able to manipulate or process 
them before the next set of samples is read in, the computer is doing real-time 
processing; as opposed to saving the data on tape and processing it later, as done 
by Reynolds and Sechrist [1970]. With high-speed access on the disk (16 msec 
access time), the background/foreground system made possible real-time collection 
and processing of partial-reflection data. Due to the complicated timing, slow 
print-out, and the noise algorithm (discussed in Section 3.4), processing of the 
data is postponed until after the file is stored on the disk. 

The background/foreground monitor system is a double monitor, multi- 
priority level, software system. The two monitors are separate software systems 
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Table 3.2 

The output of the A/D converter aTe numbers between 1 and 511. The 
input for each output is given in this table. 


1 

1 

1 

1 

1 

7 

11 

12 

14 

16 

1R 

19 

21 

22 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

35 

36 

37 

38 

39 

40 

41 

41 

42 

44 

45 

46 

46 

47 

48 

49 

50 

51 

52 

52 

53 

54 

55 

55 

56 

57 

57 

58 

59 

60 

61 

63 

64 

65 

65 

66 

67 

67 

68 

68 

69 

70 

70 

71 

72 

72 

73 

74 

75 

76 

76 

77 

78 

79 

80 

80 

> 81 

82 

83 

83 

84 

85 

86 

86 

87 

88 

89 

89 

90 

91 

92 

93 

94 

95 

95 

96 

97 

98 

99 

100 

101 

102 

103 

103 

104 

105 

105 

106 

107 

108 

108 

109 

110 

110 

111 

112 

113 

1 13 

114 

115 

116 

116 

117 

118 

119 

120 

120 

121 

122 

123 

124 

125 

125 

126 

12 7 

128 

129 

129 

130 

131 

132 

132 

133 

134 

135 

135 

136 

137 

138 

138 

139 

140 

141 

141 

142 

143 

143 

144 

145 

146 

147 

140 

148 

149 

150 

151 

152 

153 

153 

154 

155 

156 

157 

158 

158 

159 

160 

161 

162 

162 

163 

164 

165 

165 

166 

167 

167 

168 

169 

169 

170 

171 

172 

172 

173 

174 

174 

175 

176 

176 

177 

178 

179 

179 

180 

181 

181 

182 

183 

184 

185 

186 

186 

187 

188 

189 

190 

191 

192 

192 

193 

194 

195 

196 

197 

198 

198 

199 

200 

201 

202 

203 

204 

204 

205 

206 

207 

208 

209 

209 

210 

211 

212 

213 

214 

214 

215 

216 

217 

218 

219 

219 

220 

221 

222 

223 

224 

224 

225 

226 

227 

228 

229 

229 

230 

231 

232 

233 

234 

234 

235 

236 

237 

238 

239 

240 

240 

241 

242 

243 

244 

245 

245 

246 

247 

248 

249 

25 0 

251 

251 

252 

253 

254 

255 

256 

25 7 

257 

258 

259 

260 

261 

262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

271 

272 

273 

274 

275 

276 

277 

278 

279 

280 

281 

282 

283 

284 

285 

286 

286 

287 

288 

289 

290 

291 

292 

293 

294 

295 

296 1 

297 

299 

300 

301 

302 

303 

304 

305 

306 

307 

308 

309 

310 

311 

312 

313 

314 

315 

316 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 

328 

330 

331 

332 

333 

334 

335 

336 

338 

339 

340 

341 

342 

343 

344 

346 

347 

348 

349 

350 

351 

353 

354 

355 

356 

357 

35 8 

359 

361 

362 

363 

364 

365 

366 

368 

369 

370 

371 

372 

374 

375 

376 

377 

378 

379 

381 

382 

383 

384 

385 

387 

388 

389 

390 

391 

393 

394 

395 

396 

397 

399 

400 

401 

402 

403 

404 

406 

407 

408 

410 

411 

412 

414 

415 

416 

418 

419 

420 

422 

423 

424 

426 

427 

42 8 

430 

431 

432 

434 

435 

436 

438 

439 

440 

442 

443 

444 

446 

44 7 

448 

450 

451 

452 

454 

455 

456 

458 

460 

461 

463 

465 

466 

468 

470 

471 

473 

474 

476 

478 

479 

481 

483 

484 

486 

488 

489 

491 

493 

509 

494 
51 1 

496 

498 

499 

501 

502 

504 

506 

507 
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sharing the same hardware with programs operating in the foreground system 
having priority. Each system has 8 automatic priority (API) levels and 
a mainstream level. There are four hardware levels which have highest priority. 
The software levels are labeled 4, 5, 6, 7, and 0 where 4 is the highest and 0 
is the mainstream, the lowest. When a program initially starts running in 
either background or foreground, it begins on mainstream. Certain commands 
require a special subroutine called a real-time subroutine and is designated 
a priority level from 0 to 4 and stops all operation on lower priority levels 
(background is lower than foreground) until it exists from the level or per- 
forms an I/O operation. 

With this system the partial-reflection collection and processing programs 
as mentioned could operate in real time, but due to several problems in the 
processing of data, the data could not easily be saved except in 
processed form. The solution used is to collect one file of data and process 
that file while the next file of data is being collected. After each file is 
collected, the operator is told what the next attenuator setting is. The 
collection program also checks the setting of the switches, on the console to 
allow the operator to control parts of the collection program. Switch 0 acts 
as an on/off switch which causes collection to stop collecting and wait in a 
loop if set to 1. Switch 1 allows the background system to share the collec- 
tion and processing storage device (1 disk) if the switch is set to 1. This 
sharing is necessary if the collected files are to be stored on DECtape, 

Switch 3 allows the processed data which are printed out onto the teletype to also 
be punched onto paper tape if the switch is set to 0. This option is presently 
used to allow for later plotting of the data using a programmable calculator. 
Switches 2 and 5 are not used at the present. The rest of the switches are used 
for determining the length of each file (default length is 513 pairs of sets of 
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26 numbers) . The time of day is determined by using the clock within the computer 
to give the time in hours and minutes. 

The flow diagram of the programs is shown in Figure 3.12. The programs are 

loaded into the computer and the computer's clock is set to the time of day. The 
operator is given the option of calibration of the receiver. The linearization 
table is stored on a disk and some initial information is read in. If the table 
read in is erroneous the operator must re-do the calibration procedures. The 
collection is started on priority level 6 and processing waits for the first file 
to be collected. After collection of the number of sets of samples set on the 
console switches and the operator changes the attenuator setting, the second file 
is collected while the first is processed and printed out. This process continues 
until stopped by the operator. Information used to calculate the noise threshold 
as described in Section 3.4 is transferred to the processing program after each 
file is collected and is not stored on the disk. The processing program therefore 
must remain faster than the collection or this information will be lost. 

The processing of files involves rejecting sets of samples that are too 
noisy (discussed in Section 3.4), summing the squares of unsaturated data, sub- 
tracting off the sum of the squared acceptable data, and taking the square root. 
The resulting data are two sets of 21 samples, one of ordinary modes (A^) and 
one of extraordinary mode (A^J radio waves. This process is done in the main 
processing program PROC (given in the Appendix) . The electron densities are 
calculated in CALC2 which is discussed in Section 3.5. The results are typed 
out on the teletype in tabular form as shown in Table 3.3. 

The first line of the print-out of processed data is the heading. This 
gives the time the collection of the file stopped, the date, the reason for the 
run, and the attenuator setting for the file. The next line contains the 
noise threshold and the square of the multiplying constant used in the 
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Figure 3.12 A diagram of the control flow of the partial-reflection programs. The programs operate 
on the API level of the preceding program unless otherwise stated. The collection and 
processing program operated in parallel with the collection programs operating on API 
levels 5 and 6, whereas everything else operates serially. 




















TABLE 3.3 
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Results of CALC2 


1 215 

8-14-73 

DAILY 

RUN 


10DB 

MAX. ALLOW. NOISE® 

19.8 

MULT. 

CONST.® 

9.610 

0 -NOISE 
X-NOISE 

AV.U) 
AV. < 1 > 

19.9 

12.3 

C 2 ) 
C 2 ) 

11.6 

7.5 


513 SAMPLES 

44 REJ 

. (NOISE) 



iiE J . 

C N. + SAT. 

> HEIlSHT 

AV. AO AV. AX AX/AO 

ED 

44 

60.3 

6.9 

4.1 
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44 

61.5 

9.9 

2.3 
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44 

63 » 3 

3.3 

-2.1 

0.30 

0 • 

44 

64.5 

3.8 
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44 
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44 
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22.5 

34.9 
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33.6 

44.9 

1.47 

199. 
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33.8 

43 . 0 
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44 

73.5 

36.2 

39.0 

1 .08 

339. 

44 
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54.2 

49.8 
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298. 

44 

7 6.5 

59.1 

45. 1 

3.76 

343. 

44 

78.0 

58*9 

34.1 

0.58 

53 2. 

44 

79.5 

57.5 

27.3 

0.48 

416. 

44 

81 .0 

56*5 

23.2 

0.36 

682. 

44 

82.5 

52.7 
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44 

84.3 

46.6 
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0. 23 

623. 

44 

85.5 

99.0 
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3.12 
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44 
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maximum noise criterion discussed in Section 3.4. The next two lines are the 
ordinary and extraordinary mode noise before (number 1) and after (number 2) re- 
jections due to excessive noise. The next line gives the number of pairs of sets 
of 26 samples collected and the number of these pairs rejected due to saturation. 
The first column of the table is the number of rejections due to both saturation 
and excessive noise for each height. The next column gives the height of the 
reflected signals for each row. The next two columns give RMS of the ordinary 
(Aq) and extraordinary signals. The fifth column gives the ratios of extra- 

ordinary partial reflections to ordinary partial reflections from the fourth 
and third column respectively. The last column gives the electron density for 
between the heights. The last electron density is given as zero since only one 
height is available to calculate it. 

The present method of collection and processing of partial-reflection data 
is fast, efficient, and easy to operate, but two problems needed to be removed. 

The increase of input/output operations have increased timing errors which are 
discussed by Birley and SeahHst [1971], and the A/D converter sometimes fails to 
respond to read commands . 

The A/D converter transfers data to the computer using multicycle block 
transfer as described by Birley and Sechrist [1971]. The process is a three 
cycle operation for each word transferred. After each transfer, the A/D con- 
verter interface is tested for synchronization. If the timing between the inter- 
face and the I/O processor is altered, transfer is stopped resulting in a timing 
error. With the present system, this error can result from hardware malfunction 
or excessive I/O operation occurring. If the latter is the reason, the problem 
is only temporary and can be remedied by issuing another read. Care is taken to 
keep the computer synchronized with the pulser. If the error is a hardware 
problem, the condition will not clear up and collection must stop. The error 
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will usually occur when data are being collected, processed data are being printed 
out, and a tape is being copied onto the disk in background, all simultaneously. 

The second problem has to do with the A/D converter's interface refusal to 
transmit data. The problem has been traced to failure in the A/D converter inter- 
face logic. The collection program will issue an A/D converter read, but not 
receive control back and no data are transferred. This problem occurs only with 
the background/ foreground system and it occurs infrequently (once in about every 
10,000 read commands). One solution is to issue a double read, but the problem 
could still occur. The solution used is for the processing to check for this 
stoppage, restart the collection in an orderly fashion if it has stopped and to 
ring the teletype bell to let the operator know of the stoppage. This solution 
does not prevent the failure of the A/D converter interface to transfer data, and 
the problem will have to be removed for faster ratio of collection, but presently 
the operator need not be concerned with this problem. The rest of the 
data is unaltered by this problem. 

3.4 Noise Bejection 

The partially reflected radio waves from the D region are usually small in 
amplitude on the order of 10 to 1000 mvolts at the output of the 80 dB gain re- 
ceiver. Noise amplitudes vary between 30 to 1000 mvolts. For the purpose of the 
noise algorithm, noise is considered to be any interference which is part of the 
receiver output signal that is not attributed to the partially reflected waves 
from the vertically transmitted pulse. This noise is divided into two types: 
background noise and noise bursts. Background noise is noise caused by the re- 
ciever (14±3 mV) and general atmospheric noise which is always present (40±10 mV). 
Noise bursts are caused by lightning and other radio transmitters, and the ampli- 
tude of this noise is dependent on the location of the source. Lightning noise 
will usually last for the duration of one encode pulse while noise due to other 
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transmitters will last for at least 1/2 second which is several encode pulses 
(see Figure 3.8) and the noise will be increased usually by 10 to 1000 mvolts. 

Both types of noise are rejected in the processing program PROC (FORTRAN IV) as 
shown in the block diagram of this program in Figure 3.13. 

Data are collected in pairs of sets of 26 numbers. Each set contains 5 
noise samples and 21 samples of partially reflected signals. Each pair contains 
a set of ordinary mode samples and a set of extraordinary mode samples. In PROC 
a noise threshold is determined and the square of this multiplied by five is 
compared to the sum of the squares of the five noise samples of each set. This 
method of comparison is faster than comparing the RMS of the noise as set up by 
Birley and Seahvist [1971] since square root operations take approximately 1 
msec and squaring takes 70 usee on the PDP-15, and the squaring need only be 
done once per file. If the noise of either mode is greater than the noise 
threshold, both sets of 21 signal samples are rejected and the next pair of sets 
are tested. If the noise of both modes is less than this threshold, the noise 
of both sets are considered acceptable and saved for later processing. The 
partially reflected signals with acceptable noise for each mode are checked for 
A/D converter saturation (.997 volts receiver output) at each height. If either 
of the two samples (one of each mode) is saturated at a height the two samples are 
rejected; otherwise the data are considered acceptable. This processing of 
pairs of 26 samples continues until the end of the file is reached. After the 
file of collected data has gone through this processing, the average of the sum 
of the squared acceptable noise for each mode is subtracted from the average 
of the sum of the squared acceptable partially reflected samples of the same mode 
at each height, and the square roots are printed out as shown in Table 3.3 and 
as described in Section 3.3. 
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Figure 3.13 A flow chart of the processing program PROC 
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Originally, the noise threshold was determined by the operator typing in 
a value chosen by him as seen in the program PROC73 in the Appendix. This was 
later changed to an automatic determination based on the attenuator setting 
used as given at the beginning of a run. This method did not account for the 
day-to-day variation in noise nor in an erroneous attenuator setting. The 
noise threshold value is presently determined by the following equation: 

45 

M = ( * ( £ AO/45) 2 (3.1) 

where M = maximum allowable noise value 
K = arbitrary constant 

N = certain noise samples collected as explained in the following 
paragraph . 

In the collection programs RSUB and LIN, the maximum and sum of each 
group of 45 noise samples are stored, and the maximum values are compared. 

The sum of the group with the lowest maximum value is transferred to the 
processing program PROC and is used in equation (3.1). The constant K has 
been chosen by trial and error, and values between 2.5 and 3.5 seem to give 
the best results (equation 3.1 is being used). 

Other algorithms have been tried, but none seem to give any obvious 
improvement in the resulting electron densities. One method is to split 5 
noise samples collected with each set of data into 2 for comparison with the 
noise threshold value and 3 subtracted from the reflected signals. This method 
works on the theory that the noise within the 5 noise samples is not the same 
amplitude as the noise within the 21 data samples for each set of 26 data 
samples, but is statistically the same over the number of samples collected 
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for one file. With the present system, when the number of rejections due to 
noise is large, (greater than 200 out of 513 pairs of sets of samples), the 
noise within the noise gate is restricted to a lower level than the noise in 
the data frame. Therefore, the noise in the data frame would not be completely 
subtracted off; as it would be with splitting the noise samples. The applica- 
tion of this technique using 4 noise samples showed no improvement in the 

results. Two possible causes are too few noise samples being used and the noise 
samples being too close together. 

Another method has been developed and tested by D. R. Ward [private com- 
munication] . A CW signal is inputed into the receiver along with the received 
data from the antenna. The noise and partially reflected signals are each 
defined as A cos0; where A is the amplitude and 9 is the phase. The noise is 
assumed to be random while the partially reflected signals are assumed to have 
only a small variation between two sets of samples. Using an algorithm developed 
by D. R. Ward [private communication], the phase and the amplitude of the noise 
portion of each signal average to zero while the phase and amplitude of the 
signals do not. This method is used to reject the noise from the partially 
reflected signals at each height. This method fails to reject interference 
caused by other transmitted signals since this type of noise does not have 
random phase. D. E . Ward [private communication] has obtained useful electron- 
density profiles from the method but generally found no improvement over the 
present system. Further study and development of either method may improve the 
processing and should not be discarded. 

3.5 Converting A x /A q Ratios to Eleatron-Density Profiles 

The partial-reflection programs assume a constant collision frequency for 
e<±L.h height with seasonal variation. The values used were determined from the 



49 


following equation [Birley and Seahrist, 1971]: 


= *P (3.2) 

where K = constant - 7.3 x 10 5 
p = pressure in pascals 

= collision frequency in sec” 1 

The pressures used are from the mean atmospheric model from COSPAR International 
Reference Atmosphere (1965) with seasonal variations given by U. S. Standard 
Atmospheric Supplements (1966). Using these pressures, experimentally the values 
calculated for K vary by as much as 2 x 10 5 [ Lodato and Meahtly 3 1971] . The 
seasonal variations in the collision frequency (Figure 3.14) can vary by as much 
as 20%. This 20% variation in can cause the calculated [e] to vary by a 
factor of 1.2. The electron densities are calculated using the refractive index 
equation given by Sen and Wyller [1960] and several approximations as discussed 
by Pimat and Bcndhill [1968]. The resulting equation given by Reynolds and 
Seahrist [1970] is: 
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e = mV 2 /2kT 
[e] = electron density 
e = electron charge = 1.6 x 10“ 19 C 
m = electron mass = 9.1 x 10 ”^ 1 kg 
e 0 = permittivity of free space = 8.85 x 10 -12 F m -1 
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<*) = angular frequency of the transmitted wave 
0 )^ = gyro -frequency of the electron 
= lower height 
h ^ = higher height 
Ah = 7? 2 - 

k = Boltzmann constant * 1.38 x 10 -2 ^ J® K -1 

T = temperature 
V - electron velocity 
R Q = ordinary mode reflection coefficient 
R = extraordinary mode reflection coefficient 

This equation required a set of collision frequency constants which are given 
in the program CALC (FORTRAN IV). The ratio (/?//?), / {R /R ) and FD 

CC O rl 2 * CC O rl 

(equation (3.4)) are calculated in ELDEN (FORTRAN IV). CAL 2 (called by PROC) 
uses these values (which vary only with v^) as constants for each pair of 
heights to calculate the electron densities according to Equation (3.5) 

[e] - ln(RATI02 x (yy ^ C 3 - 3 ) 

where RATI02 = Cys,k / CV*A • 

2 1 

This method is used to reduce the amount of core memory required and in- 
crease speed of execution of the program. A new CALC2 can be obtained by revis- 
ing the collision frequencies and running the program CALC which writes the pro- 
gram CALC2. The electron densities are printed out as shown in Table 3.3 and 
described in Section 3.3 
3.6 Equipment Testing 

The equipment needs to be tested periodically to determine if it is in 
operating order. The transmitter is tested by observing and keeping a log of 
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the voltage and current at various locations via meters and an oscilloscope. 

The antennas are tested by transmitting and receiving signals at various times 
during the day. At noon the extraordinary signal should be absorbed and at 
night the ordinary signal should be absorbed. By transmitting and receiving 
ordinary and extraordinary signals as described in Progress Report 73-1 
[Edwards , 1973], the phase and attenuation of each antenna of each array can 
be set and checked for possible damage. This process is also a partial check 
for the transmitter and receiver. A spot check of 30 dB difference in ordinary 
and extraordinary reflections from the E region at noon is done on a daily 
basis. 

The program CHECK (FORTRAN IV) has proved valuable in checking the receiver 
and the analog to digital converter. CHECK performs a modified dump of the A/D 
converter as read by the computer. If the number 31 is typed, the output is in 
the form of partial-reflection data (ordinary and extraordinary pairs) , patterned 
after the new encode pulse shown in Figure 3. 8. if any other number is typed in 
an average of that number rounded to the next higher multiple of 50 is printed 
out. The 31 pairs of samples are printed out in millivolts only, while the 
averages are printed out in millivolts and as represented in the A/D converter. 
This program has had many applications; it showed the blanking gate oh a new 
receiver to be too long. It was used to calibrate the A/D converter using an 
input from a standard source. Table 3.4 shows the accuracy of the A/D converter 
as the standard voltage source was varied from 1.0 volts to .1 in .1, .01, and 
.001 volt increments. It was used in comparing the paper punch system set up 
by Reynolds and Seahrist [1970] with the computer storing method presented in 
Section 3.3. CHECK has also been used to determine the number of samples 
required to have less than 10% error due to noise (at least 100 samples are 
required) . The program is easy to operate and has become important in testing 
and checking the receiver and the analog to digital converter. 
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Table 3.4 

The output of the A/D converter using a 
calibrated input source 




ADC 



Output 


Input 

Average 


Voltage 

Voltage 

511.204 


998.444 mV 

1000 mV 

460.558 


899.527 mV 

900 mV 

409.625 


800.050 mV 

800 mV 

358.528 


700.250 mV 

700 mV 

307.057 


599.720 mV 

600 mV 

256.020 


500.040 mV 

500 mV 

205.252 


400.883 mV 

400 mV 

154.082 


300.941 mV 

300 mV 

102.787 


200.756 mV 

200 mV 

51.076 


99.758 mV 

100 mV 

46.349 


90.525 mV 

90 mV 

40.843 


79.772 mV 

80 mV 

35.208 . 


68.766 mV 

70 mV 

30.844 


60.242 mV 

60 mV 

25.769 


50.330 mV 

50 mV 

19.976 


39.016 mV 

40 mV 

15.022 


29.339 mV 

30 mV 

10.200 


19.922 mV 

20 mV 

4.830 


9.433 mV 

10 mV 

3.857 


7.533 mV 

9 mV 

3.233 


6.315 mV 

8 mV 

3.010 


5.880 mV 

7 mV 

2.847 


5.560 mV 

6 mV 

2.443 


4.771 mV 

5 mV 

2.054 


4.012 mV 

4 mV 

1.404 


2.743 mV 

3 mV 

0.659 


1.286 mV 

2 mV 

0.125 


0.244 mV 

1 mV 

0.010 


0.020 mV 

0 mV 
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3.7 Future Development 

Several improvements are being made to the system. A new receiver is being 
made using a linear detector and new RF and IF stages to reduce the receiver 
noise. Figure 3.10 shows a comparison of the input versus output between the 
new receiver and the old one. With no input signal, the noise level of the new 
receiver is 2.5 mV and the level of the older receiver is 14 mV. The circuitry 
and discussion of it are given in the Aeronomy Progress Report 73-1 [Edwards , 
1973]. 

A digital input/output device is presently being sought which would improve 
the calibration time and free the operator for other tasks as well as simplify 
the operation of the system. The purchase of such a device would also reduce 
the amount of paper presently required. 

Another asset would be a line printer. One could reduce the processing time 
by at least half and allow for more sophisticated processing (with possibly 
better noise rejection) if such a line printer were purchased. 

As mentioned by Birley and Seohrist [1971], an increase of transmitter 
power is also needed. This would improve the signal-to-noise ratio and give 
better data below 70 km. 

The noise problem should be studied more carefully. Perhaps a combination 
of the method discussed in Section 3.4 would improve the results. Another 
possibility would be to reject extremely low values of reflected signals. 

An additional program to transfer collected data to tape would be helpful. 
The original programs set up by Birley and Seohrist [1971] saved data on tape 
for future processing. With the present system, collected data can be stored on 
tape by using a system program called PIP. This requires knowledge in operation 
of the computer, and the transferring of files can get complicated. 



4. EXPERIMENTAL RESULTS 


This chapter describes the results from partial- reflection data which was 
collected and processed by the computer on July 9, 10, and 11, 1972. A solar 
eclipse occurred on July 10, 1972. The obscuration function shown in Figure 1.1 
shows the first contact to be at 1319 CST and the last contact to be at 1536 CST 
with 60% of the solar disk obscured. The data were collected from 1200 to 1700 
CST to show the effects of the solar eclipse on the electron density and 
collected between the same times on July 9 and 11 to be used as control data. 
Data were collected in blocks called files. Each file of data, consisting of 
1026 sets of 26 numbers, was collected and stored on DECtape every 3.8 minutes. 
The signal prior to entering the receiver was attenuated with four attenuator 
settings (0, 10, 20, and 30 dB) . Each file was collected beginning with 

the lowest attenuator setting of 0 dB with each subsequent file collected at 
the next attenuator settings; 10, 20, and 30 dB, respectively. This process 
was then repeated. This process was used to obtain the very small echoes as 
well as the very large ones. The files of data are divided into approximately 
15 minute intervals, corresponding to the four attenuator settings. 

The data between 1400 and 1430 on July 9 was lost due to an erasure of 
the disk before it could be processed. These data have been interpolated. The 
data from July 10 between 1200 and 1300 was erroneous and therefore has been 
eliminated from the results. The computer results were processed further 
combining the files with different attenuator settings. 

4 . 1 Reduction of Data 

Individual results shown in Figure 4.1 show valid electron densities but 
are limited height range; therefore, multiple attenuator settings were used to 
obtain usable data over a greater range of heights. The computer processes 
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JULY 11,1972 
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Figure 4.1 Comparison of electron-density profiles on July 10 and 11, 1972. 

The data were taken at 1432 CST with the attenuator set at 30 dB 
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only one file at a time; therefore, further processing was necessary to com- 
bine four files corresponding to the four attenuator settings into one set of 
results. Three methods have been developed to accomplish this. The first 
method was originally used but problems developed in determining acceptable 
data and method two was used. Using method two, some acceptable data were 
being ignored and the 20 and 30 dB settings were found to give similar results. 
Therefore, method three was developed to utilize much of this acceptable data 
that were being ignored. 

1. In method one, the results with the lowest attenuator setting (0 dB) 
were used for 60 km up to the height where 5% of the ordinary and extra- 
ordinary data was rejected due to saturation (see Section 3.4). The 
electron densities for the higher heights were obtained from the next 
higher attenuator setting under the same restrictions of saturations. 

This process continued until the last electron density was obtained. 

The results of this method seemed to be satisfactory except for above 

81 km and below 66 km. 

2. Method two is the same as method one, but accounts for inaccuracies in 
the receiver by rejecting electron densities that used A S /^ Q ratios 
that were less than .09. Electron densities were rejected also if the 
signal to noise ratio was less than 1. These two revisions eliminated 
much of the results below 65 km and above 85 km. 

3. Method three is similar to method two except for the way the multiple 
attenuators are combined. The electron densities are considered accept- 
able if the A 3 /A q ratios for both heights are greater than .08, the signal 
to noise ratio is above 1 for both heights, and the rejections due to 
saturations were less than 5% for both heights used to calculate the 
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electron density. If more than one attenuator setting had acceptable elec- 
tron densities for between two heights, then the median of the acceptable 
electron densities was used. Using these three methods, the computer re- 
sults were combined to give one electron^density profile for every 4 attenua- 
tor settings. Using either average or medians, electron densities of dif- 
ferent heights or of different times were combined as discussed in 
Section 4.2. 

4.2 Electron-Density Results 

The results are presented in two forms: by the total differential absorp- 

tion below each height (A 3 /A q ratios) and by electron densities. The A JA 
ratios given in Figures 4.2, 4.3, and 4.4 are plotted using a sixth order poly- 
nomial approximation of the ratio as calculated by method one. The eclipse 
shows a reduction in absorption which indicated a reduction in electron density 
as expected. The third day shows irregular absorption with a large increase 
in absorption. Referring to Figure 2.2. the increase in absorption is related 
to the X-ray flux burst. The electron density for above 75 km for the three 
days given in Figure 4.5 shows a good correlation between the large increase 
in electron density on July 11 and the burst of X-ray flux. Due to this obvi- 
ous contamination, the second control day is not used for comparison during the 
burst period. 

Figure 4.6 gives the AJA q ratios versus height. The ratios were deter- 
mined using method one and taking the median of the groups within the hour 
corresponding to the maximum obscuration of the solar eclipse (1400-1500 CST) . 

Due to the much larger absorption in the control days than during the eclipse, 
the electron densities above 81 km (approximately) are not valid according to 

method two and three, but with the eclipse day, the values should be acceptable 
up to 85 km. 
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Figure 4.3 Comparison of the A I A ratio at 75 km for July 9, 10, and 11 
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Figure 4.6 Median A / A profiles between 1400 and 1500 CST for each day. 

CC 0 




64 


Figure 4,7 gives the electron-^ density variation with time. The electron 
densities are averages between 70,5 and 78 km and between 78 and 87.5 km with 
the electron densities obtained by using method one for processing the compu- 
ter result. Figure 4.5 and 4.8 give the electron-density median for 75 to 
82.5 km and 67.5 and 75 km, respectively, as each varies with time. These elec- 
tron densities were obtained using method three. At the lower altitudes, the 
median electron densities show no effect from the eclipse while the average 
electron densities do show a slight effect. This difference, though, is mainly 
attributed to the higher heights the averages were taken from rather than to 
the method used. The highest heights show large effects due to the eclipse. 
Figure 4.7 shows a minimum electron density near maximum obscuration of the 
eclipse while Figure 4.5 shows the minimum being delayed by half an hour. 

This is attributed to the variation in the data due to the inaccuracies in the 
partial-reflection equipment. The X-ray burst shown in Figure 2.2 seems to 
have no effect at the lower altitudes. 

Median electron-density variations with height are given in Figure 4.9. 
These values are the median obtained by processing the computer results util- 
izing method two and finding the median value between 1400 and 1500 CST. 

Below 75 km the eclipse does not seem to have much effect on the electron 
density as shown in Figure 4,9, but above 75 km, the electron density decreases 
by 45 to 65%. The upper height for this comparison is 81 km due to the small 
A a/ A o ratios (shown in Figure 4.6). The electron-density profile shows some 
conformity to the expectation given in Section 2.4. 

4.3 Theoretical Applications 

Since the eclipse never reached totality, the electron production (q) can- 
not be assumed to be zero, but equation (2.4) can be used as an approximation 



ELECTRON DENSITY (cm* 3 ) 


65 


2000 

1000 

500 




1300 1400 1500 1600 1700 1800 

TIME (CDT) 


Figure 4.7 Average electron densities between the altitudes 
78.0 - 82.5 km and 70.5 - 78.0 km. 
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Figure 4.8 Median electron densities between 67.5 and 75 km. 
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Figure 4.9 Median electron-density profiles between 1400 and 1500 CST 



68 


to [e] and Equation (3.1) 

q = (NO) [NO] e" T F q (3.1) 

— 3 1 

where q = electron production rate in cm - sec 

^ (NO) = ionization cross-section of nitric oxide = 2 x 10' 18 cm 2 

[NO] = number density of nitric oxide in cm 

J » = incident lyman-alpha flux at the top of the atmosphere = 3.1 x 10 11 
-2 -1 

photons cm sec 

F = the function of the unobscured solar disk 
t = optical depth 

given by Sechrist [1966], was used to approximate the electron production rate 
and equation (3.2) was used to approximate the optical depth. 

T * a a t( V H secx (3.2) 

[O 2 ] = number density of molecular oxygen in cm - ^ 

11 - scale height 
X = solar zenith angle 

Figure 4.10 shows the variation of q during the eclipse as compared to the 
variation without the eclipse. The electron production rates were used to 
obtain theoretical electron densities with <* e ££ being chosen to give the best 
fit to the experimental results. A value of 2 x 10' 6 for « e£f was determined for 
the eclipse day between 75 and 82.5 km and 1.77 x 10 -6 for the same height range 
on the control days. For the heights 78 to 87.5 km a e££ was found to be 8.46 x 
10 . These values for a g ££ are similar to ones given by Mitra [1968]. Figure 
4.11 shows a comparison between the theoretical [e] during the eclipse and without 
the eclipse using a eff of 1.77 x 10' 6 . 



ELECTRON PRODUCTION RATE (cm -3 sec” 



TIME (CST) 

Figure 4.10 Electron production rate between 75 and 82,5 km during the eclipse 
and during the control days. The NO distribution used is from 
Meira [1971]. 
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The electron density of the eclipse was divided by average electron density 
of the control data and compared to the obscuration function as seen in Figure 
4.12. The comparison of the experimental [e] during the eclipse and the theo- 
retical [e] without the eclipse using equation (2.4) was also made and is shown 
in Figure 4.13. 

The electron density for July 9 shows a good correlation with the solar 
zenith angle (Figure 4.14) and was therefore divided by the theoretical [e] to 
eliminate the effects of the solar zenith angle and to determine the variability 
of the experimental [e] (Figure 4.15). The same comparison is made with the 
eclipse [e] (Figure 4.15) and shows a similar but greater variability. 

Generally, the eclipse electron densities show a decrease that is greater 
than expected from the equation (2.4). Other than the possibility that this is 
caused by variabilities due to inaccuracies in the experiment, there are three 
reasons why this may occur: 

1. The obscuration function of the ionization source (Lyman-a) is 
different than the uniform-disk obscuration function used. 

2. The « e ££ increased during the eclipse. This could be caused by a 
change in the hydrated-ion composition between 75 and 81 km. 

3. Loss by attachment is increased by the eclipse. 

The electron-density profiles in Figure 4.10 show good comparison with the 

profile with 40% obscuration given in Figure 2.7 and with 60% obscuration shown 
in Figure 2.6. Smith t et al . [1965] described small changes below 70 km as the 
C-layer caused by cosmic rays which disappear as the eclipse reaches totality. 
The effect can be seen up to 69 km in Figure 4.9. 

4 . 4 Summary 

Comparing Figure 4.9 Figures 2.6 and 2.7, the electron-density profiles 
of this eclipse are similar to previous eclipses for the same obscuration. 
Generally, similar conclusions can be drawn. The difficulty in interpreting the 
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Figure 4.12 The ratio of electron densities for the average of the control day 
as compared to the unobscured sun. 


to 




TIME (CST) 

Figure 4.13 The graph of the ratios of the theoretical [e] for the unobscured sun to the 

experimental [e] for the eclipse as compared to the unobscured sun. The a e ff 
used for 75 to 82.5 km is 1.77 x 10“^ and for 78 to 87.5 km is 8.46 x 10 - ^. 
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Figure 4.14 Scatter plot correlating the electron density for July 9, 1972 
between 75 and 82.5 km to the solar zenith angle. 
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Figure 4.15 The graph of the ratio of theoretical electron densities 

to the experimental electron densities for July 9, 10, 1972 






results lies in the variation of the electron density of the eclipse with time. 

In Figures 4.12 and 4.13 a small decrease in electron density precedes the obscura- 
tion of the sun. An error of 20% can be expected due to the equipment and 20% 
error can be expected in the collision frequencies. Errors due to collision fre- 
quencies will cancel in Figure 4.12 but the errors due to the equipment will in- 
crease. For Figure 4.13, the reverse is true, but there are also errors due to 
the approximations made in equations (2.4), (3.1), and (3.2). With these pos- 
sibilities of errors and observing that the ratio after the eclipse can get as 
low as .8 in Figure 4.12, the initial decrease can be interpreted as experimen- 
tal error. The errors in Figure 4.13 can be seen in the variations in Figure 4.15. 

No correlation could be seen between the X-ray flux and the electron density 
on the third day except furing the X-ray burst period. Therefore, Lyman-a is 
assumed to be the main ionization source and the theoretical calculations were 
made on that assumption. 

Of the three reasons for the large decrease in [e], the effects due to 
changes in hydrated ions is the most likely. During the day electron loss by 
attachment is insignificant above 75 km. Since the obscuration of the sun was 
only 60% which corresponded to a production rate similar to that of 65° solar 
zenith angle, the loss process would still be by recombination. 

The larger concentrations of Lyman-a on the solar disk were in the southern 

O 

hemisphere and were not obscured and the intensity of 1-8 A X-ray flux was too 
small to have any large effect. Therefore, the obscuration function of the 
ionizing source would have the same obscuration or less. This leaves the only 
possibility for the larger decrease in free electron as being due to changes in 
the hydrated ions. 
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5. CONCLUSIONS 

The solar eclipse provides a good opportunity to study several processes 
of the D region and to develop its theoretical model. Accurate interpretation 
of the eclipse data is required to determine exactly the P-region ion produc- 
tion and loss processes, the variation of formation of hydrated ions, 

and negative ion chemistry. A brief theory of the ^-region chemistry is pre- 
sented in Chapter 2 and used to analyze the data in Chapter 4. The equipment 
used in the collecting and processing of the partially reflected waves, as 
well as the refinements made in the collection process are given in Chapter 3. 

The newer partial-reflection system, discussed in Chapter 3, has been in use 
for the daily collection of data. Results from this newer system are given by 
Denny and Bowhi.ll [1973], This chapter reviews the results of the partial- 
Xeflection data taken during the eclipse and suggests further developments 
of the partial -reflection system. 

5 . 1 Review of Results 

The effect of the eclipse below 75 km is below the experimental errors. 

These errors are due to the variability of receiver gain caused by temperature 
fluctuations, the 40 psec pulse width of the transmitter, inaccuracies in the 
collision frequencies, and inaccuracies in noise reduction. In comparing the 
[e] profiles for July 9 and 10, 1972 in Figure 4.8, the beginning of the forma- 
tion of a C layer can be seen resulting from cosmic rays. From Section 2.1, 
the main ionization source between 70 and 80 km is Lyman -a since the X-ray source 
effects were not observed below 81 km except when the X-ray flux increased above 
1 x 10 erg cm sec 

The decrease between the electron density from July 9 and from July 10 is 
dependent on the height and is very marked between 79 and 81 km. Near 80 km. 
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this change in electron density is as much as 55% between the results of July 9 
and 10, which was not expected according to equation (2.4), The most probably 
answer given in Chapter 4 is that it is due to an initial large decrease in 
hydrated positive ions which are the major ions between 75 and 80 km during the 
daytime (as seen in Figure 5.1 by Krankowsky i et at-. [1972]). 

The theoretical [e] were used to compare with the experimental [e] in 
Figure 4.13 to remove any electron density variability not due to the eclipse. 

The results in Figures 4,12 and 4.13 show unexpected initial decreases in [e] 
prior to the eclipse and larger decreases than would be expected during the 
eclipse, but allowing for 20% error in these results, these variations are 
within the error limits. In general, there is good agreement with the data 
from Smith , et al. 3 [1927] and Deeks [1966]. 

5.2 Suggestions for Further Work 

The present partial -reflection system has proved invaluable in presenting 
variations in electron densities diurnally and from day-to-day as presented by 
Denny and Bowhilt [1973]. The system has several limitations, though. Either 
the signal-to-noise-ratio should be increased or the rates of data collection 
increased. Both of these changes would require alterations in the transmitter. 

By doubling the peak power of the transmitter, meaningful partial reflections 
could be obtained at lower altitudes without excessively disturbing the iono- 
sphere due to the slow pulse rate as is done in the cross modulation experiment. 

By increasing the pulse rate, more data could be collected in the same interval 
of time, allowing for a more accurate statistical evaluation of the noise. 

A new receiver has been built as mentioned in Chapter 3. The initial results 
obtained using it show an improvement in the results, but the problem of elimi- 
nating atmospheric noise remains. The main problem lies in defining the noise. 



ALTITUDE- 



PULSES PER MASS 

Figure 5.1 Rocket measurements of the positive-ion chemistry by 
Krankowsky 3 et at . [1972] . 
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A study should be done on the specific types of noise received and the algorithms 
required to reject each. This would include receiving and storing noise on 
DECtape for later evaluation of the amplitude and phase. 

A digital input/output would increase the efficiency of the collection and 
process. Presently the system requires the assistance of the operator every 
3-1/2 minutes and uses one page of computer paper for every page of data. With 
a digital input/output , the computer could set the attenuators and control other 
switching which would free the operator for other tasks. This would also improve 
the usefulness of taking differential phase measurements as described by Wiersma 
and Seohrist [1972] . 

Using a line printer for outputting the data would allow for more sophis- 
ticated and complicated processing of data. This would also be required if the 
rate of collection is increased. To collect one file of data takes 3.5 minutes, 
to process one, about 45 sec, but to print out the results on the teletype and 
paper tape takes 2.6 minutes. Therefore, the processing would not be able to 
keep up with a faster collection unless the speed of printing the results 


increased. 
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/ ADDRESSES FOR 1,2* 3,4# 5 

/ AND 10 

MINEQ /USED TO DETERMINE LENGTH OF RUN 

/CONVERTS INPUT NUMBERS TO 
/ THE EQUIVALENT 

/ BCD # OF MINUTES 

/USED TO DETERMINE MULTIPLES OF 10 
/TEMP. STORAGE FOR DB SETTING 
/TEMP. STORAGE FOR DB IN ASCI! CODE 
/COUNTER FOR # OF DB TO BE USED 
“2 /ALLOWS NO MORE THAN 2 DIGIT DB’S 

-A /ALLOWS ONLY 4 DB SETTINGS 

/* 55 TO CECK FOR CARRIAGE RETURN 
MINEQ+2 /FOR LESS THAN 10 HOUR RUN 

/ASCII SPACE AND A ZERO 
/CHECK FOR CARRIAGE RETURN 
/CHECKS FOR MULTI -DB SETTINGS 
/ASCII DEFAULT DB FOR 0 DB SETTING 
/USED TO FIND THE DEVICE # IN .IOIN 
/-12 TO FIND BCD 0 MULTIPLE OF 10 
/SETS THE DB SETTING 
/ TO ASCII CODE 
/INITIALIZES .IOIN TABLE POINTER 
/USED TO DETERMINE VALUES OF ASCII #*S 
/LOOKS FOR DEVICE # 5 CDISK) 

/LOOKS FOR THE DEVICE # 4 (DECTAPE) 

/USED TO FIND THE UNIT # FOR A DEVICE 
/CHECKS FOR DECIMAL POINT 

*32 /USED TO CHECK FOR NONNUMBER ASCII CHAR. 

/SETS MAXIMUM AMOUNT COLL. ON A DISK 
DECN /INIT. MULTIPIERS FOR BCD f'S 

H1NEQ /INIT* MULT* TO CONVERT TIME TO MIN 

/SETS UP 0 OF DISKS TO BE USED 

/CHECKS FOB UNIT l 

/CHECKS FOR UNIT 3 

/CHECK FOR ASCII "1” IN DATE 

/USED TO CHECK FOR ASCII M 5" TO "7” 

/USED TO GET SECOND CHAR. IN ASCII WORD 
/CHECKS FOR AN ASCII ZEHO 
/USED TO CHECK FOR AN ASCII M A M 
/USED TO INCREMENT 2 ASCII LETTERS 
/USED TO INCREMENT 4 ASCII LETTERS 
/USED TO INCREMENT 1 ASCII LETTER 
/CHECKS FOR AN ASCII "A" 

/USED TO MASK ALL HIGHER BITS 
/CHECKS f O r FIRST HALF ASCII "Y" 

/DEFAULT VALUE FOR TIME OF DAY <20 Hi.) 
/CHECKS FOR UNITS LESS THAN 4 
/DUMMY CODE FOR .IOIN TABLE 
/SECOND WORD FOR *IOIN TABLE 
/LOOKS FOR DEVICE 07 <P. PUNCH) 

/ USED TO GET CONSOLE SWITCH 1 


* £ jec r 

/S $5$ $*$$$$£ ********* ****$*$$ 

/ BEGINNING OF THE PARTIAL REFLECTION PROGRAMS. 

/ THIS PART INITIALIZES THE .DAT SLOTS, SETS UP NEW 
/DIRECTORIES* ALLOTS SHARE MODE* AND GETS AND INITIATED thf 
/CLOCK ro GIVE THE TIME OF DAY* ™ E 

/ 


SS SSSiSSS 


1NTIM 

* INI T 

• INIT 
LAC* 

I AC 

TTO* OUT* IN TlM 
TTI * IN* IN TIM 

DMV4 

/TTY OUT 
/TTYIN 

/GET ADOrt. OF FOREGND .DAT SLOTS 


DAC 

FDAT1 

/ ONE 
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IAC 

DAC 

IAC 

DAC 

IAC 

IAC 

DAC 

TAD 

DAC 

LAC* 

DAC* 

LAC* 

DAC* 

• INI T 

• CLEAR 
LAV 
TAD 
SPA 
JMP 
LAC 
DAC 

DAC 

■ INI T 

« CLEAR 

JMP 

DZM* 

DZM* 

LAS 

AND 

DAC 

LAV 

DAC* 

DZM* 

JMP 


F0AT5 
C0DE34 

F DAT 10 / EIGHT 

FDAT2 /SET ♦ DAT SLOT 5 EQUAL 

FDAT5 / TO • DAT SLOT 2 

FDAT3 /SET *DAT SLOT 1 EQUAL 

FDATl / TO • DAT SLOT 3 

OJTPT#OJT#INTIM /PREPARE FIRST STORAGE DEVICE 

OUTPT 

-2 /ARE TWO STORAGE 

C0DE21 / DEVICES 

/ REQUIRED ? 

NDKil /NO# OMIT FOLOWING CODE 

MDUC1l /YES# CLEAR JUMP AROlfJD INSTRUCTIONS 

N DK2 / TO USE A SECOND DISK 

NDK3 / FOR STORAGE 

DATIN2.0UT# INTIM /PREPARE SECOND DEVICE FOR USE 

DATIN2 /REMOVE ALL FILES FROM I T 

• * 3 

FDATi /CLEAR .DAT SLOTS FOR THE 

FDAT3 / SECOND STORAGE DEVICE 

/GET THE VALUE OF CONSOLE 
CODE40 / SWITCH *1 

SRPO /SET ADDR. TO THIS VALUE 

-I 

DMV3 /ALLOW SHARE 

DMV2 /SET TT0 TO ACKNOWLEDGE fC 

RDErtH /JUMP AROUND SECOND BUFFER 




BUF2 

/ . 1- - r t 

.BLOCK 

DTBLK< 

RDErtH 

LAC 

CODE!' 


DAC 

HRC 


ISZ 

SKP 

CTMl 


JMP 

OTLP 


•WHITE 

TTO# aj 


• WAIT 

rro 


JMS 

REA DM 

DT 

0 



LAW 

-5 


DAC 

CTMl 


LAW 

-3 


DAC 

CTM2 


DZM 

MIN 


DZM 

HR 

NXTi 

JMS 

CHKN 


JMS* 

.AD 


LAC* 

HRC 


DAC 

SAV2 


TAD 

HR 


DAC 

HR 


ISZ 

CTM2 


JMP 

JPAR 


LAW 



DAC 

CTM2 


LAC 

SAV2 


TAD 

MIN 


DAC 

MIN 

JPAR 

ISZ 

HRC 


ISZ 

CTMl 


JMP 

NXTI 

OTLP 

LAC 

TCA 

MMIN 


TAD 

SMA 

MIN 


JMP 

RDERR 


LAC 

TCA 

CHR 


TAD 

SMA 

HR 


JMP 

R DERR 


LAC 

HH 


DAC 

TIMR 


DAC 

TIME 


/SECOND COLLECTION BUFFER 


•TIMER 

JMS* 

JMP 


360# TOD# S 
CONTL 
« +2 


/INITIALIZE BCD 
/ POINTER 

/HAS FOUR NUMBERS BEEN 
/ HEAD 

/YE S# EXIT FROM ROUTINE 
/NO# ASK FOR TIME 

/READ IN TIME 

/CONTAINS THE ADDR. OF CHAH. READ IN 

/INITIALIZE COUNTER TO EXIT THE 

/ ROUTINE AFTER FIVE 

/INITIALIZE CO IN TER TO GET 

/ THE MINUTES 

/INITIALIZE LOC. THAT SAVE 

/ THE TIME AND MINUTES 

/GET NEXT NUMBER 

/MULTIPLY BY POWERS OF TEN 

/ TO GET BCD EQUIVALENT 

/SAVES THE MINUTES 

/SETS UP THE NUMBERS READ IN 

/ AS THE PRESENT TIME 

/IS THE NUMBER PART OF THE MIN. ? 

/NO# FIRST TWO #*S ARE THE HOURS 

/RESET COUNTER TO GET ALL THE 

/ MIN. CREST OF THE #' S> 

/GET MINUTES AND 
/ SET INTO AN ADD R. 

/ 'WHICH SAVES MIN. 

/GET NEXT MULTIPLYING # 

/OBTAINED 5 NUMBERS ? 

/NO# GET NEXT NUMBER 

/CHECK THE MINUTES 

/IS THE MINUTES GREATER THAN 

/ THE MAX. NUMBER OF MINUTES 

/ IN AN HOUR ? 

/YES# ASK FOR THE TIME AGAIN 
/NO# CHECK THE TOTAL TIME 
/IS THE TIME OF DAY # 

/ LARGER THAN THE BIGGEST # 

/ ALLOW FOR THE TIME OF DAY ? 

/YES# ASK FOR TIME AGAIN 

/NO# PUT THE TIME INTO THE 

/ ADDR. WHICH GIVE THE 

/ TIME OF DAY 

/SET UP THE TIMING R- T. SUB. 

/TRANSFER CONTROL TO CONTROL PROGRAM 




• DSA 

♦ IDLE 

SRPO 

/ADDRESS FOR CALIBRATION SURPRESSION 

MSQT 

2000 

0 

•ASCII 

#, TIME M «I5» 


C TM 1 

0 


/LOC. WHICH COUNTS 5 NUMBERS 

C IMS 

0 


/USED TO IGNORE THE HOURS 

SAV2 

0 


/LOC* TO SAVE THE MIN. #‘S 

HR 

0 


/LOG* TO SAVE THE TIME 

HRC 

• DSA 

• DEC 

DECN 

/POINTER FOR THE BCD MULTIPLIERS 

DECN 

10000 

1000 

100 

10 

I 

• OCT 


/BCD MULTIPLIERS 

CHKN 

0 


/SUB- TO SEPARATE OUT THE NUMBERS 


LAV 

•72 

/PREPARE TO LOOK AT CHARACTER 


TAD* 

DT 

/ HEAD IN 


SAD 

CODE1 

/IS CHAR. A CARRIAGE RETURN ? 


JMP 

RDERH 

/YES# CHECK FOR POSSIBLE ERROR 


ISZ 

DT 

/NO#PREFARE FOR NEXT CHAR. 


SMA 


/IS ASCII CHARACTER LESS THAN 72 


JMP 

.-A 

/NO# GET NEXT CHARACTER 


TAD 

C0DE12 

/YES# IS CHARACTER 


SPA 


/ LARGER THAN 57 ? 


JMP 

♦ •7 

/NO# GET NEXT CHARACTER 

/ 

JMP* 

CHKN 

/YES# CHAR. ISA# SO EXIT 


/jsstitmss ssmutsssttimtisitttitttsmtmstiittfsssuiimst 

• EJECT 

/FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF 


/GET VARIABLES AND PLACE ADDRESSES BELOV 

/JUMP AROUND VARIABLES 

/ADDRESS OF THE BUFFER ADDRESS 

/ADDRESS OF the WORD COUNT 

/ADDRESS OF THE FLAG 

/ZERO FLAG -- VAI T FOR READ IN 

/INSERT DUFFER ADDRESS INTO THE 

/ A/D CALL ROUTINE 

/INSERT THE WORD COUNT INTO THE 

/ A/D CALL ROUTINE 

/INSERT THE FLAG ADDRESS INTO THE 
/ A/D CALL ROUTINE 

/THE A/D CALL ROUTINE t TO INIT. READ > t 

/NUMBER OF SAMPLES TO TAKE 

/BUF. ADDRESS IN WHICH JO STORE SAMPLES 

/COMPLETION AND ERROR FLAG ADDRESS 

/R-T SUB » FOR INT. SEHV. ROUT* TO GO TO 

/RETURN TO FORTRAN PROGRAM 

/SUBR. TO WRITE LIN • TABLE ON -DAT ” TBI *' 


/ADDH* of addh- of lin. table 

/GET ADDH. OF LIN- 
/ TABLE 

/PREPARE STORAGE DEVICE 
/OPEN FILE FOR TABLE 
/PUT ADDRESS OF TABLE INTO 
/ WRITE COMMAND 
4 /PUT TABLE ON STORAGE DEVICE 


RADC 0 

JMS* 

JMP 

NBIA 0 

NB2C 0 

NB3F 0 

DEM* 

LAC* 

DAC 

LAC* 

DAC 

LAC 

DAC 

JMS 

NBA 0 

NBS U 

NB6 a 

NB7 0 

,JMP* 

TTM 0 

JMS* 

JMP 

STA 0 

REP LAC* 

DAC 
. INIT 
. ENTER 
LAC 
DAC 

* WRITE 

• WAIT 
-CLOSE 
JMP* 


-DA 
• ♦A 


NB3F 

NBIA 

NB5 

NB2C 

NBA 

NB3F 

NB6 

ADHEAD 


RADC 

- DA 

*♦2 

STA 

STA 

TBI -OUT# TTM 
TBI# STT 
STA 

• + 3 

TBI # DUMP# 0# 51 

TBI 

TBI 

TTM 


/ 


FFFFKFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFKFFFFFFFFFKF 
• EJECT 


FFFFFFFFFFFFF 


/ 

ONCA 


DZM* 

SUM4 

LAC* 

TC3 

DAC 

TC3 

LAC* 

rc3 

SPA 


JMP 

LETMON 

TAD 

C0DE24 

SPA 


JMP 

VINT 

TAD 

C0DE25 

SPA 


JMP 

ONC3 

TAD 

C0DE25 


♦ + ***+++ +♦ 44 . 4 4 44 444 * ♦♦+♦♦+++♦++++++ 


/INITIALIZE TIME OF TEAR LOC- 
/GET ADDRESS OF THE 
/ DATE 

/GET FIRST 2 1/2 CHAftAC TER5 
/ARE THEY LETTERS? 

/YES# CHECK FOR WORDS 
/NO# CHECK THE NUMBERS 
/IS THE FIRST # A ONE? 

/YES# CHECK FOR WINTER MONTHS 
/NO# CHECK FOR SUMMER MONTHS 
/IS THE # >A ? 

/NO# EXIT 
/YES 



/is THE NUMBER « 8 7 
/NO# EXIT 

/YES# SET TIME OF YEAR LOC# TO SUM. 
/EXIT 

/GET SECOND NUMBER 

/IS THE NUMBER AN ASCII ZERO C 60 ) 7 
/YES# EXIT 

/NO# SET TIME OF YEAR LOC. 

/ TO WINTER 
/EXIT 

/IS THE FIRST LETTER 
/ AN "A" 7 
/YES# EXIT 

/NO# IS THE FIRST LETTER 

/ A "D" 7 

/YES# WINTER MONTH 
/NO# IS THE FIRST LETTER 
/ A ,# F M 7 
/YES# EXIT 
/NO# IS THE LETTER 
/ A "J M ? 

/YES# LOOK AT SECOND CHARACTER 
/NO# IS THE LETTER 
/ A ”M M ? 

/YES# LOOK AT THIRD LETTER 
/NO# IS THE LETTER 
/ A ”N M 7 

/NO# EXIT 

/SET TIME OF YEAR LOC. 

/ TO WINTER 

/EXIT 

/CHECK IF SECOND LETTER 
/ IS AN M A W 

/YES# A WINTER MONTH 

/NO# SET TIME OF YEAH LOC* TO SOMMER 
/EXIT 

/IS THE THIRD LETTER 
/ A M Y” ? 

/YES# SET TIME OF YEAR TO SUMMER 
/EXIT 

/ 

/ ++*+++++ ♦♦♦ + ♦* + ♦♦♦+♦♦+♦*♦♦♦♦♦♦♦♦♦♦+♦♦♦♦ + ■► 

.EJECT 

/DDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD 
/ SUBROUTINE CLEARS I/O DEVICES FROM MEMORY TO BE ABLE TO 

/DISALLOW SHARING THE DATA COLLECTION DEVICES 


/ 

FRDT 

0 





LAC* 

DMV6 

/GET LOC. OF THE .IOIN TABLE 



DAC 

DUM2 

/SAVE THE LOCATION 



TAD 

CODE1 1 

/GO TO FOREGROUND DEVICES 



DAC 

DUMI 

/ AND SAVE THAT LOC. 



LAC* 

DUM2 

/GET THE NEQ. OF THE # FOREGND 



TAD 

C0DE34 

/ DEVICES AND DOUBLE THE § 



RAL 


/ (2 WORDS PER DEVICE USED) 



DAC 

DUM2 

/ AND USE AS THE COUNTER 


REDOl 

LAC* 

DJMt 

/GET FIRST WORD FOR THE FOREGND 

DEV 


AND 

CODE7 

/CHECK FOR THE DEVICE * 



SAD 

CODE 1 3 

/IS IT A 5 (DISK) ? 



JMP 

DKS 

/YES# CHECK WHICH DISK 



SAD 

CODEU 

/NO# IS IT A 4 (DECTAPE) 7 



JMP 

DCTP 

/YES# CHECK WHICH DECTAPE 



SAD 

C0DE39 

/NO# IS IT A 7 (P. PUNCH) 7 



JMP 

SAVNI 

/YES# SAVE IT 


DLTDS 

LAC 

CODE37 

/DELETE ALL OTHER DEVICES 



DAC* 

DUMI 

/ BY INSERTING A DUMMY NAME 



TAD 

CODE22 

/SET UP THE NEXT DUMMY 



DAC 

CODE37 

/ NAME 



ISZ 

DUMI 

/GO TO NEXT WORD 



LAC 

CODE36 

/INSERT SECOND DUMMY 



DAC* 

DUMI 

/ NAME INTO THE TABLE 


ANOl 

ISZ 

DUMI 

/GO TO NEXT DEVICE 



ISZ 

DUM2 

/HAS ALL THE DEVICES BEEN CHECK 

? 


JMP 

REDOl 

/NO# CONTINUE 



JMP* 

FRDT 

/'YES# RETURN 


DKS 

LAC* 

DUMI 

/GET FIRST WORD 



AND 

C0DE15 

/LOOK AT UNIT NUMBER 



SAD 

COD £22 

/IS IT A 1 ? 



JMP 

DLTDS 

/YES# DELETE THE DEVICE WORDS 



SAD 

C0DE23 

/IS IT A 3 7 



$KP 


/YES# CHECK FOR DOUBLE DEVICE COLL. 


JMP 

SAVNI 

/NO# SAVE THE TWO WORDS 




SMA 

JMP 

ONC3 


ISZ* 

SUM4 


JMP 

0NC3 

MINT 

AND 

C0DE26 


SAD 

CODE27 


JMP 

0NC3 


LAV 

-I 


DAC* 

SUM 4 


JMP 

ONC3 

LETMON 

TAD 

C0DE28 


SPA 

JMP 

ONC3 


TAD 

C0DE29 


SPA 

JMP 

WINTL 


TAD 

CODE30 


SPA 

JMP 

0NC3 


TAD 

C0DE29 


SPA 

JMP 

SVDIS 


TAD 

C0DE29 


SPA 

JMP 

EWD1S 


TAD 

CODE31 


SMA 

JMP 

0NC3 

■J li'JTL 

LAW 

-1 


DAC* 

SUM4 


JMP 

QNC3 

SV7D1S 

AND 

C0DE26 


SAD 

C0DE32 


JM? 

WINTL 


ISZ* 

SJM4 


JMP 

ONC3 

EV7DIS 

AND 

C0DE33 


SAD 

C0DE34 


ISZ* 

SJM4 


JMP 

DMC3 
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LAC 

NDK2 

/HAS THE LOC. BEEN CHANGED 


SAD 

MDUM1 

/ TO USE 2 DEVICES ? 


JMP 

SAVN1 

/YES. SAVE 

DC TP 

JMP 

DLTDS 

/NO. DELETE LOCATIONS 

LAC* 

DUMl 

/GET FIRST WORD AGAIN 


AND 

CODEI 5 

/CHECK THE INIT NUMBER 


TAD 

C0DE36 

/IS THE INIT 0 LESS 


SPA 


/ THAN 4 7 


JMP 

DLTDS 

/YES. DELETE IT 

SAVN1 

1SZ 

DUMl 

/NO. GO TO NEXT WORD 

y 

JMP 

ANOl 

/GO TO NEXT DEVICE 

/ddddddddddddddddddddddddddddddddddddddddddddddddddddddddddodddddddddd 

J 

• EJECT 



r 

CAL ERR 

.WRITE 

TTO.ASC. 

ERRCAL.0 /LIN. TABLE ERROR 


.WAIT 

TTO 


JM$* 

JMP 


TBFORL 

CALTB 


/RECALI BRATE 


BUF3 .BLOCK DTBLK-342 /PRQCESSI NQ f S BUFFER 

/nmtnnmxmmmmxxmnmxxxmmxsmxxmtmxxm: 


/ 

ONC3 


DBA 

CKDB 


NLTM1 


NUM2 


DBI MC 


• INIT TBI .IN. RES 


/LINEARIZATION TABLE IN 


LAC* 

DAC 

LAS 

AND 

SZA 

JMP 

•WRITE 


DB 

DB 

CODE40 

ENDDB 
TTO. ASC. 


ENDDB 


JMS 

REA DM 

0 

LAV 

-40 

TAD* 

DB0 

ISZ 

DB0 

SPA 

JMP 

ENDDB 

TAD 

CODE17 

SMA 

JMP 

NUM1 

TAD 

CODEI 2 

SPA 

JMP 

NUMi 

ISZ 

DBCNT2 

SKP 

JMP 

NUM2 

DAC 

DBP1 

JMP 

CKDB 

ISZ 

DBCNT2 

JMP 

resetd 

LAC 

DBP1 

DAC* 

DB 

TAD 

C0DE3 

CLLIRAL 

DAC* 

DBP 

JMP 

DBI NC 

LAW 

-2 

DAC 

DBCNT2 

JMP 

CKDB 

DAC* 

DB 

TAD 

CODE10 

CLLIRAL 

DAC 

DBP 8 

LAC 

DBP | 

JMS* 

• AD 

LAC 

CODEI 2 

TAD* 

DB 

DAC* 

DB 

LAC 

DBP 1 

TAD 

CODEI 0 

SVHA 
CLL 1 RAR 

TAD 

OBP2 

DAC* 

DBP 

LAW 

-2 

DAC 

DBCNT2 

ISZ 

DB 

ISZ 

DBP 

ISZ 

DBCNT1 

ISZ 

DBCNT3 

JMP 

CKDB 

LAC 

DBO 


/GET ADDR. OF DB SETTINGS FOR 
/ THE FORTRAN PROG . 

/CHECK CONSOLE SWITCH 
/ 0 \ 

/IS IT SET 7 

/YES. USE DEFAULT DB SETTING 
MSGDB.0 /DB MESSAGE 

/GET RESPONSE 
/LOCATION OF THE RESPONSE 
/SET ASCII CHAR. LESS THAN 
/ AO TO A NEG. 0 
/PREPARE FOR NEXT CHAR. 

/IS CHAR. > 40 ? 

/NO. THERE ARE NO MORE CHAR. 

/NO. IS THE ASCII 0 LESS THAN 

/ THAN 72 7 

/NO# USE ONLY* ONE NUMBER 

/YES. IS THE ASCII CHAR. AN 

/ ASCII * ( >57 ) 7 

/NO. USE ONLY ONE NUMBER 

/YES* 1$ THIS THE SECOND # ? 

/NO# CONTINUE 

/YES. PROCESS THE TWO S 

/SAWE THE FIRST 0 

/ AND GET SECOND 0 

/HAS ONE 0 BEEN OBTAINED 7 

/NO. RESET COUNTER 

/YES. GET THE NUMBER AND 

/ SET INTO FORTRAN ARRAY 

/SET THE NUMBER UP AS A SPACE 

/ AND A 0 IN ASCII FORMAT 

/USED TO PRINT OUT DB MESSAGE 

/PREPARE FOR NEXT DB 

/RESET COUNTER FOR 

/ S NUMBERS 

/TRY AGAIN 

/SAVE SECOND I IN FORTRAN ARRAY 

/SET LP NUMBER IN 

/ ASCII CODE 

/ AND STORE 

/GET FIRST NUMBER 

/SET UP THE BCD 

/ EOUI VALENT 

/ADD TO THE SECOND # 

/ AND STORE IN FORTRAN PROG. 

/GET FIRST NUMBER AGAIN 
/SET UP THE NUMBER 
/ IN ASCII 
/ CODE 

/ADD TO PRIVIOUS 0 TO FORM 
/ THE D8 SETTING IN ASCII CODE 
/RESET COUNTER FOR THE 
/ NEXT 2 NUMBERS 
/NEXT LOG- IN THE FORTRAN ARRAY 
/NEXT LOC. IN THE MACRO ARRAY 
/COUNTER TELLING # OF DB SETTINGS 
/HAS FOUR DB SETTINGS BEEN OBTAINED ? 
/NO. GET NEXT SETTING 
/CHECK— HAS ANY NUMBERS 
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DEFDB 
/- 

CALTB 

CALLP 


DLVT 

RUNTlM 

G ETHrt 


SNA 


/ 

BEEN READ IN 1 

JMP 

DEFDB 

/NO# 

USE DEFAULT SB 

LAV 

-I 

/YES# ALLOW POSSIBLY ONE MORE 

DAC 

D9CNT3 

/ 

DB SETTING 

ISZ 

DBCNT2 

/IS 

ONE # STILL UNPROCESSED 7 

SKP 


/NO/ 

CONTINUE THE EXIT 

JMP 

NUM1+8 

/YES/ PROCESS THE LAST # 

DAC* 

DBP 

/SET 

THE LAST LOC. TO A NEG* # 

LAC 

DBCNTi 

/SET 

THE FORTRAN COUNTER 

DAC* 

DBC 

/ 

TO THE # OF DB SETTINGS 

TAD 

CODES 

/IF 

THERE IS ONLY ONE DB 

HAL 


/ 

SETTING SET THE LINK 

LAC 

DMVl 

/SET 

UP JIMP AROUND COLL* READ IN 

SZL 


/IS 

THERE MORE THAN ONE DB SETTING 

DAC 

UP DAT 1 

/NO/ 

INSERT JUMP AROUND 

JMP 

*♦3 

/USE 

IMPUTED DB SETTINGS 

LAC 

C0DE6 

/SET 

IP ASCII CODE FOR 

DAC 

DBO 

/ 

ZERO DB SETTING 


LAC 

DAC 

DZM 

LAC 

DAC 

LAC* 

DAC 

« SEEK 

• READ 

• VAI T 
•CLOSE 
LAW 
DAC 
LAC 
DAC 
LAC 
TCA 
TAD* 
SPA 
JMP 
LAC* 

I sz 

ISZ 

JMP 

DZM* 

LAS 

AND 

SZA 

JMP 

•WRITE 

•WRITE 

•WRITE 

♦WRITE 

• WAIT 
JMS 
JMS 

0 

■ I N I T 

• I N I T 
LAC 
DAC 
DZM 
LAW 
DAC 
LAW 


CODE20 

MINEQC 

TCI 

CDBO 

DSP 

DBP 

M5G9-M 1 
TBl/STT 

TBI # DUMP# ST1 / 51 < 

TBI 

TBI 

-775 

DUM2 

< ST1 ♦ l 

DUM3 

ST1 


/NO/INITIALIZE MULTIPLIER 
/ TO CONVERT #’S TO MIN* 

/INITIALIZE STORAGE ADDR* FOR CONVERTING 
/INITIALIZE DB SETTINGS TO 
/ THE BEGINING 

/SET UP COLLECTION’S DB 
/ MESSAGE 

/READ TABLE IN 


DUM3 

CAL ERR 

DUM3 

DUM3 

DUM2 

CALLP 

FDAT10 

CODE40 

DLWT 

TTO/ ASC# MSG7/9 

TTO/ ASC/MS38/0 

TTQ/ASC/MSG9/3 

TTO/ASC#MSG10/3 

TTO 

DBSUB 

REA DM 

TTO/ OUT/ P START 

TTI/ IN/? START 

CODE35 

TC3 

D5T0K 

-3 

TC2 

-40 


/PREPARES COUNTER TO CHECK 
/ THE VALIDITY OF THE LIN. TABLE 
/SET POINTER TO THE SECOND 
/ LOCATION OF THE TABLE 
/GET THE FIRST NUMBER AND 
/ COMPLIMENT 

/IS THE PREVIOUS NUMBER 
/ LARGER? 

/YES# ERRONEOUS TABLE 

/GET PRENSENT LIN* NUMBER 

/GO TO NEXT NUMBER 

/IS IT THE END OF THE TABLE ? 

/NO# CHECK NEXT NUMBER 
/YES/ DELETE • DA T SLOT 10 
/GET THE # 1 DATA SWICH FROM 
/ CONSOLE 

/IS IT SET 7 

/YES* SUHPHESS PRINT OUT 


/NUMBER OF HOURS 

/SET UP NEXT DB SETTING 
/HEAD IN # OF HOURS 

/CONTAINS ADDRESS OF THE #*S READ IN 
/RESET THE »P RESTART 
/ ADDRESS 

/DEFAULT VALUE FOR THE LENGTH 
/ OF HUN C 20 HOURS) 

/CHECKS # OF HEAD IN 

/PREPARE TO USE NO MORE 
/ THAN THREE NUMBERS 
/CHECK IF CHARACTER READ IN 


TAD* RUN TIM 

DZM* RUN TIM 

ISZ RUNTIM 

SPA 

JMP FNUM 

SAD CODE16 

JMP DECPT 

TAD CODE! 7 

SMA 

JMP GETHR 

TAD C0DE12 

SPA 

JMP GETHR 

JMS* • AD 

LAC* MINEQC 

TAD TCI 

DAC TCI 

ISZ MINEQC 

ISZ DSTOR 

I SZ TC2 

JMP GETHR 


/ IS AN ASCII t OR PERI 00 
/i*ESET L0C. IN CHAR. BLOCK 
/GO TO NEXT CHARACTER 
/IS ASCII CHARACTER < 40 ? 

/YES/ EXIT 

/NO/ IS IT A DECIMAL POINT ? 

/YES/ GO TO DECIMAL POINT ROUTINE 
/NO/ CHECK FOR AN ASCII NifIGER 
/IS THE CHAR. < 72 ? 

/NO/ GET NEXT CHAR. 

/YES/ IS CHARACTER GREATER 
/ THAN 57 OCTAL ? 

/NO/ GET NEXT CHARACTER 

/YES# GET THE BINARY COOED DECIMAL 

/ EQUIVALENT OF THE NUMBER 

/ADD TO THE PftECEEDING NUMBERS 

/ AND SAVE 

/SET UP NEXT MULTIPLIER 

/INCREMENT NUMBER COUNTER 

/HAS ENOUGH #’ $ BEEN OBTAINED 7 

/NO/ GET NEXT NUMBER 



DECPT 


FNUM 


CNUM 


NO NUM 


/ 


JMP 

CNUM 

LAW 

-I 

DAC 

TC2 

LAW 

m 2 

TAD 

DSTOR 

SMA 


JMP 

GETKR 

LAC 

TCI 

CLL 


IDIV 


12 


LACO 


DAC 

TCI 

LAC 

C0DE2 

DAC 

MINEQC 

ISZ 

DSTOR 

JMP 

GETHR 

LAW 

-2 

TAD 

DSTOR 

SMA 


JMP 

CNUM 

I AC 


SZA 


JMP 

NONUM 

LAC 

TCI 

CLL 


IDIV 


12 


LACQ 


DAC 

TCI 

LAC 

TCI 

TCA 


DAC 

TC3 

DZM 

TCi 

DZM 

TC2 

DZM 

MCNTR 

LAC 

CODE! 8 

DAC 

DSTOR 

JMS 

FRDT 

.TIMER 

0# BEG I N# 6 

JMS* 

PROC 

JMP 

. + 4 

. D5A 

SUM 4 

. D$A 

ALF 

• DSA 

TIMR 

• IDLE 



/YES- EXIT 

/SET COUNTER TO GET ONLY 
/ ONE MORE NUMBER 
/CHECK THE NUMBER 
/ COUNTER 

/HAS TWO NUMBERS BEEN OBTAINED 7 

/YES# GET THE LAST NUMBER 

/NO# IF THERE WAS A NtWSER# IT IS 

/ OFF BY A FACTOR OF TEN TOO HIGH 

/THEREFORE REDUCE THE NUMBER BY A 

/ FACTOR OF TEN (ZERO IS UNEFFECTED) 

/GET THE QUIOTIENT 

/REPLACE VI TH CORRECTED t 

/SET THE BCD POINTER TO 

/ THE LAST MULTIPLIER 

/INCREMENT NUMBER COlWTER 

/GET NEXT NUMBER 

/CHECK THE NUMBER 

/ COUNTER 

/HAS 2 NUMBERS BEEN READ IN 7 

/YES# IGNORE THE FOLLOWING CODE 
/HAS EVEN ONE NUMBER BEEN 
/ READ IN 7 

/NO# USE THE DEFAULT VALUE (13-65 HR) 
/YES# THERE IS ONE NUMBER BUT 
/ IT IS THE WRONG BCD 

/ EQUIVALENT SO REDUCE IT BY 

/ A FACTOR OF TEN 
/GET THE INTEGER ANSWER 
/SAVE CORRECTED NUMBER 
/GET LENGTH OF RUN NUMBER 
/ COMPLIMENT IT 
/ AND SAVE THE NEGATIVE 

/INITIALIZE LOC* TO DETERMINE 
/ THE END OF THE RUN 
/INITIALIZE BINARY TIME OF DAY 
/SET UP MAX* STORAGE FOR 
/ STORAGE DEVICE 
/FREE BACKGROUND DEVICES 
/SET UP MAIN REAL TIME SUB* 

/START PROCESSING PROG. 

/# OF PARAMETERS +1 
/THE SEASON OF THE YEAR 
/LINE FEED 

/TIME AT THE END OF EACH FILE 


/uzniinninnummummnxnmnxHimnxnnmiuu 

END 2000 

0 

.ASCII <14>"***END OF PROCESS 2 NG***"<I5> 

S TPM 2000 

0 

•ASCII «7»<7><1S> 

•EJECT 

/RRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRRHRHRRRRRRRRRHRRRRRRRRRRRRHRRHR 

/ THIS SUBROUTINE IS USED TO READ IN CHARACTERS FROM THE 

/TELETYPE ONE AT A TIME UNTIL A CARRIAGE RETURN IS FOUND* THE 
/CHARACTERS ARE STORED IN "NBLK" (A BLOCK OF LENGTH 30 >. THE 
/ADDRESS OF "NBLK" IS PUT INTO THE LOCATION AFTER THE "JMS REA DM" 
/INSTRUCTION USED TO CALL THIS ROUTINE. THE PROGRAM CONTROL IS 
/THEN RETURNED TO TWO LOCATIONS AFTER THE CALLING INSTRUCTION. 

READM 0 


CONR 


NUMB 

READM 

READM 

TTI # I A#MTTI # 3 
TTI 

MTTI+2 

NUMB 

(177 

HOU 

NUMB 

(25 

DLT 

<15 


LAC 

DAC* 

ISZ 

• READ 

• WAIT 
LAC 
DAC* 
SAD 
JMP 
ISZ 
SAD 
JMP 
SAD 
SKP 
JMP 
LAC 
DAC 
LAC 
DAC 


CONR 

(NBLK 

NUMB 

(13 

RBGT+2 


/SET THE ADDR. OF THE BLOCK CONTAINING 
/ THE READ IN CHAR. INTO THE CALLING 
/ LOC. + 1 AND RETURN TO LOC. +2 
/READ IN ONE CHARACTER 

/GET THE CHARACTER 
/STORE IN THE CHAR. BLOCK 
/IS THE CHAR* A RU80UT ? 

/YES# DELETE PRIV1QUS CHAR. 

/NO# PREPARE FOR NEXT CHAR. 

/IS THE CHAR. A rU ? 

/YES# DELETE THE LINE 

/NO# IS THE CHAR. A CARRIAGE RETURN ? 

/YES# EXIT FROM THE READ LOOP 

/NO# GET NEXT CHAR. 

/INITIALIZE CHARACTER BLOCK 
/ POINTER 

/WRITE OUT ON TTY 
/ A LINE FEED 



JMS 

JMP* 

DLT LAC 

DAC 
LAC 
JMP 

ROU LAC 

DAC* 
LAV 
TAD 
SAD 
JMP 
DAC 
LAC 
DAC 
JMS 
JMP 

VRT 0 

•WRITE 
• VAIT 
JMP* 

RBOT 2003 

0 
0 

MTU • block 

N UMB * DSA 

NBLK .BLOCK 


VRT 
REA DM 
<N8LK 
NUMB 
<100 
.♦tl 
US 
NIWS 
-1 

NUMB 
(NUMB 
CON ft 
NUMB 
<134 
RBOT* 8 
VRT 
CONK 

TTO# I A. RBOT* 3 

TTO 

VRT 


3 

NBLK 

30 


/ONE LETTER VRITE ROUTINE 
/RETURN TO CALLING ADDR. +8 
/RESET POINTER TO WE BEGINING 
/ OF THE BLOCK 

/SET IF TO VRITE 
/ AN • CHAR. 

/SET A CARRIAGE RETURN INTO WE 

/ PRESENT LOG. OF WE CHAR. BLOCK 

/GO BACK TO VRITE OVER 
/ WE PREVIOUS CHAR. 

/IS WE POINTER AT WE BEGINING ? 
/YES# DO NOT PRINT OUT ANY WING 
/NO# RESET POINTER BACK 1 LOG* 
/PREPARE TO VRITE A "RUBOUT" CHAR. 
/SET IN OUTPUT ADDR. 

/VRITE OUT WE ONE CHAR. 

/READ NEXT CHAR. 

/VRITE ONE IMAGE ALPHANUM. CHAR. 
/RETURN 

/IMAGE A. OUT PUT 
/ MESSAGE 

/LOC. TO STORE CHAR. TO OUTPUT 
/BLOCK CHAR. IS READ INTO. 

/CHARACTER BLOCK POINTER 

/BLOCK TO STORE ALL THE READ IN CHAR. 


/rrrrrHRRRHRRRRRRRRRRRRRRRRHRRRRHRRRRRRRRRRRRRRRBRRRRRRRHRRRRRHRRR 
/TTTTTTTT T TTTTTTTTTTTTTTTTTTTTTTTT TTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTT 
✓ THE R.-T. SUBROUTINE THAT KEEPS TRACK OF THE TIME OF DAY. 

/THE SUBROUT INE IS INITIATED BY A ".TIMER" WHILE STORING THE PRESENT 
/TIME IN LOCATION "TIME" AND THE MINUTES * TEN IN "MIN". WE ROUTINE 
/RESTARTS IT SELF STORING THE (F DATED TIME IN "TIME" AND THE UPDATED 
/MINUTES TIMES TEN IN "MIN" USING BINARY CODED DECIMAL. THE BINARY 
/REPRESENTATION OF THE NUMBER OF TENTHS OF MINUTES THAT HAS PAST IS 
/GIVEN IN LOCATION "MCNTR". 


TOD 

0 



DAC 

SVAC 


.TIMER 

360# TOD# 5 


ISZ 

MCNTR 


LAC 

MIN 


SAD 

MMIN 


JMP 

CMIN 


ISZ 

MIN 


ISZ 

TIME 

RESTOR 

LAC 

SVAC 


.RLXIT 

TOD 

CMIN 

DZM 

MIN 


LAC 

TIME 


SAD 

CHR 


JMP 

CTIME 


TAD 

NXHR 


DAC 

TIME 


JMP 

RESTOR 

CTIME 

DZM 

TIME 


JMP 

RESTOR 

/TTTTTTTTTT TTTTTTTTTTTTTTTTTTTT 


.EJECT 


P START 

DZM 

BEGIN 


LAC 

MDUM1 


DAC 

START 


LAC 

IDCOU 


SPA 



JMP 

♦ *3 


LAC 

< JMP RESTAR 


DAC 

START 


•TIMER 

0#BEi3IN#6 


JMP 

DUM7 + 1 

✓ MMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMt 

/ MAIN PROGRAM 

/ 



BEGIN 

0 



DAC 

SAVEAC 

START 

NOP 


Cl 

• INI T 

OUTPT# OUT# RESTAR 


JMP 

INIT 

ft DO 

JMS 

ADREAD 


tnsam 



BUF 



ERS 



500000+R5U8 


/SAVE AC FROM LOWER APL LEVELS 
/HESTART CLOCK FOR ANOWER 6 SECONDS 
✓INCREMENT BINARY COUNTER 
✓CHECK THE MINUTE ADDR. 

✓IS IT SET TO 59.9 MINUTES 7 
/YES# SET UP THE NEW HOUR 
✓NO# INCREMENT THE MINUTE COUNTER 
✓ AND TIME OF DAY COUNTER 

/RESTOR AC FOR LOVER LEVEL PROGRAM 
/EXIT FROM PROG. AND API LEVEL 
/RESET MINUTES TO EERO 
✓CHECK TIME OF DAY 

/IS THE TIME AT THE END OF THE DAY ? 

/YES# START A NEW DAY 

/NO# CHANGE TO NEXT HOUR 

/ AND SAVE THE TIME 

/EXIT 

/RESTART TIME TO A NEW DAY 
/ AND EXIT 


/RESET R.-T. SUB. BEGIN TO BE REENTERED 
/PREPARE TO REINITIALIZE THE STORAGE 
/ DEVICE AND BEGIN A NEW FILE 

/HAS THE COLLECTION FINISHED COLLECTING 

✓ THE PRESENT FILE ? 

✓YE5# LEAVE AS SET UP 
✓NO# PREPARE TO CLOSE THE 

✓ OLD. FILE 
/RESTART THE COLLECTION 
/RESTART THE PROCESSING PROGRAMS 


/ENTERANCE TO THE MAIN REAL TIME SUB. 

✓ USED TO CONTROL PROGRAM FLOW 
✓DT OUT 

/GO TO INITIALIZING ROUTINE 
/PREPARE TO READ FROM A-D CONVERTER 
/# OF DATA NUMBERS TO READ 
/ADDR. TO STORE DATA #' S 
/COMPLETION AND ERROR TIMING FLAG 
/PRIORITY LEVEL ♦ R- T SUBft. TO EXEC. 



HET3 

LAV 

-1 

/LOAD -I INTO MEMORY TO KEEP TRACK OF 


DAC 

CNT2 

/ THE SAMPLE THAT HAS BEEN COLL. 


LAC 

SAVEAC 

/RESTORE ACCUMULATOR 

END! 

.KLXIT 

BEG I N 

/RELINQUISH CONTROL TO LOWER PRIORITY 

ONCE 

NOP 



.TIMER 

9# CHECK# 5 

/SET CLOCK TO WAIT 9/60 SECONDS 

PKO 

JMS 

ADREAD 

/A-D CONV. READ FOR X- SAMPLES 


TNSAM 


/THE VARIABLES USED 


buf 


/ ARE THE SAME ONES USED FOR 


Ert2 


/ THE O-SAMPLE AND ARE 


500000+RSUB 

/ EXPLAINED ABOVE 


LAC 

SAVEAC 


EN02 

.RLXIT 

BEG IN 

/RELINQUISH CONTROL 

HET2 

LAC 

TRANF 

/GET DECTAPE TRANSFER FLAG 


SZA 


/ AND TEST IT 


JMS 

DTRANS 

/FLAG SET - TRANSFER DATA 

A 

CLAICLL 


/CLEAR AC TO READ CONSOLE SWITCHES 


LAS 


/GET # FROM CONSOLE'S DATA SWITCHES 


HAL 


/PUT BIT 00 OF AC IN L 


SNL 


/IS THE LINK A ZERO ? 


JMP 

RDO 

/YES# L=3 COLLECT DATA 


•TIMER 

1 20 # W T 1 # 6 

/NO# STOP COLL. AND GIVE TIME TO A 


.RLXIT 

BEGIN 

/ LOWER LEVEL BEFORE RECHEKIN3 SWITCH 

WT1 

0 



DAC 

SAVEAC 



LAC 

MDUM1 

/PUT CLOCK BACK INTO 


DAC 

ONCE 

/ OPERATION 


DZM 

CNT 

/TELLS PROC. COLL. HAS STOPPED 


LAC 

(JMP A 

/SET UP TO RECKECK CONSOLE SWITCHES 


DAC 

START 

/ AND PUT INTO START 


LAC 

SAVEAC 



.TIMER 

0 # BEG IN# 6 

/RETURN TO R.-T. SUB. BEGIN 


.RLXIT 

WT1 


/ END 

OF MAIN PROGRAM 



/MMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMrtMMMMMMMMMMMMMMMMMMMM 

BUT .BLOCK DATBLK /BUFFER TO STORE AD SAMPLES 

NAME .SIXBT M DATAF 1 DAT'* /NAME OF FILE TO STORE DATA 

/ INITIALIZING ROUTINE 


LAC 

(BUFl 

DAC 

POINT 

DAC 

TSUF 

LAC 

IDCOU 

DAC 

CN1 

DZM 

IDCOU 

DZM 

TKANF 

DZM 

CNT 

LAC 

MDJM1 

DAC 

ONCE 

LAW 

-2 

DAC 

CTT6 

LAC 

(777 

DAC 

MAX A 

LAW 

MNCC 

DAC 

CTT2 

LAC 

MONTH 

TAD 

TCI 

DAC 

DBG 

TAD 

TC2 

DAC 

TC2 

TAD 

TC3 

SMA 


JMP 

EXTIP 


/POINTER IN DECTAPE BUFFER 
/NAME OF DECTAPE BUFFER IN USE 
/SAVE THE ID # (THE # 

/ OF FRAMES PER FILE) 

/ID NUMBER 

/DECTAPE TRANSFER FLAG 

/COUNTER FOR CLOCK 

/PUT CLOCK INTO OPERATION 

/SET BUFl AS THE FIRST 
/ BUFFER TO BE USED 
/INI T. NOISE MAX. 

/ LOCATION 

/INITIALIZE COUNTER FOR MAXIMUM 
/ ALLOWABLE NOISE 
/GET THE LENGTH OF TIME REQUIRED 
/ COLLECT THE PRIVIOUS FILE 
/SAVE IT 

/ADD TO ALL OTHER PRIVIOUS TIMES TO 
/ COLLECT THE OTHER FILES# SAVE AND 
/ COMPARE TO THE MAXIMUM TIME 
/ARE THE FILE TIMES LARGER ? 

/YES# STOP COLLECTION 


/ DECTAPE FILE ROUTINE 
/ STORE DATA IN FILE ACCORDING 
/ 


UP 

■ FSTAT 

OUTPT#NAME 

UR 1 

SKP 



JMP 

update 

WRITE 

LAC 

TCA 

MCNTR 


DAC 

TCI 


.DLETE 

OUTPT#NAME 

C 3 

.ENTER 

OUTPT#NAME 


JMS 

DTRANS 


JMS 

DTRANS 


LAC 

MDUM2 


DAC 

HEPL 


JMP 

RDO 


TO RESPONSE 

/CHECK FOR COLLECTED FILE 
/(REPLACED BY "SZA">IS FILE PRESENT? 
/YES# ACKNOWLEDGE THE PRESENCE 
/GET THE BEGINNING TIME FROM 
/ THE BINARY CLOCK COUNTER 

/ TO DETERMINE THE COLL. TIME 

/DELETE FILE IF PRESENT 
/OPEN FILE 
/WRITE DUMMY BLOCK 
/TWICE 

/TELLS PROCESSING THAT COLLECTION 
/ HAS STARTED COLLECTING A FILE 
/RETURN 
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UPDATE 

LAC 

TIME 

/SETS THE TIME TO WE END 


TAD 

(S 

/ OF THE COLLECTED FILE AND 


DAC 

TIMR 

/ AND ROUNDS OFF TO WE NEAREST MIN. 


.WRITE 

TTO#ASC#MSG1«0 

/FILE PRESENT 


.WAIT 

TTO 


UP DAT 1 

NOP 


/REPLACED BY "JMP BLK2” FOR 1 DB SET* 


• WRITE 

TTO. ASC.M5G2. 34 

/KEEP IT? 


• WAIT 

TTO 


BLK1 

JMS 

REA DM 

/READ RESPONSE 

COM 

0 


/ADDRESS OF WE RESPONSE 

3 ETCH1 

LAC* 

COM 

/GET READ IN CHARACTER 


DEM* 

COM 

/ AND ZERO THE LOC. 


SAD 

(116 

/IS CHARACTER A "N"7 


JMP 

STAQN 

/YES 


tad 

(•72 

/CHECK IF NUM. IS LESS 


SMA 


/ THAN 7 2 OCTAL 


JMP 

BLK3 

/NO. CONTINUE 


TAD 

(12 

/YES. IS NUM. GREATER THAN 


SPA 


/ 57 OCTAL? 


JMP 

BLK3 

/NO. CONTINUE 


DAC* 

COM 

/YES. SAVE 


ISZ 

COM 

/NEXT CHARACTER 


JMP 

3ETCHI 

/REPEAT 

BLK3 

LAC 

L ETC HQ 

/INITIALIZE THE BLOCK 


DAC 

COM 

/ CONTAINING THE INPUT CHAR. 


LAC* 

COM 

/GET FIRST CHAR. 


SNA 


/X 5 IT A ZERO ? 


JMP 

CONDB 

/YES. IGNORE IT 


DAC* 

CHG 

/PASS TO PROC THE DB CHANGE 


TAD 

<•1 

/OFFSET THE # BY -1 


CLLfRAR 


/DIVIDE BY TWO AND SAVE REMAINDER 


TAD 

(60 

/SET UP AS ASCII CHAR. 


SWKA 


/PHEPARE FOR MESSAGE 


DAC* 

DBP 

/SAVE NUMBER 


lac 

( 140 

/ASCII FOR ZERO 


SZL 




LAC 

(152 

/ASCII FOR NUHGER 5 


TAD* 

DBP 

/ADD TO OTHER DB SETTING 


cllirar 


/SET UP AS CHAR* 4*5 IN ASCII WORDS 


OAC* 

DBP 

/SAVE NEW DB SETTING 

C0ND8 

ISZ 

COM 

/GET NEXT CHAR. 


LAC* 

COM 

/SET UP NEXT CHAR. 


SAD. 

(15 

/IS CHAR. A CARRIAGE RETURN ? 


JMP 

• *3 

/YES. EXIT 


SZA 


/NO. IS IT A ZERO 7 


DAC* 

CMULC 

/NO. SET UP MUL. CONSTANT CHANGE 


DZM* 

COM 

/CLEAR CHAR. 

BLK2 

JMS 

DBSUB 

/SET UP NEXT DB MESAGE 


LAC 

CN2 

/GET STORAGE ALREADY USED 


TAD 

CNI 

/ADD STORAGE SIZE OF LAST 


DAC 

CN2 

/ FILE AND SAVE 


tad 

TEM 1 

/ADD IT AGAIN AND CHECK- - 


TAD 

DSTOR 

/ WILL ANOTHER FILE OF THE SAME 


SMA 


/ LENGTH OVERFLOW THE STORAGE ALLOW? 

MDUM3 

JMP 

RPTI 

/YES. RESTART THE COUNTING 


ISZ 

NAME+1 

/NO. INCREMENT NAME 


ISZ 

MSG 1 ♦ 3 

/INCREMENT TELETYPE 


ISZ 

MSG 1 + 3 

/ MESSAGE TWICE 


JMP 

WRITE 

/RETURN TO NEW FILE 

STAGN 

LAC 

CNI 

/REMOVE THE SIZE OF 


TCA 


/ THE LAST FILE 


TAD 

CN2 

/ FROM THE 


DAC 

CN2 

/ STORAGE COUNTER 


LAC 

DSC 

/REMOVE THE AMOUNT OF 


TCA 


/ TIME USED 


TAD 

TC2 

/ BY THE PRECEEDING 


DAC 

TC2 

/ FILE 


JMP 

WRl TE 

/RECOLLECT FILE 

/DBDBDBDBDBDBDBDBDBDBDBDBDBDBDBDBDB DSDBDBDBD3DBDBDBDBDBDBDBDBDBDBDBDB 

/ 

SUBROUTINE TO CHANGE THE DB MESSAGE TO BE PRINTED OUT 

/ 

DBSUB 

0 




ISZ 

DBP 

/GET NEXT DB 


LAC* 

DBP 

/ MESSAGE 


SMA 


/IS IT THE END OF THE DB MESSAGES ? 


JMP 

• *4 

/NO. USE THIS MESSAGE 


LAC 

( DBO 

/YES. REPEAT THE 


DAC 

DBP 

/ FI R5T DB 


LAC* 

DBP 

/ SETTING GIVEN 


DAC 

MSQ 2* 1 l 

/INSERT MESSAGE 


JMP* 

DBSUB 

/RETURN 
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/ 

/ DBDBDBDBDBDBDBDBDBDBDBDBDBD9DSDB DBDBDBDBDS D8DBDBDBD8DB DBDBDBDBDS 
MSQ1 HSG2-MSG1/2#1 000 

0 

•ASCII <7><14><7><7>“r* 

• ASCII M PRESENT"* I 5* 

MSGS 2000 
0 

•ASCII "SET ATTENUATOR TO 00DB AND C.R."<15» 

/ 

/msstssssisi$$s$ssmsssssss$$$ssm$$$$sssms$ssis$mMS$si$si 

/ PACKING ROUTINE 

/ AN ID IS ASSIGNED AND 26 SAMPLES ARE PACKED TWO PER WORD 
/ A TOTAL OF 14 WORDS ARE PUT INTO A DEC TAPE BUFFER FOR EACH CALL 
/ HIGHEST ORDER BIT IS LOST AND NE3 . NUMBERS ARE SET TO EERO 
/ ALTERNATES BETWEEN TWO DEC TAPE BUFFERS- 


/ THE 

BUFFERS 

ARE LOCATED 1 

/ 

rsub 

0 



DAC 

SAV 


ISZ 

CNT 

MDUM1 

NOP 



DZM 

DUM2 


DZM 

DUM3 


LAW 

NPFC 


DAC 

CTT I 


I AC 
DAC 

CTTS 


LAW 

-1001 


SAD 

ER2 


JMP 

ERS1 

iNOISE 

ISZ 

IDCOU 


LAC 

IDCOU 


DAC* 

POIlNT 


ISZ 

POINT 

DATA 

LAW 

-NDPC 


JMS 

PAC 


LAC 

DUM3 


TAD 

DUMS 


DAC 

DUM5 


LAC 

DUM2 


TCA 

TAD 

DUM4 


SMA 

JMP 

.♦3 


LAC 

DUM2 


DAC 

DUMA 


ISZ 

err* 


JMP 

SKP 4 

COMP 

LAW 

MNCC 


DAC 

CTT2 


LAC 

UJM4 


TCA 

TAD 

MAX 4 


SPA 

JMP 

SKP 3 


LAC 

DUM4 


DAC 

MAX4 


LAC 

DUMS 


DAC 

SUM 4 

SKP 3 

DZM 

DUM4 


DZM 

DUMS 

SKP4 

LAC 

<JMP ONCE 


ISZ 

CNT2 


LAC 

(JMP RET 2 


DAC 

START 


LAC 

SAV 


• TIMER 

0* BEJ I N* 6 


.RLXI T 

RSUB 

PAC 

0 



DAC 

COUNTP 


LAC 

(BUF-I 


DAC 

0POINT 

FAKING 

JMS 

LIN 


LAC 

DUM1 


SWHA 

DAC* 

POINT 


JMS 

LIN 


LAC 

DUM1 


/SAVE AC 

/INCREMENT COUNTER 
/GUARDS AGAINST A -1 IN "CNT" 
/INITIALIZE THE VARIABLES USED 
/ TO DETERMINE THE NOISE 

/USE FIRST 5 NOISE SAMPLES TO SET 

/ MAX. ALLOWED NOISE FOR PROC. 

/SKIP NOISE AFTER THE 5TH 
/ NOISE SAMPLE 

/IS THERE A DISK TIMING 
/ ERROR? 

/YES* PRINT MESSAGE 
/NO* INCREMENT ID NUMBER 

/ASSIGN ID NUMBER 

/ TO BE PLACED ON STORAGE DEVICE 

/INCREMENT POINT TO NEXT LOCATION 
/COUNTER SO THAT 13 WORDS ARE PACKED 
/GO TO PACKING ROUTINE 
/GET SUM OF THE NEW 4 NOISE SAMPLES 
/ AND ADD TO THE 72 NOISE SAMPLE GROUP 
/USE NEW SUM 

/GET NEW MAX • FOR THE 4 NOISE SAMPLES 
/COMPARE NSW MAX. 

/ VI TH OLD MAX. 

/IS THE NEW MAX. LARGER? 

/NO# KEEP THE OLD MAX. 

/YES* REPLACE THE OLD MAX* 

/ WITH THE NEW MAX. 

/HAS 72 NO I SE SAMPLES BEEN COLLECTED? 

/NO* PREPARE TO COLLECT MORE SAMPLES 

/RESET COUNTER FOR 

/ NOISE DETERMINATION 

/YES# JET MAX. OF NEW 72 NOISE SAMP . 

/ AND COMPARE IT ■/! TH THE MAX. 

/ OF THE PREVIOUS SET 

/IS THE NEW OWE LESS THAN THE OLD OJT? 
/NO* PREPARE FOR ANOTHER COLLECTION 
/YES# REPLACE THE OLD MAX • WITH 
/ THE N£‘W MAX. 

/REPLACE THE OLD SUM WITH 
/ THE NEW SUM 

/RESET NOISE SA-MPLING 
/ LOCATIONS 

/PREPARE to REENTER R-T SUO. "3EGIN" 

/ BY INSERTING INTO LOC. START 

/ A JUMP STATEMENT DEPENDING ON 

/ CNT2 C WHICH TYPE OF SAMPLE) 

/ENTER R-T SUB. "BEGIN" 


/STORE WORD NUMBER BEING PACKED 

/GUF POINTER 

/.JO TO TABLE ROUTINE 

/ROTATE TO LEFT HALF 
/STORE 
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TAD* 

POINT 

/PACK INTO PREVIOUS ONE 

DAC* 

POINT 

/STORE IN BUFI 

ISZ 

POINT 

/MOVE POINTER UP ONE WORD 

ISZ 

COUNIP 

/ONE WORD HAS BEEN PACKED 

JMP 

PAKIN3 

/1 3 WORDS HAVE NOT BEEN PACKED 

LAC* 

POINT 

/1 3 WORDS HAVE BEEN PACKED 

ISZ 

CTT0 

/END OF BUFFER 7 

JMP* 

PAC 

/NO# CONTINUE COLLECTION 

LAW 

-DATSIR 

/YES# RESET THE BUFFER 

DAC 

CTT8 

/ COUNTER 

LAC 

BUF21 

/SET POINTER TO SECOND 

DAC 

POINT 

/ BUFFER— BUF2 

LAC 

BUF1 1 

/GET FIRST BUFFER ADDRESS 

ISZ 

CTT6 

/JUST FINISHED FILLING BUFI 7 

JMP 

.♦5 

/YES# PREPARE TO STORE BUFI 

DAC 

POINT 

/NO# SET POINTER TO BUFI 

LAW 

-2 

/RESET COUNTER WHICH DETERMINES 

DAC 

CTT6 

/ WHICH BUFFER TO TRANSFER 

LAC 

BUF21 

/PREPARE TO TRANSFER 

DAC 

TBUF 

/ BUF2 

DAC 

TRANF 

/SET TRANSFER FLAG 

JMP* 

PAC 

/CONTINUE COLLECTION 

.WRITE 

TTO# ASC#MSG 5# 0 

/TIMING ERROR 

LAC 

MDUM1 

/PUT CLOCK OPERATION BACK 

DAC 

ONCE 

/ IN PROGRAM 

ISZ 

CNT 

/PREPARE TO REJECT O-X PAIR 

LAC 

(JMP RDO 

/REJECT NEXT X-SAMPLE 

ISZ 

CNT2 

/WAS THE LAST AN 0- SAMPLE 7 

LAC 

(JMP ONCE 

/NO# REJECT THE FORMER 0- SAMPLE 

DAC 

START 

/YES# REJECT NEXT SAMPLE 

•TIMER 

0 # BEG 1 N# 6 

/RETURN TO COLLECTION 

•RLXIT 

2000 

0 

•ASCII 

RSUB 


"TIMING ERROR" 

A 

in 

V 


/SSS$5*SSSSSSSSm$$SS£i££S£S$S$SSS$$$£$S$$$s$$£SSS£SS$£$$$SS$$$$$S 

/ SUBROUTINE TRANSFERS 1 BLOCK OF DATA FROM A DESIGNATED BUFFER 
/ TO THE STORAGE DEVICE BEING USED* 

/ 


0 

LAC 

DAC 

•WRITE 

WAIT 

DEM 

JMS 

LAC 

DAC 

JMP* 


OUTPT# DUMP / 0 #252 

QUTPT 

TRANF 

CKCNT 

(NOP 

ONCE 

DTRANS 


/PASS NAME OF BUFFER 
/ TO •WRITE 


/CLEAR TRANSFER FLAG 
/CHECK CKCNT ROUTINE 
/PUT CLOCK BACK INTO OPERATION 


/ +♦ ♦ + ♦' + ♦+ + ♦♦♦♦ + + ++ + + ♦♦♦♦♦ «■♦♦ + -*■ + ♦ + + ♦♦♦♦ ♦♦♦ + ♦ + + ♦ + ♦♦4.4.4. + + + + + + + + + + 

/ SUBROUTINE TO CHECK WHICH DATA CONSOLE SWITCHES ARE SET* 

/ TOO LARGE SETTINGS AND TOO SMALL SETTINGS ARE GUARDED 
/ AGAINST. THE DEFAULT SETTING IS 2000, 


<17777 

(-7 

(2000 

(-10000 

(2000 

(10007 

TEMt# 


RESTAR 

TEM1 

CN2 

DSTOR 

RESTAR 

CKCNT 


/LOAD DATA SWITCHES FROM CONSOLE 
/PUT AC BIT 1 INTO LINK 
/TO DISALLOW SHARING 
/IS AC BIT 1 SET ? 

/YES# DO ALLOW 

/ SHARING 

/RELOAD DATA SWITCHES 

/IGNORE TOP 5 BITS 

/are THE SWITCHES SET TO 

/ LESS THAN SEVEN 7 

/YES# USE DEFAULT SETTING 

/NO# ARE THEY SET TO GREATER 

/ THAN 10000 (DISK OVERFLOW) ? 

/YES# USE THE DEFAULT SETTING 

/NO# RESET THE AC BACK 

/TEMPAR1LY STORE NUMBER 

/COMPARE THE DATA SWICHES 

/ TO THE ID NUMBER 

/IS THE ID LESS 7 

/ID GREATER THAN SWITCH SETTING 

/WILL THE SIZE OF 

/ THIS FILE 

/ OVERFLOW 

/ THE DISK ? 

/YES# CLOSE FILE 

/ID LESS THAN SWITCH SETTING 



/ 0000000 0000000 000 0 00000000000040 0000400 044 44040440044440 04 00 §4044 

/ ROUTINE TO CLOSE FILE AND SET CJP PARAMETERS FOR 

/ COLLECTION AND PROCESSING 
/ 


hestak 

LAC 

TftANF 

/COMPLETE TRANSFER IF NECESSARY 

MDUM2 

5ZA 




JMS 

DTHANS 


C 6 

•WRITE 

OUTPT * DUMP# 

SEVN* 2 /WRITE END OF FILE ID 

C2 

« WAI T 

OUTPT 


C 7 

.CLOSE 

OUTPT 



LAC 

MDUM1 

/INITIALIZE THE LOCATIONS 


DAC 

START 

/ M START M FOR COLLECTION 


DAC 

BEPL 

/ AND M R£PL" FOR PROCESSING 


LAC 

SUM A 

/GIVE THE SUMMATION OF THE LOWEST 


DAC* 

SUMS 

/ NOISE TO THE PROCESSING 


DZM* 

CHG 

/RESET DB CHANGER 


LAC 

MDUM2 

/IGNORE THE FIRST 


DAC 

UPl 

/ FILE ONLY 


JMP 

START 

/GO TO START 


/ 


/&&&&&&&&&&&&&&& 


/<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
/ CLOCK INTERRUPT ROUTINE FOR AUTOMATIC O-SAMPLE START 


/ 

CHECK 

0 




DAC 

SAV 

/SAVE AC 


DZM 

BEGIN 

/ZERO R-T SUB* TO AVOID POSS. ERROR 


LAC 

CNT 

/EXAMINE COUNTER 


RAH 


/LOWEST BIT OF CNT IN L 


SZL 




JMP 

BKUP 

/L= 1 -CNT ODD-HALF FRAME DURING 9/60 SEC 


JMP 

NORM 

/L=0 -CN T EVEN-FULL FRAME DURING 9/60 SEC 

8KJP 

DZM 

CNT 

/CLEAR HALF FRAME THAT WAS TAKEN t 


RAL 


/RESTOR LINK 


LAW 

-1 

/RESET IDCOU BACK ONE FRAME 


TAD 

IDCOU 



DAC 

I DC O'J 



LAW 

-1 

/RESET BUFFER COUNTER 


TAD 

CTT8 

/ BACK ONE 


DAC 

CTT8 

/ FRAME 


LAW 

- DPR 

/RESET POINT BACK HALF FRAME 


TAD 

POINT 



DAC 

POINT 



LAC 

MDiJMl 

/PUT CLOCK BACK INTO OPERATION 


JMP 

EXT 1 

/PREPARE TO EXIT 

NORM 

RAL 


/RESTORE THE LINK 


LAC 

(JMP PKO 

/STOP CLOCK FROM OPERATION 

EXT 1 

DAC 

ONCE 



LAC 

( JMP RDO 

/PREPARE TO COLLECT 


DAC 

START 

/ AN 0- FRAME 


LAC 

SAV 

/RESTORE AC 


« TIMER 

0*3E3 IN*6 

/HETUilN TO POINT AT WHICH INTERRUPT 

/ 

.RLXIT 

CHECK 

/ OCCURED 


/>>>>>>>>>>> >> >>> >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>>> 

/ m ! m m i m m m m ■ ! * ! ! n 1 1 r m ( ? h r t ! ! m f m f t ? m r f ! ! t m r m m 

/ SUBROUTINE TO DO A TABLE LOOKUP FOR DATA LINEARIZATION* 

/ ALSO SAVES SUM AND MAXIMUM OF EACH SET OF S NOISE SAMPLES* 

/ SETS NEGATIVE NUMBERS TO ZERO* AND IF NECESSARY JUMPS AROUND 

/ "SKAR M NUMBER OF DATA NUMBERS BETWEEN THE NOISE AND DATA POINTS* 

/ 

lin y 


Eiil 


ISZ 

BPOINT 

LAC* 

BPOINT 

AND 

<1777 

TAD 

(-I0J3 

SMA 


JMP 

£R1 

TAD 

CST1+1300 

OAC 

DJM 1 

LAC* 

DIJM1 

SKP 


CLAICLL 


DAC 

DUM1 

ISZ 

CTT5 

JMP 

* + 5 

LAC 

CSKAR 

TAD 

BPOINT 

DAC 

BPOINT 

LAC 

DUM 1 

ISZ 

cm 

JMP 

DOMOS 

LA 7 

- 1 

DAC 

CTT 1 


/DATA STARTS AT 0-J F 
/GET INPJT DATA WORD 
/MASK ANY EXTRA BITS 
/CHECK FOR NEG. #♦ S 

/N EG 0 FOUND 
/LOCATE ? IN TABLE 
/GET ADDRESS OF NUMBER 
/LOAD LINEARIZED # INTO AC 

/SET NEG 0 TO ZERO 

/STORE LINEARIZED 0 

/IS THIS IKE FIFTH NOISE SAMPLE ? 

/NO* SKIP AROUND CODE 

/YES* SKIP AROUND "SKAR M DATA 

/ OF THE FRAME BEING 

/ COLLECTED 

/RESTORE THE LIN* DATA $ 

/IS IT A NOISE SAMPLE ? 

/TES* P ROCK SS NOISE SAMPLE 
/RESET THE 7WO COUNTERS 
/ FOR THE DATA OF THE 
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JMP* 

LIN 

DONOS 

TAD 

DUMB 


DAC 

DUM3 


LAC 

TCA 

DUM1 


TAD 

SMA 

DUM2 


JMP* 

LIN 


LAC 

DUMI 


DAC 

DUM2 


JMP* 

LIN 


/ RETURN 

/ADD TO THE OTHER 4 

/ NOISE SAMPLES 

/IS THE NOISE SAMPLE 
/ GREATER THAN THE 
/ OTHER 4 
/NO# RETURN 
/YES* SET THIS SAMPLE 
/ AS THE MAX. 


/ 

/?????7?777??7???7??7??7????????????7?????????????7????77????????77 
/&&&& &£&&£&&& A &&&&£& A &&&£ AS 
/ ROUTINE TO REINITIALIZE VARIABLES AND IF NEEDED TO 

/ SWITCH STORAGE DEVICES* 


/ 

RPT1 


NDK2 


NDK21 


A FI 
EXTIP 


EDC 


LAC 

AFl 

/RESET FILE 

DAC 

NAME* 1 

/ 

NAME 

DZM 

TEM1 

/ZERO LOG. 

DZM 

CN2 

/RES 

ET DEVICE STORAGE COUNTER 

JMP 

NDK2I 

/DO 

NOT SWITCH STORAGE DEVICE 

LAC 

(OUTPT 

/CHANGE TO SECONQ STORAGE 

ISZ 

CTT3 

/ 

DEVICE BY CHANGING THE 

LAC 

(0UTPT2 

/ 

.DAT SLOTS IN THE COMMANDS: 

DAC 

C 3 

/ 

• EN TER 

DAC 

C2 

/ 

.WAIT 

DAC 

C 5 

/ 

.WAIT 

DAC 

C7 

/ 

.CLOSE 

TAD 

(1000 



DAC 

WRITE 

/ 

♦ DLETE 

DAC 

Cl 

/ 

♦ INI T 

TAD 

(2000 



DAC 

UP 

/ 

. FSTAT 

TAD 

C 1000 



DAC 

C4 

/ 

• WRI TE 

DAC 

C6 - 

/ 

•WRITE 

LAC 

nasal 

/RESET TELETYPE 

DAC 

MSG 1*3 

/ 

MESSAGE 

LAC 

MDUM4 

/IGNORE THE FIRST FILE ON 

DAC 

UPI 

/ 

NEXT STORGE DEVICE 

LAC 

C TT 3 

/CHECK DEVICE SWITCHING CONTROLLER 

SMA 


/IS 

IT STILL NEG.7 

LAW 

-2 

/YES* RESET IT 

DAC 

CTT3 

/NO# 

LEAVE IT ALONE 

JMP 

Cl 

/NO* 

GO TO INIT. 

. SIXBT 

"AF1 M 

/REINITIALIZES "NAME" AND "FIL" 

DZM 

TC2 

/PREPARE TO COLLECT A NEW SET OF DATA 

LAC 

TIME 

/GET TIME FOR LAST FILE 

TAD 

(5 

/ROUND OFF TO THE NEAREST MINUTE 

DAC 

TIMR 

/STORE FOR FORTRAN PROGRAM 

LAC 

( JMP HP T 1 

/PREPARE TO START 

DAC 

START 

/ 

COLLECTION AT FILE 1 

ISZ 

NAME* l 

/TELLS PROC. TO PROCESS LAST FILE 

.write 

TTO# ASC* EDO# 0 

/END OF COLLECTION 

.WAIT 

TTO 



DZM 

TCI 

/INI T. COLL. TIME COUNTER 

LAW 

-1 

/TELLS PROCESSING THAT COLLECTION 

DAC 

CTT1 1 

/ 

IS FINISHED COLLECTING 

• RLXl T 

BEGIN 



3000 




a 

.ASCII 

<1 I><1 l>< ll><ll><tl><7><7><7><7»<lS> 


/ 

/ AAAAAAA&AAAAAAAAAAAAAAAAAAAAAAAAAAA&AA&AAAA&AAAAAAAAAAAAA&AAAAAAAA 
. EJECT 

/ + + + + + +++ + + + + + + + + + + * + +**♦*♦+***♦♦♦♦*♦+ ♦♦♦+*♦♦♦ ********** ********** 

/ BFKM15 V3A SERVICE ROUTINES FOR THE HP S610A A TO D 

/ CONVERTER. these routines permit input of any specified 
/ NUMBER of SAMPLES INTO A CORE BUFFER. INPUT MAY 3£ OVER- 
/ LAPPED WITH PROGRAM EXCUTIOM# AND CONTROL MAY BE RELINQUISHED 
/ TO LOWER PRIORITY PROGRAMS 'WHILE DATA TRANSFER TAKES PLACE. 

/ MACRO- l 5 CALLING SEQUENCE: 

/ JMS ADREAD 

/ NUMBER OF SAMPLES REQUIRED 

/ BUFFER ADDRESS 

/ COMPLETION FLAG ADDRESS 

/ REAL-TIME SUBROUTINE ADDRESS# PRIORITY LEVEL IN BITS 0-3 

/ (EXAMPLE: 500000+HTSUBA ) 

/ (RETURNS HERE IMMEDIATELY) 

/ IF THE 4TH WORD AFTER THE JM5 IS 0# NO REAL-TIME SUBROUTINE 

/ WILL BE ACTIVATED. NOTE! THE PRIORITY CODE FOR MAINSTREAM IS 1 
✓ THE COMPLETION FLAG IS CLEARED BY THE CALL TO ADREAD* 

/ AND SET TO *1 FOR NORMAL COMPLETION OR -1001 IF A DATA 
/ TIMING ERROR OCCURS. 

/ 
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ADWCR*26 
ADCAR-ADWCR+1 
♦SCOM-100 
A DWI -703724 
ADSO*70370i 
ADST-703721 
ADCO-703704 
ADCT-70 3744 
/ 

/ ENTRY POINT FOR A 
/ 

A DREAD 0 


/A-D WORD COUNT 

/AND CURRENT ADDRESS REGISTERS 
/MONITOR’S COMMUNICATION AREA 
/A-D CONVERTER WRITE INITIALIZE 
/SKIP ON WORD COUNT OVERFLOW 
/SKIP CN DATA TIMING ERROR 
/CLEAR OVERFLOW FLAG 
/CLEAR TIMING FLAG 


-D INTERFACE INITIALIZATION 



JMP 

TCA 

INSET 


DAC* 

(ADWCR) 


ISZ 

ADREAD 


LAW 

-1 


TAD* 

ADREAD 


DAC* 

( ADC AR ) 


ISZ 

ADREAD 


LAC* 

ADREAD 


DAC 

I NFL AG 


DZM* 

I NFLAG 


ISZ 

ADREAD 


LAC* 

ADREAD 


DAC 

INSUB 


ISZ 

ADWI 

ADREAD 

/ 

JMP* 

ADREAD 

/ the 

FOLLOWING CODE IS 

INSET 

LAC* 

( • SCQM+55 ) 

ADSVA 

DAC 

• 

SAV 

LAC* 

( .SC0M+5I ) 

HEALTP 

DAC 

. 


LAC 

ISA 

(400010 


JMS* 

ADSO 

ADINT 

DBK 

ADSVA 


LAC 

( LAC* ADfr 


DAC 

ADREAD* 1 


JMP 

ADREAD* 1 


/REPLACED BY "LAC* 
/SET WORD COUNT 


A DREAD* 1 


/BUFFER ADDRESS -1 
/ TO CURRENT ADDRESS REG. 

/GET FLAG ADDRESS 

/CLEAR FLAG 

/GET REAL-TIME SUBROUTINE ADDRESS 

/POINT TO RETURN LOCATION 
/INITIALIZE INTERFACE 
/RETURN 

(TED ONLY ONCE 

/GET ENTRY POINT ADDERSS OF .SETUP 

/ENTRY POINT OF REAL TP 

/RAISE THE API 
/ LEVEL 

/CALL .SETUP TO CONNECT 
/ ADI NT TO 

/ THE API 

/DE3REAK FROM API LEVEL 

/MODIFY INSTRUCTION 
/ AND JUMP TO IT 


/INTERRUPT SERVICE ROUTINE. EXECUTED IMMEDIATELY AFTER COMP! FT, n„ 
' OF DATA TRANSFER. DETERMINES STATUS OF A-D INTERFACE SFT^ 

X COMPLETION FLAG AND ACTIVATES REAL-TIME SUBROUTINE? 

/ RUNS AT API LEVEL 0. ^ 


/ 

/ 

ADINT 


0 

DBA 

DAC ADS\ 

ADST 

SKPfCLAJ I AC 


-100 I 
INFLAG 


C .SCOM+103 


JMS* REAL TP 

ADXIT LAC (404000 

ISA 

LAC ADSVA 

DBK 

JMP* ADINT 

/ 

/♦♦+ +♦♦+♦♦♦+*+♦♦ +♦ + + 4 . +++++++++ 

.EJECT 


/PAGE ADDRESSING MODE 
/SAVE AC 
/TIMING ERROR? 

/NO# + 1 TO AC 
/YES# ERROR CODE 
/SET FLAG 
/CLEAR 

/ INTERFACE FLAGS 
/HA I SE TO API 
/ level 0 OR I 

/REAL-TIME SUBROUTINE ADDRESS 

/BYPASS MONITOR CALLS IF ZERO 

/ACTIVATE REAL-TIME SUBROUTINE 
/REQUEST AN API INTERRUPT 
/ AT SOFTWARE LEVEL H 

/RESTORE AC 

/SET TO LEAVE HARDWARE API LEVEL 
/RETURN TO INTERRUPTED PROGRAM 

4 4 * *'*’ ** + 4 * +++++ + ♦♦ + + ♦ + *♦♦ +++++++ + + + 


. EJECT ' 

/xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

PROCESSING S MACRO PROGRAMS. THEY INITIALIZE THE 
' DEVICE, WIT FOR FILE TO BE COLLECTED, CHECH FOR WWN?™ 

X COLECTION STOPAGE. SWITCH STORAGE DEVICES IF NECE^W, ?ELL 
/ SSL™* PROGRAMS HAVE REACHED THE END OF RUN, GlvIS TO 
/ BACKGROUND# AND READS IN DATA AND UNPACKS IT. 
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/ 

CHNG 

0 




JMS* 

• DA 

/LOAD PARAM. ADDR. IN SUMS AND CHG 


JW 

.♦A 

/SKIP OVER PARAM# LIST 

SUMS 

0 



CHG 

0 



CMULC 

0 


/LOC. TO CHANGE MUL. CONSTANT 


I sz 

CTT10 

/RESTART THE DEVICE ? 


SKP 


/NO 


JMP 

AGN 

/YES* RESET PARAMETERS TO BEGIN* OF DEV' 

DU M7 

I SZ 

FJL+1 

/CHECK IF COLLECTION IS 


JMS 

CON TR 

/ FINISHED COLLECTING NEW FILE 


LAC 

FILM 

/IS the collection 


CMA 


/ RECOLLECTING 


TAD 

NAME+1 

/ THE FIRST FILE 


SPA 


/ AGAIN ? 


LAW 

-l 

/NO* SET THE COUNTER TO Jl*IP AROUND 


DAC 

CTTI0 

/YES* SET COUNTER TO A POS. # 


LAW 

-l 

/PREPARE TO READ 


DAC 

CTT9 

/ DATA 


LAC 

BUF3 1 

/RESET BUF3 POINTER WITH 


DAC 

POINTS 

/ THE ADDR. OF BUF3 


LAC 

CISZ SWITC 

/PREPARE TO READ TWO 


DAC 

LCA 

/ DUMMY BLOCKS 

C13 

• INIT 

DATIN*IN#0 


HT 

. FSTAT 

DATIN# FIL 

/CHECK FOR PRESENTS OF NEW FILE 


SNA 




JMP 

ERR3 

/FILE NOT PRESENT 

Cll 

.SEEK 

DATIN# FIL 

/PREPARE TO READ 


JMP* 

CHNG 

/RETURN TO FORTRAN PROGRAM 

FIL 

•SIXBT 

"DATAF0DAT" 

/NAME OF DATA FILE USED BY PROC. 

EKR3 

.WRITE 

TTO# A SC* MSG 6#0 

/FILE NOT FOUND 


JMP 

CMULC+ I 

/ERASED FILE ? LOOK FOR NEXT FILE 

MS36 

2000 




0 

.ASCII 

"FILE NOT F0UND"<I5> 

AGN 

LAC 

AF1 

/INITIALIZE DATA* FILE 


DAC 

FIL+1 

/ NAME 

N DK3 

SKP 


/DETERMINES 9 OF STORAGE DEVICES 


JMS 

SWlf 

/CHANGE STORAGE DEVICES 


JMP 

DUM7 + 1 

/NO* CONTINUE PROCESSING 

ENDPK 

.WRITE 

T TO* A SC# END# 0 

/END OF PROC. MESSAGE 


.WAIT 

TTO 



JMS 

REA DM 

/WAIT FOR RESPONSE AND PUT 

ANSE 

0 


/ THE ADDR. OF IT HERE 


LAC* 

anse 



SAD 

(120 

/IS THE RESPONSE A ,, P M ? 


JMP 

RSTR 

/YES* SET UP TO RESTART EVERYTHING 


.IDLE 


/GIVE COMPLETE CONTROL TO BG* 

RSTR 

.clear 

DATIN 

/CLEAR COLLECTION DEVICE 


•WRITE 

TTO# ASC*MSG2*0 

/DB MESSAGE 


• WAIT 

TTO 



JMS 

DBSUB 

/SET UP NEXT DR SETTING 


JMS 

REA DM 

/WAIT FOR REPLY 

dmrply 0 


/THE REPLY 


DZM 

MONTH 

/ZERO BINARY TIMER 


.TIMER 

9* BE j IN# 6 

/RESTART COLLECTION 


JMP 

AGN 

/RESTART PROC. 

/ subroutine 

USED TO WAIT FOR FILE TO BE COLLECTED AND STORED 

/ WHILE GIVING TINE TO BACKGROUND 

CONTR 

0 



WT2 

LAC 

NAME+1 

✓COL ECTI ON* S FILE NAME 


SAD 

FIL+1 

/IS PROC.'S FILE NAME THE SAME ? 

MDUM4 

SKP 


/YES# GO TO WAIT 


JMP* 

CONTR 

/NO, RETURN TO PROCESS FILE 


ISZ 

CTT 1 i 

/END OF THE 


SKP 


/ RUN 7 


JMP 

ENUtPii 

/YES*EXI T 


• TIMER 

30*WAI T1 * 0 

/YES# RELINQUISH TIME TO 93. 


LAC 

SAVAC 

/RESTORE AC 


.IDLE 


/WAIT FOR CLOCK INTERRUPT 

/HEAL 

TIME SJ3.--JSED TO ALLOW 

TIME FOR BACKGROUND 

W A I T 1 

0 




DAC 

SAVAC 

/SAVE AC 


DZM 

AI T 1 

/ZERO H-T SUB » ENTRY PI. TO ALLO; REENTK 


JMS 

CKCOL 

/CHECK FOR COLLECTION STOPPAGE 


JMP 

wra 

/CHECK FOR END OF COLL. 

/ 

SUBROUTINE TO CHECK FOR UNWANTED COLL. STOP AT E 

CKCOL 

0 




LAW 

-2 



TAD 

CNT 

/HAS COLL. ENDED ALL COLLECTING 


SPA 


/ FOR TODAY? 


JMP* 

CKCOL 

/YES# RETURN 


LAC 

TST1 

/NO* HAS COLLECTION STOPPED 



REPL 


TCA 

TAD 

SZA 

JMP 

LAC 

DAC 

ISZ 

• WRITE 
•TIMER 


IDCOJ 

SET! 

MDUM1 

ONCE 

CTT7 

TTO# ASC# STPM#0 
0# CHECK# 5 


/ READING 

/ IN DATA? 

/< REPLACED BY "NOP” WHEN COLL. 
/NO* RESET TESTER AND RETURN 
/YES# FREE TIMER 
/ OPTION IN COLL. 

/TELLS IF STOPPAGE OCCURED 
/RING BELL 
/RESTART COLLECTION 


IS DONE) 


SET 1 


/ 

SUIT 


/ 

PP7 


LAC 

IDCOU 

/SET 

ID # INTO 

DAC 

TST I 

/ 

TESTER FOR STOPPAGE 

JMP* 

CKCOL 

/CHECK FOR COLL. TO BE FINISH 

U TO 
0 

LAC 

CHANGE THE 

.DAT SLOT #' 

S TO CHANGE STORAGE DEVICES 

(DATIN 

/SWITCH STORAGE DEVICES 

ISZ 

CTT4 

/ 

BY CHANGING THE 

LAC 

(DATIN2 

/ 

.DAT SLOT IN COMMANDS! 

DAC 

Cll 

/ 

• SEEK 

DAC 

Cl 3 

/ 

• INI T 

DAC 

Cl 4 

/ 

• WAIT 

DAC 

CIS 

/ 

•CLOSE 

DAC 

RSTR 

/ 

•CLEAR 

TAD 

(3000 



DAC 

RT 

/ 

• FSTAT 

TAD 

(1030 



DAC 

LBA 

/ 

• READ 

LAC 

CTT4 

/IS 

DEVICE ON -DAT SLOT "DATI N" 

SMA 


/ 

TO BE PROCESSED? 

LAW 

-2 

/YES# RESET DEVICE CONTROLLER 

DAC 

CTT4 

/MAKE ANY CHANGE IN CONTROLLER 

JMP* 

SWIT 



0 

LAS 


/GET 

THE CONSOLE DATA SWITCH 

AND 

<43030 

/ 

NUMBER 3 

SNA 


/IS 

IT A 1 ? 

JMS* 

WRPP 

/NO, 

PRINT DATA OUT ON PAPER TAPE 

JMP* 

PP7 

/RETURN TO PROC 


/xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

. EJECT 

/iXZZX%%XZXZXZZZ*XXX*XZXXXlXXXt%ZXZZXZZ%Z*ZXZX%ZiXZ%%X%XXXXZZXi%XZ*lZX 
/ •READ, DUMP MODE FROM DEC TAP E ON A VARIABLE ♦ DAT SLOT 
/FILLS 252 DEC WORD BUFFER AMD OUTPUTS 26 
/WORDS TO ARRAY I DAT EVERY TIME CALLED. 

/THESE ARE UNPACKED FROM 16 WORDS OF THE BUFFER. 


/ I DAT* 
/ 

/ 

/ N EG F : 
/ 

DUMPT 


A2 

FLAG 


WORD 1 
WORD 2-6 


I . D. 0 

NOISE SAMPLES 


WORD 

7-27 DATA 




SET IF A NEGATIVE NUMBER WAS IN 

THE DATA 

0 

JMS* 

.DA 

/PICKUP ADDR OF 

AD DH 

JMP 

0 

0 

• + 3 

/OF ARRAY 



/SET ON NEG # 


LAC* 

A2 

/GET ADDR. 


DAC 

A2 

/ OF ARRAY 


LAW 

-NDPC 

/SET COUNTER OF 

DATA TO BE 

DAC 

COUNT 

/ PROCESSED 


ISZ 

CTT9 

/GET POINTER 


JMP 

LBB 

/NO# CONTINUE WITH PRESENT 

LAW 

-DATSTR 



DAC 

CTT9 




LBA .READ DATI N# DUMP# BUF3# DTBLK 

C 1 4i .WAIT DAT IN 

LCA ISZ 5WITC 

JMP LBA 

LAW -3 

DAC SWITC 

LAC < JMP LCB 

DAC LCA 

LCB LAC BUF31 

DAC P0INT2 

LBB LAC* P0INT2 

SAD SEVN 

JMP ENF 

LAC* POINTS 

DAC* A2 

ISZ A3 

LOOP ISZ POINT2 

LAW -2 

DAC TC10 

LAC* PQINT2 


/GET 1 BLOCK OF DATA 
/INITIALLY READ TWO DUMMY BLOCKS 


/RESET CONTROL TO READ 
/ TWO DUMMY BLOCKS 
/HEAD ONE BLOCK OF DATA 
/ AT A TIME 
/GET ADDRESS OF BUF3 
/ POINTS TO BUF3 

/GET THE ID (FIRST WORD IN DATA SET) 
/END OF FILE ID? 

/YES, RESET PARAM* S AND CLOSE FILE 

/GET ID AND PUT 

/ INTO THE FORTRAN ARRAY 

/GO TO NEXT A DDR. OF THE ARRAY 

/GO TO NEXT DATA WORD 

/PREPARE TO UNPAC 

/ TWO DATA WORDS 

/GET DATA WORDS FROM BUF3 




SWHA 


/FIRST WORD IN LEFT HALF 

UNPLP 

AND 

(777 

/SAVE ONE DATA WORD 


SNA 


/CHECK FOR NEG. NIMBKR 


isz* 

FLAG 

/SET IF NEG • NUMBER POUJD 


DAC* 

A2 

/LOAD # INTO FORTRAN ARRAY 


LAC* 

POINTS 

/GET DATA WORD AGAIN 


ISZ 

A2 

/GO TO NEXT LOC. IN ARRAY 


ISZ 

TCI 0 

/UNPAGED TWO WORDS? 


JMP 

UNPLP 

/NO* LOOP AROUND 


ISZ 

COUNT 

/YES* HAS 34 DATA WUHU* Wfcfcw uwrnun^u? 


JMP 

LOOP 

/NO* REPEAT UNPACKING PROCESS 

0UT3 

ISZ 

POINTS 

/YES* GO TO NEXT ID 


JMP* 

DUMPT 

/RETURN 

suite 

777775 



POINTS 

• OSA 

9UF3 


/- * 

/ END 

OF FILE ROUTINE 


ENF 

LAC 

SEVN 

/SET LAST ID TO 


DAC* 

A2 

/ 130050 DECIMAL 

C15 

•CLOSE 

DATIN 

/CLOSE FILE 

RT2 

JMP* 

DUMPT 

/RETURN TO PROC. PROGRAM 


/mxxiuxuimutsixtxnxxxnxnxxHxuumxuutumuxmxxnxz 

• END INTIM 


CCCCCCCCCCCCCCCCCCCCC--CQNTL~CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C PROGRAM SETS UP CALIBRATION AND THE HEADING FOR THE PRINT 

C OUT OF THE PROCESSED DATA IN PROC. THE PROGRAM CALLSl 
C TBFQRL--CALIBRATION PROGRAM 

C DLOGF-- - INITIAL MACRO PROGRAM 

c cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc ccccccc 
c 

SUBROUTINE CONTLUSURP) 

INTEGER DB< 4) # DBC 
REAL DATE( 2) * REASC 5) 

COMMON /STAT/ DB* DATE, REAS. DBC, DBS, NC4 
DATA NR/ 120784/ 

REWIND 4 

C INITIALIZE THE DB SETTINGS USED 
DBC >3 
DBC 2) 8 10 
DB(3>=25 

X F< I SURP • NE.0 )G 0 TO 5 
C GIVE PRECALIBRATION SETUP 
URITEC 6*100) 

FORMAT ( 46H TURN OFF PULSER AND ENCODE PULSE POWER SUPPLY/) 

UP CALIBRATION AND LINEARIZATION TABLE 
CALL TBFORL 
FOR AND GET THE DATE 
VRIT£<6*10I>NR 
FORMAT < 5H DATE* Ad ) 

READC 4* 20 1 ) DATE 
FORMAT ( 2A5 ) 

FOR AND GET THE REASON FOR THE RUN 
VK1 TEC 6* 1 02 >NK 

FORMA T< 1 6H REASON FOR DATA*A2> 

READC 4* 202 ) REAS 
FORMAT <5A 5) 

C PREPARE FOR COLLECTION AND PROCESSING 
CALL DLOGFCDB* DBC* DATE) 

RETURN 

END 


1 00 

C SET 

C ASK 
5 

101 

201 

C ASK 
132 
202 


CTTTTTTTTTTTTTTTTTTTTTTTTTT— •TBFORL- — TTTTTTTTTTTTTTTTTTTTTTTTTTTTTTT 
C THIS PROGRAM SETS UP A LINEAR APPROXIMATION TO THE INPUT OF 

C THE RECEIVER VERSUS THE OUTPUT OF THE A/D CONVERTER AS READ BY THE 
C COMPUTER. FROM THl SLI NEAR APPROXIMATION* A TABLE IS FORMED BY 
C USING INPUTS FROM .5 TO 511.5 INCREMENTED BY t* THE OFFSET OF *5 
C IS USED FOR BETTER ACCURACY IN THE ROUNDOFF ERROR OF THE A/D CONVER- 
C TER* ALL THE VALUES OF ZERO IN THE TABLE ARE CHANGED TO 1 SINCE 
C 0 IS USED TO DESIGNATE NEGATIVE NUMBERS DURING THE COLLECTION. 

C THE PROGRAM CALLS THE SUBROUTINES! 

C RADC---READS #'S FROM THE A/D CONVERTER (MACRO) 

C LI NAP -•CONVERTS INPUTS TO OUTPUTS FOR THE FORMATION 

c OF THE LINEARIZATION TABLE 

c TTM----WRITES LIN. TABLE ONTO A STORAGE DEVICE 

CTTTTTTTTT TTTTTTTTTTTTTTTTTT TTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTT 
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c 


c 

7 

105 


304 


C 


C 


C 


C 


20 

C 

1 1 

C 

100 

200 

C 


SUBROUTINE TBFORL 

INTEGER START! 512)* DLPO* STAT* I AS< 44 > 

REAL Kl 

COMMON /TA/ S(43>*TU(44>* TU0<44> 

DATA C#CS»CN*NR/1HC# LHS# 1HN* 1 20704/ 

DETERMINE IF CALIBRATION IS NEEDED AND WHICH PRINTOUT TO USE 
WRITE! 6* 105>NR 

FORMAT! 1 8H WHICH CALI BRAT ION/ 2 4H !S- SHORT* R- REGULAR# OR 
137H C-COMPLETE PRINT OUT OR N-NO CALIB.)#A2) 

READ! 4* 204 > CAL 
FORMAT! A1 ) 

IER R«0 
DLPO-0 

SET UP THE WANTED CALIBRATION 
I FC CAL • EQ . CN > RETURN 
IF! CAL » EQ « C S) DLPO«- 1 
IF! CAL • EQ « C ) DLPO* 1 

NUMBER OF A/D NUMBERS TO READ PER DB SETTI NG-NAVt *NAV 
NAVI * 5 
NAV-500 
AV=NAV1*NAV 
ICT=*0 

NUMBER OF ATTENUATOR 5ETTINGS-NT1 
N T 1 = 42 
NT2-NT1-I 
NT3-NT1-2 

INPUT SIGNALS--FHOM 5DB (OF 1000 > TO INFINITY 
TU!NTI ) = 562 « 34 
TU(NT2) = 501 • 19 
TU(NT3>=446.68 
TU!NT1-3)»396.1 I 
TU(NTl-4>=354.82 
TUCNT1-5>=316.23 
TU(NTl-6>=281 .84 
TU (NT l “7 ) =251 . 19 
TU(NTl-8>=223*87 
TU(NTI-9)=199. 53 
TU(NTl-10 >=177.83 
TIKNTl-il > = 158.49 


TU(NT1-12>=141*25 
TU(NTl^l 3>*125.90 
TU(NT1-'14) = 1 1 2.20 
TU(NT1-15>=I00.00 
TU!NT1-16>= 89.13 
TU!NT1-17>= 79.43 
TU!NTl-i8>= 70.80 
TUCNT1 -I 9 > = 63.10 
TU (NT 1 -20 ) = 56.24 
TU (NT 1 -2 1 > = 50.12 
TU!NT1 -22 )= 44*67 
TU (NT 1 -2 3 >= 39.81 
TU (NT1 -24 )= 35.48 
TU!Ni'l-25>= 31.62 
TU(NT 1 -26 > = 28.18 
TU (NT 1 -27 ) = 25.12 
TiJCNT 1 -28 >= 22.39 
TU(NT 1 -29 > = 19*95 
TU(NT 1 “30 >= 1 7 • 7 B 
TU (NT 1-31 >= 15.85 
TU!NTl-32>* 14.13 
TU CNT1 -33 >= 12.59 
TJ(NTl-34)= 11.22 
TU!NTl-35>= 10.00 
TU CNT1 -36> = 8.91 

TU (NT 1 -37 ) = 7.94 

TU (NT 1 *38 ) * 1 *08 

TU(NTl-39)= 6*31 

TU (NT! -43 ) = 5.62 

T U ! 1 > = 0*00 

TU0(NT1 >=512. 

SET UP MESSAGES FOR TELLING .-/HICH ATTENUATOR SETTING TO DO 
DO II 1 = 1 # N T 3 
IAS C I ) = I + 5 
I AS ! NT2 > =99 

DO LOOP TO INPUT ALL THE OUTPUTED SIGNALS 
DO 2 1 1 = 1 * N T2 
1 2= NT 1 - 1 1 
DUM = 0 

rfRIIE (6*100) I AS ( 1 1 )j NR 

FORMAT (8H SET TO *I2*23HDB ATTENUATION AND C.R*#A2> 
HEAD! 4* 200 ) F 
FORMAT (F5. 1 > 


IF THE INPUTED NUMBER IS TWO DIGITS RESTART THE 
I F(F*GT.10. >30 TO 23 


SETTINGS 
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C ISSUE AN A/D CONVERTER HEAD "NAVI ” TIMES 
DO 10 J1-1,NAV1 

C READ "NAV" NUMBERS FROM THE A/D CONVERTER 
CALL RADC (START, NAV, STAT) 

400 I F( STAT*EQ*0 >G0 TO 400 

DO 1 J* 1 # NAV 

IF*START<J>*GT*511 > STAR1(U>=0 
C STORE THE INPUTED NUMBERS IN A REAL VARIABLE 
1 DUM* DUM*FL0AT<5TART*U) > 

10 CONTINUE 

AVERAGE THE OUTPUTED NUMBERS AND GO TO THE NEXT SETTING 
TUO< I 2 ) *DUM/AV 
DO 3 11*1, NT3 

12 = 11+1 

ICT=1CT+1 
1 3 a NTi •! 1 

SET UP THE SLOPES OF EACH LINE SEGMENT APPROXIMATION 
SCIl >=*TJ*I2)-TU<I1 ) )/<TUO*I2)-TUO(Il )> 

IF<ICr*EQ*3)30 TO 3 

WHITE OUT STRAIGHT LINE APPROX1MAT1 ON TABLE 

WHIT EC 6, 102) 1 1 , S< 1 1 >, 1 1 » TJC I 1 > , I 1 , TUO ( I 1),IAS(I3> 

I CT**l 

I F < DLP 0 « G T * - 1 >ICT=0 

POSSIBLE ERROR CONDITIONS FOR THE APPROXIMATION JUST FORMEO 
IF* SC II )#LT* ,3. OR* SCI 1)*GT*7* > IEHH*I ERR+ 1 

102 FORMAT <4H S* , 1 2, 2H) = , F6 * 3, 5X, 3HTUC , I 2, 2H > «, F8* 3, 5X, 

1 4HT UOC, 1 2, 2H>=, F8* 3, 3X, 1 2,2HD8) 

C WRITE OJT LAST VALUES OF THE TABLE 
WHl TEC 6, 1 01 3 TUCNfl >, TUO*tNTl > 

101 FORMAT ( 1 9X ,7HTU<42)=# F8* 3, 5X, 8HTU0C 42 ) = , F8 • 3,7H 5D0///) 

1ERR2-0 

C FINAL OUTPUT VALUE FOR THE LINEARIZATION TABLE FORMATION 

X F= 5 1 1*5 

C GET THE INPUT VALUE AND STORE IN "XF" 

CALL L I NAP * X F, N T 1 ) 

IF <XF«LE*0* )XF=512 

C NORMALIZATION FACTOR OF THE OUTPUT VALUES OF THE LIN* TABLE 

X I * 1 023 ♦ /XF 

C INITIAL OUTPUT USED TO DETERMINE THE INPUT VALUE 
Xl=5li *5/1024* 

C FIRST INPUT VALUE OF THE TABLE 
START U> -1 
DO 4 1=2,512 

C NEXT OUT PUT VALUE USED 
X3=Xl*FLOAT<2*I-l > 

c GET THE INPUT VALUE AND STORE IN M X3 M 
CALL LINAP <X3,Nri> 

X3=*X3+Kl+l *)/2« 

I FCX3.LE* 1 ♦ >X3=1 *1 

C STORE INPUT VALUE IN INTEGER LIN. TABLE 
5TAKT<l>*IFIX<X3> 

C ERROR CONDITION FOR LlNEAKl ZATI ON TABLE 
4 I F* START* I >+l *LT* START* 1-1 > >IERR2*I ERR2+1 

C LAST VALUE OF THE TABLE 
START* 5 1 2 > = 51 1 
I F( DLPO • LT * 1 ) GO TO 6 
C NEW PAGE 

WRITE *6, 103) 

103 FORMAT * 1 HI > 

C WRITE LINEARIZATION TABLE ON TELETYPE 
WRITE *6,104) C S T A RT *I>,I=1,5I2) 

104 FORMAT *10(15,2X3) 

C NEW PAGE 

WRITE* 6, 103) 

C WRITE TABLE ON A STORAGE DEVIVE (DUMP MODE) 

6 CALL TTM* START) 

C WRITE OUT ANY ERROR AND ALLOW RECALI BRATI ON IF NEEDED 
I FCIERK2.NE.0 3GO TO 9 
I F< I ERR -EQ.0) RETURN 
WRITE* 6, 1 06 > 

106 FORMAT * / /47 H ♦♦♦♦CHECH CALIBRATION FOR POSSIBLE ERRORS*^* /// > 
GO TO 7 

9 WHITE* 6, 1 07 > 

107 FORMAT* 37H +++ ERROR- -BAD CALIBRATION TABLE + + + > 

GO TO 7 

RETURN 

END 
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C *+**********+**+*4**tS^ H0U TiNE LINAP****+**+******+«***** 

C LI NAP TRANSFORMS OUTPUT VOLTAGES INTO INPUT VOLTAGES OF 

C THE RECEIVER* THE CALIBRATION DATA IS CONTAINED IN SUB- 

C ROUTINE VALUE. U 

c ♦*♦***♦«**♦♦♦*♦♦**«****♦♦**♦***««**♦♦♦♦*♦*♦*••**♦*♦«* 

C 

C INPUT AND OUTPUT! 

C A IS THE OUTPUT VOLTAGE THAT IS TRANSFORMED INTO 

C INPUT VOLTAGE 

C NC5 IS THE NUMBER OF DB SETTINGS 

C 

SUBROUTINE LINAP<A#NC5> 

COMMON /TA/ 5<43)#TU(44># TUQC44) 

N« 1 

C DIVIDE THE STRAIGHT LINE APPROXIMATION INTO A AREAS 
NC54«<NC5*2)/4 
NC53»NC5*3/4 
NC52b<NC5+1 ) /2 

C FIND WHERE THE INPUTED NUMBER LIES 
I FtA.GT *TUQ<NC54) >N=NC54 
I P(A«GT* TUO<NC52) >N a NC52 
I F< A.GT. TU0(NC53) )NoNC53 
C SET THE UPPER LIMIT OF THE SEARCH 
K-N+NC54 

C SEARCH FOR THE CORRECT LINE SEGMENT 
DO 5 I «N# K 
J«I + 1 

I FCA.QT*TUO(I ) >AND*A*LE. TUOC J) ) GO TO 10 
5 CONTINUE 

C LINE SEGMENT COULDN’T BE FOUND 
A=0* 

RETURN 

C GET THE VALUE OF THE CORRESPOND! NG OUTPUT 
10 Ao (A-TUOC I ) )*S< I > + TU( I ) 

KETU HN 
END 


CPPPPPPPPPPPPPPPPPPPPPPPPPPPP*-PHOC--PPPPPPPRpppppppppppppppppppppp pp 
C THIS PHD'IHAM f>PP»ATPS nwrt 1.1 t tu r>i r, r? < 


THIS PROGRAM OPERATES ALONG WITH DL03F1 . IT TAKES DATA READ 
BY A SUBROUTINE AND CHECKS THE NOISE AND SATURATIONS 
FOR GOOD DATA. IT SUMS THE SQUARES OF THE DATA AT EACH HEIGHT 
AND THE ACCEPTABLE NOISE# SUBTRACTS THE ACCEPTABLE NOISE SQUARED 
FROM THDATA# TAKES THE SQUARE HOOT OF THE DATA PRESERVING THE 
SIGN# TAKES THE SQUARE ROOT OF THE NOISE# AND PRINTS OUT THE RESULTS 
ON THE TELETYPE. THE RESULTS MAY ALSO BE PRINTED OUT ON PAPER TAPE 
IF NEEDED. THE PROGRAMS USED BY PROC ARES 

CHNS INITIALIZES STORAGE DEVISE AND WAITS FOR FILE 

TO BE COLLECTED PLUS COMMUNICATING WITH COLL. 
DHD7 3- -STORES DATA IN REAL ARRAY ( FORTRAN > 

CALC2 - “•'CALC ULATES ELECTRON DENSITIES < FORTRAN > 
CKCOL**- CHECKS FOR LEM VAN TED STOPS IN COLLECTION 
PROGRAM (MACRO) 

PP7 — CHECKS DATA SVICHES ON THE CONSOLE TO ALLOW 


— fl ^ HL.L, 

^ OH DISALLOW RESULTS TO BE PUT ON PAPER TAPE 

CPPPPPPPPPPPPPPPPPPPPPPPPPPPPpppppppppppppppppppppppppppppppppppppppp 


SUBROUTINE PROCC IA#PLF# I TIME) 

COMMON / PPC / A0C21 )#AX<2I ) # AVAO< 21) , AVAX ( 21) , I TI M < 4 > # 
lX0<2|)#IRJ(2n#BN0(5>#BNX(5)#EL(21 ), 

2HBMX# AVNO# BMO# A VNX# BMX # I D# IK 

COMMON /STAT/ I DB C 4) # RDATE ( 2) # REAS( 5) # I DBG# BMXNS# NCA 

SRBMX* 3.1 

RBMX* SRBMX* SRBMX 


IWIN=-1 


IEQU=0 
I SUM* 1 
NC4M 

C INITIALIZE PAPER PUNCH 
VRI T E(7 # 1 500) 

C INITIALIZE ALL VARIABLES NEEDED TO BE INlT. 
10 SNO&0 * 

SNX-0. 

IR«0 
IRN«0 
BMO«0 . 

BMX a 0 * 

DO 100 I * 1 # 2 1 
AVAOC I >~0. 

A VAX ( 1 > “0 ♦ 

AQ( I >*0 • 

AX ( I ) *0 * 

I RJ( I ) *0 


100 



109 


DO 110 1-1*4 
BNOCD-0* 

110 BNXCZ1-0. 

C INITIALIZE STORAGE DEVICE AND 
C PREPARE TO READ DATA 

CALL CHNG( 1BMX* I DBCH* IMCC) 

I F( IMCC#LE«0 • 0R*IMCC#QT*9)G0 TO 80 
CMC- FLOAT < I MCC > ~4 . 5 
RBMX-R8MX*CMC*CMC* . 4 

2 0 1FUBMX*LT«1 )I8MX-816 

IMCC-0 

C NOISE CRITERION 

c 1 BMX-TKE SOM OF THE SET OF 45 NOISE SAMPLES WHICH HAS THE 
C MIN. MAXIMUM OF ALL THE MAXIMUMS OF EACH SET OF 45 NOISE SAMPLES 
C THE AVERAGE NOISE FOR THE FIRST 5 IN EACH FRAME HAS 
C TO BE LESS THAN SRBMX* IBMX/45# WHERE RBMX IS THE SUPPLIED CONSTANT. 
C FOR SPEED RBMX* ( IBMX/45* )**2*5# IS COMPARED TO THE SUM OF THE 
C SQUARES OF THE NOISE# 

DUMX- FLQATC I BMX > / 4 5# 

I F( DUMX* SRBMX • GT« 500# ) DUMX- 500 • / SRBMX 

BMXN5-RBMX* DUMX* 5 . * DUMX 

KEOFO-0 

KEOFX-0 

ID-0 

C GET ORDINARY MODE DATA 

30 CALL DRD7 3( AO* BNO# I ERR# I Dj KEOFO) 

IFfKEOFO.EQ.l >30 TO 80 
C GET EXTHAORDINARY MODE DATA 

CALL DRD7 3* AXj BNX* I ERR* I D# KEOFX > 

IF<KEOFX*EQ. 1 >G0 TO 00 

C SET UP CHECK FOR REJECTION BECAUSE OF NOISE CRITERION 
BMEANO-0* 

BMEANX-0* 

DO 120 1-1*5 

BMEANO=BMEANO*BNO< I >*BNO( I > 

120 BMEANX«BMEANX* 8 NXCI>*BNXCn 

I F<BMEANO#G T# BMXNS# OR#BMEANX# GT# BMXNS>GO TO 50 
C NOISE USED TO SUTRACT FROM DATA SAMPLES 
BMO=BMO+BMEANO 
BMX-BMX+BMEANX 

c SUM OF THE SQUARE OF THE UNSATURATED DATA AT EACH HEIGHT 
DO 140 1=1*21 

35 I FC AOf I)*GE*510«# OH* AXC I)*GE.510#>GO TO 40 

AVAOC 1 3 = AVAO< l > +AOC I ) *AQ< I ) 

AVAX<I>=AVAX< I >*AX< I )*AX(I > 

GO TO 140 

C REJECTIONS DUE TO SATURATIONS OF DATA 
HQ IRJCI >=IRJ( I >+l 

140 CONTINUE 

GO TO 60 

C REJECTIONS FROM NOISE CRITERION 
50 IK»IH*1 

C SET UP AVERAGE NOISE USED IN REJECTION CRITERION 
60 SNO- 8 MEANO+SNO 

SNX-BMEANX+SNX 
GO TO 30 
80 I D= 1 D /2 

BID- ID* 5 

C MAXIMUM ALLOWABLE NOISE 

BMXNS-BMXNS/C DUMX*5# * SRBMX > 

C RMS OF ALL NOISE SAMPLES 
AVNO- SQRTC SNO/BI D> 

A VNX- SORT C SNX/BI D> 

C NUMBER OF ACCEPTABLE NOISE SAMPLES 
RN=5*<ID-IR> 

DO 150 1-1*21 

C NUMBER OF REJECTIONS AT EACH HEIGHT 
iaj<n«iRjm+iR 

c NUMBER OF ACCEPTABLE DATA AT EACH HEIGHT 

hsam-id-irju > 

c AVERAGE SUM SQUARED OF ACCEPTABLE DATA FOR EACH HEIGHT MINUS THE 

C AVERAGE SUM SQUARED OF THE ACCEPTABLE NOISE 

AVQC-AVAOC I > /RSAM-BMO/RN 
AVXC-AVAXC I >/RSAM-BMX/RN 

C THE RMS OF THE ACCEPTABLE DATA AT EACH HEIGHT (PRESERVING THE SIGN) 
AvAO ( l > = < ABS( AVOC ) /AVOC )* SQRTC ABS(AVOC) > 

AVAXC I >=CABS<AVXC ) /AVXC )* SQRTC ABSC AVXC) > 

EL< 1 ) -0 # 

XO< I >=0. 

Z F(AVAO< I ) «LE«0 .0 .OR. AVAXC I >#LE#0.0)GO TO 150 
XO ( I > -A VAX ( I > /AVAO C I > 
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no 


1 50 CONTINUE 

C THE RMS OF THE ACCEPTABLE NOISE 
BMO»SOAT(BMO/RN> 

BMX-SQRTOMX/RN) 

CALL CALC2IXQ. 1. 20. EL. I A) 

C GET THE TIME OF DAY 
DO IS5 t • l» A 
I1-5-I 

1 55 ITIM!I>-ITIME/!10**I1) 

C WRITE THE HEADING ON THE TELETYPE 

WRITE! 6. 1050 ) I TIMjRDATE, REAS. I D8!NC4> 

1050 FORMAT! 1 HI* 41 1 . AX* 2A5. 3X» SA5. 3X» I 2* 2HDB) 

WRITE { 6*1100) BMXN5* RBMX 

1100 FORMAT! //10H MAX. ALLOW. NOl SE»* F7 < I * 1 6H MULT. CONST. •* F7 • 3/> 

WR1TEI6* 1200) AVNO.BMO* AVNX.BMX 
1200 FORMAT ! 1 6H O-NOISE AU.!1) *F8.1*7H !2> *F6.l/ 

116H X-NOISE AU.!1) *F0.1*7H !2> *F8.l> 

WRITE16* J300>ID*IR 

1300 FOKMAT!//IX.I4.BH SAMPLES* SX* I 5* 1 2H REJ. !NOI SE) //3X* 

14HREJ./48H !N.+SAT.) HEIGHT AV. AO AV. AX AX/AO 
HT°5 8.5 
DO 160 I “ l * 2 1 

C CHECKS FOR COLLECTION STOPPAGE 
CALL CKCOL 
HT«HT+1.5 

WRITE! 6* 1400)IRJ!1>*HT*AVAO<I).AVAX!I>»XO! I )*PLF*EL!I> 

1 A00 FORMAT! AX* I A* 4X. F5. 1 . 3X* F6.1*2X*F6.1* 2X*F5. 2*A3*F6.0) 

| 60 CONTINUE 

C ALLOWS RESULTS TO BE SAVED ON PAPER TAPE 
CALL PP7 

C NEXT ATTENUATOR SETTING 
NC4-NC4+1 

IF(NC4.3T*IDBC)NC4=I 
IF< IDBCH.GT.0 >IDB(NC4)«5*IDBCH 
GO TO 10 

1500 FORMAT! 1H > 

RETURN 
END 


□ ************ + + **** ****** SUBROUTINE DRD7 3 ************************ *** 

DHEAD HEADS 21 SAMPLES OF SIGNAL AND 5 SAMPLES OF NOISE 
FROM DECTAPE* THE OUTPUT VOLTAGES HAVE BEEN TRANSFORMED INTO 
INPUT VOLTAGES* THE PROGRAM USES SUBROUTINE DUMPT (MACRO) TO 
HEAD DATA FROM STORAGE DEVICE. 

************ ************ *********** ************************ ********* 

SUBROUTINE DHD7 3( A# 8MEAN* 1 ERR* I D#KE0F) 

DIMENSION A C 21 I DATC 27 ) j 3MEAN (SI 

KEO F=0 

N-5 

N 1 = N* 1 

N2*N+2 

N3 = N 1 +2 1 

GET ONE SET OF DATA (26 NUMBERS) 

CALL DUMPTC I DAT* NE3 F) 

C CHECK ID CONSECUTIVE 

I KC I D- I DAT< 1 ) + l > 1 -I j 15* 10 
C CHECK FOR END OF FILE 
10 IFCIOATC1 >.NE. 130053) GO TO 20 

C TELL PROC IT’S THE END OF THE FILE 
KE0F=1 

return 

15 10-IDATC1) 

C SET DATA SAMPLES INTO A REAL ARRAY 
40 DO 42 MIN=N2,N3 

MKVE"MIN-NI 
A(MFVE>-I DAKmIN) 

42 CONTINUE 

C SET NOISE SAMPLES INTO A REAL N JMBEH ARRAY 
DO 130 U= 1 * N 
JEL=U+1 

1 30 DMSAN< J> = IDATCJ£L> 

RETURN 

C THE ID IS ERRONEOUS* IGNORE THE REST OF THE DATA 
20 VRI T£< 6* 1 JU) ID* I DAT C l ), IDAT<2> 

103 FOliMAT ( 44H ID UAS NOT CONSECUTIVE AND N0T=1 3.3-350 * I D=* 3( I 7 # 3X > ) 

KEO F= 1 
RETURN 
END 
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CCCCCCCCCCCCCCCCCCCCCCCCCCC— -CALCS— -cccccccccccccccccccccccccccccc 
c CALC2 IS A LIST OF CONSTANTS CALUATEO FROM THE PR 03 RAMS 

C CALC AND ELDEN AND CONTAINS THE FUNCTION THAT CALCULATES THE 
C ELECTRON DENSITIES FOR THE PARTIAL- REFLECTI ON PROCESSING PROGRAMS* 
C THE PROGRAM CALC tf RITES THIS PROGRAM* 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

SUBROUTINE CALC2 ( ARRAY* LL*LH#FD* IA> 

DIMENSION ARRAYC21)# RATI 02(21 >*FD<31 ) 

C GET THE PREDETERMINED CONSTANTS FOR THE RIGHT SEASON 
1FC IA >200 * 300* 100 
C CONSTANTS FOR THE SLWMER 
100 RATI02C D* 1.0731152 

RAT I 02 ( 2 )■ 1.0776633 
RATZ02C 3)- 1 *0841143 
RATI 02C 4> = 1.0909986 
RATI 02 ( 5> = 1.0990518 
RATI02C 6) = 1 .090 1243 
RATI 02 ( 7>« 1.0880302 
RATI02C 8 > = 1*0864129 
RAT 102 < 9)» 1.0768397 
RATIO2(10)= 1.0590323 
RATI 02C ll)^ 1.0495014 
RAT 1 02 C 1 2 > = i .0 328141 
RATI02C13)= 1.0262468 
RAT I 02( 14)° 1*0166022 
RAT 1 02C I 5 > 2 3 * 5 * * * * * * 12 13 14 * * 17 18 19 20 1.0117923 
RATI 02 C 1 6 ) =* 1 .0076428 
RATI 02C 17)° 1.0044388 
RATI02C 18) = 1.0029128 
RATI02(19)= 1.0015084 
RATI 02 (20 >» 1 .000844 3 
FD< 1>» 0 • 170 327 E-0 3 
FDC 2>» 0.243443E-03 
FDC 3 >■ 0.329813E-03 
FDC 4)- 0.427651E-03 
FDC 5> = 0 • 532098 E-0 3 
FDC 6)* 0 *627 51 3E-03 
FDC 7> = 0 *6997 30 E-0 3 
FDC 8 ) = 0.744879E-03 
FDC 9 ) = 0.7 53246E-0 3 
FDC 10)= 0.722645E-03 
FDC 1 1 ) = 0. 6 60 679 E- 03 
FDC 12)* 0 . 579908 E-0 3 
FDC 1 3> * 0.492934E-03 
FDC 1 4>* 0 . 40 577 2 E-0 3 
FDC 1 5 >■ 0.328260E-03 
FDC l 6 > = 0.258602E-03 
FDC 17) = 0.201748E-03 
FDC 1 8 > = 0*1 55429E-0 3 
FDC 19)= 0*11841 5E-0 3 
FDC 20 ) = 0.9I5612E-04 
GO TO 400 

C CONSTANTS FOR THE WINTER 
200 RATI02C 1 > = 1 .068 2487 


I .0880302 
1.0864129 
1.0768397 
1.0 590 323 


1 .0076428 
1.0044388 


HATI02C 
AATI02C 
RATI02C 
RATI 02 ( 
RAT 1 02 C 
RATI 02 C 
RAT I 02 C 
RATI02C 
JSATI02C 
RATt 02 C 
RATI02C 
RATI 02 < 
RATI02C 
HATI02C 
RATI 02 C 
RATI02C 
RATI02C 
RATI02C 
RATI 02 C 
FDC 1> = 
FDC 2> = 
FDC 3 ) = 
FDC 4> = 
FDC S) = 
FDC 6>» 
FDC 7> = 
FDC 8 ) = 
FDC 9> = 
FDC 10 ) = 
FDC 1 l )= 


2 ) * 1. 0890572 

3) * 1.0699695 

4> = 1.0979864 

5) = 1.0692204 

6>» 1 .0854059 

7>= 1.0770541 

8> = 1.0745673 

9 > = 1.0665843 

10)= 1.0531266 

lt> = 1.0458097 

12) = 1.0329590 

13) = 1.0286678 

14) = 1.0167225 

!5)= 1.0127123 

1 6> = 1.0076428 

17) = 1.0046618 

18) = 1.0023713 

19) = 1.0019426 

20) = 1.0007090 
0.237897E-0 3 
0 » 31 99 60 E-0 3 
0.410827E-03 
0.502861 E-03 
0 . 592S94E-03 
0.662S25E-03 
0.720342E-03 
0.750909E-03 
0 .7 510 34E-0 3 
0.720 29 3 E-0 3 
0 .663644E-03 
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FD<12>« 0.568398E-03 
FD(I3)« 0.S00S27E°03 
FD()4)» 0.410427E-03 
FD<15> 8 0.330916E-03 
FDC 1 6>*= 0.258602E-03 
F0< 17 )° 0 •200210 E-03 
FDC 1 0>o 0.156383E-03 
FD< 19 >** 0*1 19 390E»0 3 
FDC20 0.899351 E-04 
GO TO 400 

C CONSTANTS FOR EQUINOX 
300 RAT 1 02 C 1>= I *0670 392 

RATI 02 < 25= 1 *0803050 
RATI 02 C 3>« 1*07 80856 
RAT I 02 < 4> = 1.0912947 
RATI 02 C 5>« 1.0914300 
RATI 02 C 6>= 1 .0917337 
RATI 02 ( 7 )= 1.0897206 
RATI02C 8> = 1.0825763 
RATI 02C 9>= 1.0714581 
RATI 02( 10)= 1.0590323 
RAT I 02C 1 1 ) = 1.0482702 
RAT 1 02C 12)= 1.0351881 
RAT I 02C 1 3) a 1.0240435 
RAT 1 02 < 1 4) = 1.0176469 
RATI02C 1 5)*» 1 .0117923 
RATI 02 (16)= 1.0070148 
RAT I 02 C 17 )- 1.0048179 
RATI 02 < 1 8 > ° 1.0027836 
HAT 1 02( 19 ) = 1 .0015986 
RATI 02(20)= 1 .0007970 
FD< 1 >= 0.188515E-03 
FDC 2>o 0.261225E-03 
FDC 3)= 0.347 697 E-0 3 
FDC 4)* 0*4421 69E-0 3 
FDC 5 > s 0.542672E-03 
FDC 6)= 0 . 633402 E*0 3 
FDC 7 >= 0 • 705094E-0 3 
FDC 8 >= 0.747132E-03 
FDC 9 ) = 0*7 527 8 1 E-0 3 
FDC 1 0 ) = 0.722645E-03 
FDC 11)= 0 • 661 807 E-0 3 
FDC l 2 ) = 0.579908E-03 
FDC 1 3 > = 0*4929 34E-0 3 
FDC 1 4 >= 0.408J05E-03 
FDC 15)= 0.320260E-03 
FDC 1 6 >= 0.26I505E-03 
FDC 1 7 >- 0.236349E-03 
FDC 18)= 0.160195E-33 
FDC 1 9 >= *0.1 24908E-0 3 
FDC 20 ) = 3*997 9 59 E-04 
4 30 DO U I=LL# LH 

IFCAHHATC 1 ) *EQ.O. -QR.ARRAYC I ♦ 1 ) .EQ.3. )G0 TO 50 
C THE FUNCTION FOR THE CALCULATION OF ELECTRON DENSITIES 
FDC I ) sALOj C (ARRAY < 1 > /ARRAY C I M > ) * RATI 02 < 1 > > /FDC I ) 
30 TO 10 
50 FDC I ) =0 • 

1.) CONTINUE 

RETURN 


CPPPPPPPPPPPPPPPPPPPPPPPPRP~--WHPP-*-PPPPPPPPPPPPPPPPPPPPPPPPRPPPpppp 
C WRPP PUNCHES THE PROCESSED PARTIAL REFLECTION DATA 
C ON PAPER TAPE 

CPPPPPPPPPPPPPRPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPP 

c 

SUBROUTINE URPP 

COMMON /PPC/ AOC 21 ) * AX C 21 > # A VAOC 21 ) # AVAXC 2l)#ITlM<4)# 

1X0C 2 1 ) . I RUC 2 1 ) * BNOC 5 ) # BNX< 5) # EL <21 )# 

2RBMX# AVNO. BMO# AVNX# BMX# I D# I R 

COMMON /STAT/ I DB< 4) » RDATEC 2) * REASC 5) , I DBC* BMXN S# NC4 
WRITEC7# 10 50 ) I TlM#RDATE# REAS# 1 DBCNC4) #BMXNS# RBMX# 

1 AVNQ# BMO# AVNX# BMX* 1 D# I R 

1 0 50 F0RMAT<1H1.4I l*4X#2A5# 3X* 5A5# 3X# I 2# 2HDB/F7 . 1# F7 . 3# F8 . 1. 

1 F8 * 1 # F8 . 1 * F8 . 1 * I 4# I 5) 

DO 160 1 8 i # 2 1 
CALL CKCOL 

WRI TEC7 * 1400) IRJCI >*AVAO< I ) # AVAXC I >*XQ< I )#EL( I ) 

1 400 FORMA T< IX# I4#F6.1#F6. 1#F5.2*F6.0> 

1 60 CONTINUE 

RETURN 
END 



C ************ ************-pR0C7 3-***************+******************** 

C PR0C73 EVALUATES COLLECTED PARTI AL- REFLECTION DATA AMD 

C PRINTS OUT THE ELECTRON DENSITY. PR0C73 USES THE FOLLOWING PROGRAMS! 

C HEAD SETS IP AND PRINTS THE HEADING (FORTRAN) 

C DINIT-*" INI TI ALICES THE STORAGE DEVIVE (MACRO) 

C FSTAT— -LOCATES THE DATA FILE (MACRO) 

C SEEK FINDS THE FILE ON THE STORAGE DEVICE (MACRO) 

C DRD73--SETSSAMPLES INTO THE REAL / ARRAY (FORTRAN) 

C CALC2-— CALCULATES THE ELECTRON DENSITY (FORTRAN) 

C************* ************************************************ ****** 

c 

integer datin 

DIMENSION FNAMC 2 ># AOC 21 )#AXC2l ) *AVAO( 31 )#AVAX( 31 )* 

1X0(21 >*IftJ( 21 )*BN0<5)#BNX(5>*EL(2l ) 

WRI TEC 6 * 1 0 5 ) 

105 F0RMATCA8H TYPE IN SEASON— (1) FOH SUMMER* (*D FOR VINTER 

11SH* (0) OTHERWISE) 

READ (A* 200 ) I A 
200 FORMAT < 12) 

DATIN«2 

10 CALL HEAD< 0 ) 

C INITIALIZE VARIABLES 

12 SNO»0 • 

SNX a 3 * 

IR -0 

I RW»0 
BMO»0. 

BMX=0. 

DO 16 1 « 1 / 2 1 
AVAO C I )«0 • 

AVAX ( I ) «0 * 

16 IRJC I >«0 

DO 17 I»1*A 
BNOC I >=*0 • 

17 BMX( I ) *0 • 

C INITIALIZE TAPE STORING THE DATA 
CALL DINIT 

c get the data file name 

VRI TE< 6# 20 > 

20 FORMAT ( 1 SH WHICH DATAFILE) 

READC A* 33 >FNAM 
30 FORMAT (2A5) 

C CHECK THE VALIDITY OF THE NAME GIVEN 

13 CALL FSTAT< DATIN# FNAM#LGG> 

I F( LOG ♦ NE.0 >30 TO A 3 

WHI TEC 6 * 35> FNAM 

3 b FORMATC 6 H FILE #2A5*19H NOT FOUND ON DAT 2) 

30 TO 10 

C FIND LOCATION OF FILE ON THE TAPE 
40 CALL SEEK (DATIN* FNAM) 

C GET THE MAXIMUM ALLOWABLE NOISE 
WRIT E< 6 j 57 ) 

57 FORMATC 1 AH MAXIMUM NOISE) 

HEAD ( A* 56) BMXNS 


5 6 FOAmA^C FI 0 • 0 > 

I BMXNS* GE. SI 0 • >BMXNS“400. 

C FOR SPEED USE THE SQUARE OF THE MAX. ALLOW. NOISE TIMES 5 
DUMA** BMXNS* 5* 

BMXNS “BMXNS* DUMA 
1 9 KEOFO«0 

KEOFX = 0 
I D“0 

C GET ONE SET OF 26 ORDINARY SAMPLES 
48 CALL DRD7 3(A0*3N0#I ERR* 1 D* KEOFO) 

I F< KEO FO • EQ • 1 ) G 0 TO 50 
C GET ONE SET OF 26 EX TRAORDI NARY SAMPLES 
CALL DRD7 3< AX* 3NX* I ERR* I D* KEOFX ) 

IFCKEOFX.EQ.l >G0 TO 49 
C GET THE SUM SQUARED OF THE NOISE 
0MEANO=0. 

BMEANX=0. 

DO A40 I»l*5 

BMEAN0*3MEAN0+BN0< I ) * BN 0(1 ) 

440 BM£ANX=BMEANX*BNX< I >*9NX( I ) 

C CHECK FOR SETS OF SAMPLES THAT ARE TOO NOISY 

I FCBMEANO. GT. BMXNS. OR. BMEANX.GT.aMXNS)GO TO 510 
C SUM THE SQUARED NOISE SAMPLES FOR THE 
C LAST FOUR NOISE SAMPLES PER 25 TOTAL SAMPLES 
BMO=BMO+BMEANO-BNO(I)*BNO( 1 ) 

BMX«BMX*8MEANX-BNX ( 1 >*BNX< 1 ) 


730 DO A7 I « 1 # 21 

IFCAOC 1 > *GE. 510. .OR.AXC I > .GE. 510 . >G0 TO A6 
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C SUM OF THE SQUARED GOOD AO AND AX SAMPLES 
AVAOC I )«AVAOC I )+AO( I )+AO( I > 

120 AVAXCI>-AVAXCl>+AX<n*AX<I> 

GO TO 47 

C THE TOTAL i OF REJECTIONS DUE TO SATURATION 
C PLUS NOISE ABOVE THE GIVEN MAXIMUM 

46 1RJ(X)-IRJ<I)+1 

47 CONTINUE 
GO TO 520 

510 IR-IR+1 

520 SNO-BMEANO+SNO 

SNX-BMEANX+SNX 
GO TO 48 

49 ID-ID-1 

50 I D=I D/2 
BI D= 1 D* 5 

BMXNS=9MXNS/DUM4 

C THE RMS OF ALL NOISE SAMPLES TAKEN 
A VNO“ SQRTC SNO/B I D ) 

AVNX*SQRTC SNX/BI D> 

C THE NUMBER OF THE ACCEPTABLE NOISE SAMPLES 
RN«4*< ID-IR) 

DO 52 1=1*21 

C NUMBER OF REJECTED DATA AT EACH HEIGHT 
IRJCI>oIRJ<I)+lH 

C NUMBER OF ACCEPTABLE DATA AT EACH HEIGHT 
RSAMoID-IRJCI ) 

C ACCEPTABLE NOISE IS SUBTRACTED OFF 
AVOCoAVAOO >/R5AM-BM0/RM 
AVXC“AVAXCI)/RSAM-BMX/KN 
C RMS OF GOOD DATA WITH THE SIGN PRESERVED 

AVAOt I ) o ( ABSC AVOC ) /A VOC > ♦ SQ RT< AB S < AVOC > > 

A VAX C I > = C ABS(AVXC) /AVXC > *5QHT< ABS< AVXC) > 

ELCI )«0 
XOCI >=0 

IFCAVAO(I).LE.0.0.OR.AOAXCI).LE,0*0>GO TO 52 
XOC I >«AVAX( I > /AVAOC I > 

52 CONTINUE 

C HMS OF ACCEPTABLE NOISE 
BMO=SQKT<BMO/RN) 

BMX » SQ RT< BMX /HN ) 

C GET ELECTRON DENSITIES 

CALL CALC2CX0# 1 # 20# LL* I A > 

C 'WRITE THE HEADING 
CALL HEADU) 

WRITE (6# 100) BMXNS# A VNQ# BMO# AVNX* BMX 
100 FORMAT (25H MAXIMUM ALLOWABLE NOI SE = »F6» 1//I6H Q-NOISE AV.Cl) j 

l F8 * 1 * 7 H (2) * F8 • 1 / 1 6H X-NOtSE AU.(l) #F8.1*7H (2) # F8*1> 

[ WKI TEC 6* 54) I D# ! R 

54 FORMAT C //IX# I 4# 8H SAMPLES# 5X# l S* 1 2H REJ. C NOI SE ) //8H REJECTS* 

1 2X# 6HHEI GHT# 2X# 6HAV* AQ*2X#6HAV* AX# 2X# 5HAX/A0# 4X# 2HED) 

HT= 58*5 

C WRITE OUT TABLE 
DO 53 1=1*21 
HToHT+1 .5 

WKITEC6* 58)1 HJ< I >#HT#AVAO< I )#AVAXCI ),XO< I )*EL< I ) 

58 FORMAT! 3X* 14# 3X» F5* 1 # 3X# F6* 1 # 2X* F 6*1# 2X# F5# 2# 3X# F6* 0 ) 

53 CONTINUE 
30 TO 10 
STOP 
END 


CHHHHHHHHHHKHHMHHHHHHHKHHH HEAD — 

HEAD SETS UP AND PRINTS OUT 
HEADING FOR A NEW PAGE. 


-MHHHHHHHHHHHHHHHHHHHHHHHHHHHHKH 
ONE LINE OF INFORMATION AS A 


CHHHHHHHHHHHHHHHHHHHHHMHHHHHHHHHHHHHHHHrlHHHHMHHHHHMHHHHHHHHHHHHHHHM 

C 


SUBROUTINE HEAD (J) 
REAL R EAS (12) 


I F (J ) S,6,lfl 
t* VRITE(C,25) REAS 

FORMAT (IH1 /I X, 1 2A5) 
R ETUR N 


<5 VRITE(S,!(W> 

10?I F0RMATOHI/27H SIVE ! LINE OF INFORMATION) 

PEAD(4,205t> REAS 
2™ F0RMAT( I2A 4) 


RETURN 


END 
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♦ TITLE READ DATA IN DUMP MODE 

/muuui nmnnmmnmtuninxmmummmnmmmsx 

/ .READ# DUMP MODE FROM DECTAPE ON A VARIABLE .DAT SLOT 

/FILLS 252 DEC VORD BUFFER AND OUTPUTS 26 

/WORDS TO ARRAY I DAT EVERY TIME CALLED- 

/THESE ARE UNPACKED FROM 18 WORDS OF THE BUFFER. 


/ I DATS 

WORD \ 


I • D. # 


/ 

WORD 2- 

6 

NOISE SAMPLES 

/ 

WORD 7* 

27 

DATA 


/ NEGFt 


SET IF 

A NEGATIVE NUMBER WAS IN THE DATA 

/ 

-GLOBL 

DIN 1 1# DUMP T# • DA 


DUMP *4 




/TYPE OF I/O MODE 

DATIN»2 




/.DAT SLOT TO READ DATA FROM 

NSAM-32 




/# OF SAMPLES PER SET 

N SAMpa.NSAM/2* l 



/SIZE OF ONE SET PACKED DOUBLE 

DTBLK»374 



/SIZE OF ONE BLOCK OF STORAGE 

datstr® 

dtblk/nsamp 


/# OF SETS PER ONE BLOCK OF STORAGE 

NDPC-NSAM/2 



/NUMBER OF STORED DATA PAIRS PER SET 


.IODEV 

DAT IN 



DINIT 

0 





♦ I NI T 

DATIN# 0 

#DINIT 

/INITIALIZE DEVICE STORING THE DATA 


LAW 

-1 


/PREPARE TO READ 


DAC 

CTT9 


/ IN ONE BLOCK OF DATA 


LAC 

BUF31 


/RESET THE BUFFER POINTER WITH 


DAC 

POINT2 


/ THE ADDR. OF BUF3 


LAC 

(ISZ 

SWI TC 

/PREPARE TO READ TWO 


DAC 

LCA 


/ DUMMY BLOCKS 


JMP* 

DINIT 


/END OF INITIALIZATION 

DUMPT 

0 





JMS* 

• DA 


/PICKUP ADDh OF ADDR 


JMP 

* + 3 


/OF ARRAY 

A2 

0 




FLAG 

0 



/SET ON NEG # 


LAC* 

A2 


/GET ADDR. 


DAC 

A2 


/ OF ARRAY 


LAW 

-NDPC 


/SET COUNTER OF DATA TO BE 


DAC 

COUNT 


/ PROCESSED 


ISZ 

CTT9 


/GET POINTER 


JMP 

LBB 


/NO# CONTINUE WITH PRESENT SET OF DATA 


LAW 

-DATSTft 


/RESET COUNTER TO THE NUMBER 


DAC 

CTT9 


/ OF SETS PER BLOCK OF STORAGE 

LBA 

♦ READ 

DATIN# DUMP# BUF3 j 

> DTBLK /GET 1 BLOCK OF DATA 

C14 

♦ WAIT 

DAT IN 



LCA 

ISZ 

SWITC 


/INITIALLY READ TWO DUMMY BLOCKS 


JMP 

LBA 




LAW 

-3 


/RESET CONTROL TO READ 


DAC 

SWITC 


/ TWO DUMMY BLOCKS 


LAC 

(JMP 

LCB 

/READ ONE BLOCK OF DATA 


DAC 

LCA 


/ AT A TIME 

LCB 

LAC 

BUF31 


✓GET ADDRESS OF f3UF3 


DAC 

POINTS 


/ POINT2 TO BUF3 

LBB 

LAC* 

POINT2 


/GET THE ID < FIRST WORD IN DATA SET) 


SAD 

SEVN 


/END OF FILE ID? 


JMP 

ENF 


/YES* RESET PARAM' S AND CLOSE FILE 


LAC* 

P0INT2 


/GET ID AND PUT 


DAC* 

A2 


/ INTO THE FORTRAN ARRAY 


ISZ 

A2 


/GO TO NEXT ADDR. OF THE ARRAY 

LOOP 

ISZ 

POINTS 


/GO TO NEXT DATA WORD 


LAW 

-2 


/PREPARE TO LIMPAC 


DAC 

TC10 


/ TWO DATA WORDS 


LAC* 

P0INT2 


/GET DATA WORDS FROM 8UF3 


SWHA 



/FIRST WORD IN LEFT HALF 

UNPLP 

AND 

(777 


/5AVE ONE DATA WORD 


SNA 



/CHECK FOR NEG. NUMBER 


ISZ* 

FLAG 


/SET IF NEG. NUMBER FOUND 


DAC* 

A2 


/LOAD # INTO FORTRAN ARRAY 


LAC* 

POINT2 


/GET DATA VORD AGAIN 


ISZ 

A2 


/GO TO NEXT LOC. IN ARRAY 


ISZ 

TC 1 0 


/UNPAGED TWO WORDS? 


JMP 

UNPLP 


/NO# LOOP AROUND 


ISZ 

COUNT 


/YES# HAS 34 DATA WORDS BEEN UNPACKED? 


JMP 

LOOP 


✓NO# REPEAT UNPACKING PROCESS 

0UT2 

ISZ 

POINTS 


/YES# GO TO NEXT ID 


JMP* 

DUMPT 


/RETURN 

SWITC 

LAW 

-3 



COUNT 

LAW 

-NDPC 



CTT9 

LAW 

•DATSTH 



SEVN 

37 6002 




TCI 0 

0 




P0INT2 

♦ DSA 

BUF3 



BUF3I 

. DSA 

BUF3 



BUF3 

•BLOCK 

DTBLK 
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/ END OF FILE ROUTINE 


ENF 

C 1 5 
RT2 
/ 


LAC 

DAC* 

.CLOSE 

JMP* 


SEVN /SET LAST ID TO 

A2 / 130050 DECIMAL 

DATIN /CLOSE FILE 

DU1PT /RETURN TO PROC* PROGRAM 


/mxxmxiunnmnmxnxnmmmmxiiimiumxmmxnm 


• END 


C***+* *********** ************* CALC ****************************** 

c FROM GIVEN COLLISION FKEQUENCIESjCALC ALONG WITH ELDEN 

c calculates the constant values used in the electron density 

C EQUATION GIVEN 3Y PIHNAT IN AERQNOMY REPORT 29 AND WRITES THE 

C PROGRAM CALC2 WHICH CALCULATES THE ELECTRO DENSITIES FOR THE 

C PARTIAL-REFLECTION PROGRAMS. 

C ****** *************************************************** ********* 

DIMENSION ARRAY <21 >*P<21 >,H< 3)#CF< 3) . EL ( 20 ) j CALC2C 2> 

DATA CALC2U > , CALC2C 2 ) /5HCALC2# 4H SRC/ 

IDAT=1 

LABL=0 

C COLLISION FREQUENCY PROFILES 

C COLLISION FREQUENCY PROFILE FOR THE SUMMER 

100 P< I >o 192* 3 

P<2>*156.9 
P<3)=127.5 
P<4>* 132.7 
P<5>=82.37 
PC 6)® 66*25 
P < 7 > = 52 ■ 53 
P<8)=4I .66 
PC9>=32.S1 
Pn«>=*25.04 
PC11)»20.1 
PCI 2>=15. 53 
PC 1 3) » I 1 .89 
PC14)=9.057 
PC 1 5 > = 6.3 17 
PCI 6) «5 • 399 
PCI7>=3.827 
PC 18)=2*S62 
PC 19)=2.124 
PC2U)sl . 563 
P<21 )=1 .180 

C WRITE THE PROGRAM HEADING ONTO TAPE 
CALL EN f ERC I DAT# CALC2 > 

WHITECOAT, 10 ) 

1 3 FORMAT < 6 9H CCCCCCCCCCCCCCCCCCCCCCCCCCC C ALC2---CCCCCCCCCC 

1 CCCCCCCCCCCCCCCCCCCC /59H C CALC 2 IS A LIST OF CONSTANTS 

2 CALUATED From THE PROG HAMS/ 6 3H C CALC AND ELDEN AND CONT 
3 AIMS THE FUNCTION THAT CALCULATES THE/69H C ELECTRON DENSITIES 
4 FOR THE PARTIAL-REFLECTION PROCESSING PROGRAMS. /43H C THE PRO 
5GHAM CALC WRITES THIS PRO 3 RAM. /69H CCCCCCCCCCCCCCCCCCCCCCCCCCC 
6CCCUCCCCCCCCCC CCCCCCCCCCCCCCCCCCCCCCCCCCC/ 37 H SUBROUTINE 
7 CALC2CARRAY, LL.LH, FD, IA>/39H DIMENSION ARHAYC 21 > j RATI 02 
8(21 ) # FDC 21 >/57H C GET THE PREDETERMI NED CONSTANTS FOR THE 
9 RIGHT SEASON/ 1 9H I FCI A > 233 * 330i 1 03 / 

128H C CONSTANTS FOR THE SUMMER) 

30 TO 430 

C COLLISION FREQUENCY PROFILE FOR THE WINTER 

200 PC1J-133. 5 

P<2) = 107 *8 
P<3>=87. 12 
P<4>=70.04 
P<5>«56.33 
P ( 6 > = 45 • 28 
P (7 ) = 36. 55 
P(g)=29.32 
P(9)=23. 52 
P <1 3 > =* 1 8 • B3 
PCI 1 >=14.97 
P < 1 2 > = ! 1.99 
PCI 3> = 9 . 561 
P < 1 4) =7 *541 
PU 53 = 6.038 
PCI 6>=4.748 
PCI7)=3.758 
P< 1S)=2.941 
P< 19)=2.321 



1 1 


c 

300 


1 2 

C 

400 

401 

C 


C 

20 

C 

405 

C 


25 

410 

G 

30 

420 

C 


C 

40 


P<20)»1.810 
P C 21 )• l *431 
VRITECIDAT,U> 

FORMAT (UH 00 TO 400/ 

1 20H C CONSTANTS FOR THE WINTERS 
30 TO 400 

COLLISION FREQUENCY PROFILE FOR THE EQUINOX 

PU >"160*2 

P(2>*130.2 

P ( 3>«105*3 

PC4)-S4.90 

PC5)-68*25 

P<6>-54.75 

PC7)-43*57 

PC8)«34«31 

P<9>»27.07 

PU0)>ei .32 

P < 1 1 )■ 1 6*62 

PC 12)“1 3* 26 

PC 1 3> = 10. 33 

PC 14>"8*062 

P ( 1 5 )° 6*246 

PC I 6)o4*835 

PC 17 > = 3 • 7 58 

P ( 18 > »2*9 1 5 

PC 19)«2* 260 

PC20 )«1 *733 

P (21 )*1 .359 

WRI TEC I DAT# 12) 

FORMAT! 1 1 H GO TO 400/ 

12SH C CONSTANTS FOR EQUINOX) 

SET THE COLLISION FREQUENCIES TO THE RIGHT ORDER OF MAGNITUDE 
DO 401 1=1,21 
PCI>»PCZ)+CI0,**5) 

CONTINUE 

STATEMENT LA0LE FOR THE NEW PROGRAM 
LABL-LABL+100 
Jl-Jl+4 
K=0 

DO 20 1*1,20 
K-K+l 

CFC 1 >=P(K> 

CF(2)=P(KM > 

CALCULATE CONSTANTS FOR THE ELECTRON DENSITY EQUATION 
CALL ELDENCR, OF* ARRAY < I >,AHKAY(I + 1 ) # EL C I > ) 

ARRAY C 1 ) = ARRAYC 2 ) /ARRAY Cl) 

WRITE FIRST CONTSTANT WITH A STATEMENT LABLE 
WKITECI DAT, 40 5>LABL, AHRAYC 1 ) 

FORMAT C 1 H ,I3,12H HATI02C 1>«,F10,7> 

WRITE THE REST OF THE RO AND RX CONSTANTS 
DO 25 1=2,20 

ARRAYC I )= ARRAY ( I + 1 >/AHHAY< I ) 

WRI TEC I DAT ,410) I » ARRAY Cl) 

FORMAT C 9H RATI 03C , 1 2, 2H >*, FI 0 ♦ 7 ) 

WRITE THE CONSTANT DENOMINATORS 
DO 30 1=1,20 
WRITE! I DAT, 420 ) I,ELU> 

FORMAT C 5H FDC, I2,2H>«, El 3*B> 

CALCULATE the REST OF THEC0NSTANT5 
IFC J1 .LT*5)<30 TO 200 
I FC Jl .LT. 10)QO TO 300 

WRITE THE REST OF THE PROGRAM TO CALCULATE ELECTRON DENSITIES 
WRI TEC I DAT, 40) 

FORMAT C 1 6H 400 DO 10 I=LL*LH/48H I FC ARRAYC I ) ♦ EQ . 0 . • OR. 

1 ARRAY < I+l ).EQ«0*)GO TO 50/S9H C THE FUNCTION FOR THE CALCULA 
2TI0N OF ELECTRON DENSITIES/ 

351 H FDC I )=ALOGC C ARRAY ( I > /ARRAY C I ♦ 1 > >*RATI 03C I ) >/FD< I > / 

410H GO TO 10/12H 50 FD< I >=0./12H 10 CON TI NUE/8H RETURN) 

CALL CLOSEC I DAT) 

STOP 

END 


Cl 1 1 « * ! f 1 1 1 > f 1 1 1 1 » i « t » -S UBR OUTI NE FLDEN- 1 1 1 1 f 1 1 1 1 til lit 1 1 1 1 II 1 1 1 1 It III 
C DURING DATA PROCESSING THERE ARE ONLY 2 VARIABLES 

C FOR EACH HEIGHT (AO AND AX). THE EQUATION FOR THE 
C ELECTRON DENSITY AS GIVEN BY BIRLY (I97|> ISi 

C ED=LNC((AX(I>/O0n))/(RXCI)/R0<!)))/(CAX(2)/A0<2))/CRX(2)/R0(2))>) /FD 
C WHERE LN IS THE NATURAL LOG AND I AND 2 ARE HEIGHT I AND 2 
C SUBROUTINE ELDEN CALCULATES THE CONSTANTS RX,RO,AND FD 

C FOR EACH HEIGHT. 

CtfittitiiiitiiiMMMtiiiiiiMitiiitttttitiittittiitumtmimiiii 
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SUBR OUTI NE ELDE N (AXBYA 0 p 6 NU ,R XRO 1 ,R XROS , FD ) 

DIMENSION AXBYA0(3) # RXBYR0(3)|RX(3) ,R 0(3) ,GNU(3) f ftATI 0 (3) 

C APPROX INTEGRAL PARAMETERS 

A4=2,3983474E-2 
A3r | , 1 2875 1 3E+ 1 
A2rl.l394160E+2 
A I =2 ,4653 1 1 5E+I 
B 6r 1 .80 64 1 28E-2 
B5= 9. 38 7 73 72 
B4 - J .492I254E+2 
B3r2.8958085E+2 
B2 = l ,20495 12 E*2 
B I =2,465681 9 E-H 
D3r| .1630641 
D2 = l .6901 002 E-f 1 
01=6.6945939 
E5 = 4. 3605 732 
E4 = 6. 4093464m 
E3 = 6 . *92 05 05 E+ 1 
E2=3. 535525 7E+! 

El =6 . 63 1 4497 
AXBYA0(3) r 0 

C G NU ( 3 ) IS MEAN COLLISION FREQUENCY AT THE INTERMEDIATE HEIGHT 

C CALCULATE C INTEGRALS AT BOTH HEIGHTS AND FOR AVERAGE GNU 

DUM = GNUU)+GNU(2) 

GNIK3)=0.5*DUM 
DO 22 X=I,3 
Or (2*5961 4E+7) /GNU(K> 

Xr7,3886E+ 6/GNU ( X) 

CT Nr 0* <0* (O CO+A1 >+A2>+A3>+A4 

CTD=O*C0* (0*C0* (0* CO+B 1 >+B2 5 +B3 )+B 4)+B5 )+B6 

CTO=CTN/CTD 

CTXNr X* (X* (X* (X+A 1 ) +A2 ) +A3 ) +A 4 

CTXD rX* (X* (X* (X* (X* (X+B l ) +B2 ) +B3 )+B4 )+B5) +R6 

CTX =CTXN/CTXD 

CFOr (0* ( O (0+D1 )+D2)+D3 ) /CO ( 0+ ( 0* CO* < 0+ El >-KF2 ) +E3 )+E4 ) +E5 > 
CFX=(X* (X* (X+Dl ) +D2) +D3 ) /(X* ( X* (X* (X* CX+ El ) + F2 ) +E3)+E4)+E5) 

C CALCULATE RATIOS 

RX<X>rSQRT< (X*CTX> **2+ (2 .5*CFX)**2> 
F0(K)rSQRT((0*CT0)**2+(2.5*CF0)**2) 

RXBYRO(K)rRXCX) /R0(K) 

RATIOCK) rAXPYAO(K) /RXBYRO(K) 

22 CONTINUE 

C CALCULATE FD FROM FINAL VALUES OF DO LOOP 

F0r(5.*3.| R24E+3*CF0) / ( 4 .*3 , 0F+ p* GNU C3) ) 

FX = (5.*3.1 B24E+3*CFX) /< 4 . *3 . 0 E+P* GNU (3) ) 

FD r ( FX -F 0) * 3 ♦ 0 E+ 9 
RXROIrRXBYRO(J) 

P XR 02 rR XBYR 0 (2 ) 

RETURN 

END 


C* ************************* ~ Off CK „********************************* 

C PROGRAM READS IN THE NUMBER OF SAMPLES ASKS FOR BY OPERATOR. 

C IF THE NUMBER IS ZERO , 3J NUMBERS ARE READ I N A ND SET UP AS PARTIAL 
C REFLECTION DATA IS <1 ,E. 5 NOISE SAMPLES AND 21 DATA POINTS A 5 EXTRA) 
C DATA IS PRINTED OUT IN THE FORM OF ONE NUMBER PER HEIGHT AFTER EACH 
C GROUP OF 26 SAMPLES ARE READ IN. THIS HAS BEEN USE TO CHECK THE 
C RECI EVER AND A/D CONVERTER AGAINST THE REYNOLDS SYSTEM AND TO SEE 
C IF EVERYTHING IS OPERATING AS IT SHOULD. IE THE NUMBER READ IN IS 
C NOT ZERO, THAT NUMBER OF SAMPLES ARE READ FROM THE A/D CONVEFTER 
C AND AN AVERAGE OF ALL THE NUMBERS ARE TAKEN AND PRINTED OUT. 

©************************ i**.****************^*^*^,^^^^^^ **** 

C 

DIMENSION IA(50),RAI (50), RAP (50) 

MA X -5 0 
131 =3t 

WRITE (6,110) 

110 FORMAT ( ] PH ADC CHECK) 

DEFAULT VALUE FOR THE 0 OF SAMPLES = 31 
NS =13 1 

READ 0 OF SAMPLE TO BE READ FROM A/D CONVERTER 
READ (4,210) IDV 
10 FORMAT(1 5) 

IF (IDV.NE.S) NS=I OV 
I F ( NS , NE .1 3 1 ) GO TO 50 
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FOR 31 NUMBERS READ IN, THE FORM USED IS 2 SETS OF 31 SAMPLES 
AS IN THE PARTIAL REFLECTION COLLECTION 
11:0 
ICH :0 
11=11+1 

CALL INPADCIA , NS, ICH) 

IF<ICH.EQ.0> SO TO 6 
iFdi.oT.nno to m 

CONVERSE OR 1 ALflORTH OM FOR A/D CONVERTER NEG. *’S TO COMPUTER 

NEGA I F <IA (I nOT^S ! 1 )I A (I > =3072+(409S+327SS)*7+IA (I ) 
RA1U)=FL0AT«A(I>)/.511 
00 TO 25 

DO L00p T F0R SECOND SET OF NUMBERS READ IN 
DO 15 1=1,131 

IF Cl A (I) ,6T ,51 1 >IA (l):3072+<4096+3276B)*7+IA<I ) 

IFCI^EO J5?0p!! .EO.l nVRlTE(6, 105) 

IFCI .E9.11)WRITF<6,10«) 

R A2(I) =FLOAT (IA(I))/.5|1 
WRITE OUT THE NUMBERS IN AN ORDERLY WAY 

WRITEC6, 1 20) HT,R A I (I) ,RA2 ( I) „ 

0 FORMAT C3X, F4. 1 ,4HKM ,F5.0,4HMV ,F5.0,2HMV> 

1 FORMAT (SH NOISE! 

? FORMAT (25H 

6 FORMAT (5H DATA) 

THE FOLLOW! NO DUMPS THE AVERAGE OF THE A/D CONVERTER NUMBERS 
AND ALSO GIVES THE VALUE IN MILLIVOLTS 
I NS: (NS+ MAX-1) /MAX 
TNS=I NS* MAX 
D 0 60 U = 1 , 1 NS 
I CH =0 

CALL INPAD CIA, MAX, ICH) 
i IFCICH.EO.0) GO TO 20 

I FC I A (I ) 1gT A .5 1 1)1 A CD =3072+14096+32762) *7+1 A C I ) 

A V=AV+ FLOAT (IA CD ) 

I CONTINUE 

A V=A V/TNS 

AVV-AV/.512 

><* FORMAT ( 9H AVERAGF- t F7 .3, 1 2H VOLTAGE^ f Fft«3 f 3M MV) 

GO TO 5 

STOP 

FND 


/ 

/ 

/ 
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.TITLE A/D CONVERTER SERVICE ROUTINES JO* M ; -FG. 

PFKM15 V3 A SERVICE ROUTINES FOR THE HP ( IfA , f 
1 ONVFftTFP THFBF ROUTINES PERMIT INPUT OF ANY SPECIFIED 
Sp OF SAMPLES INTO A COPE BUFFER. INPUT MAY BE OVER- 
LAPPED WITH PPOGPAM FYCUTION, AND COWTPOL MAY NOU I SHED 

•0 LOWER PRIORITY PROGRAMS WHILE DATA TRANSFER TAKES PLACE. 
MACPO-15 CALLING SFQUFNCE: 
jms I NPAD 

NUMBER of SAMPLES RFQUIPED 
BUFFER ADDRESS 

PFAL-TI MF^SUBROUTI NF^ADDPFSS , PRIORITY LEVEL IN BITS *-Z 
(FXAMPLF: 500PW0+PTSUBA) 

(RFTURNS HEPF IMMFDIATFLY) ^ ^ IS H, NO REAL-TIME SUBROUTINE 

WILL BF ACTIVATFD . NOTFs THF “RIORITY CODF FOR MAINSTREAM IS 1 
1 ?HF COMPUHON FLAG IS CLFAREP BY THE CALL TO INPAD. 

AND SET TO +1 FOR NORMAL COMPLFTION OR -100! IF A DATA 
TIMING FRROP OCCURS. 


/ 

ApWCR :26 
ADCAR :ADWCP+1 
.SCOM = 100 
ADVI : 703724 
ADSO=70378| 
ADST =70372 1 
ADCO:703784 
ADCT:703744 


/A-D WORD COUNT 

/AND CURRENT ADDRESS REGISTERS 
/MONITOR'S COMMUNICATION AREA 
/A-0 CONVFRTER WRITE INITIALIZE 
/SKIP ON WORD COUNT OVERFLOW 
/SKIP ON DATA TIMING ERROR 
/CLEAR OVERFLOW FLAG 
/CLEAR TIMING FLAG 
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ENTRY POINT FOR A -0 INTERFACE INITIALIZATION 


I WF 
I NWC 
I NFLAG 
I NR 


.GLOBL INPAD,. DA 
0 

JMS* .DA 

JMP .+4 


jmp 

TCA 

INSET 

/REPLACED BY ’'LAC* I NWC 1 

DAC* 

(ADWCR) 

/SET WORD COUNT 

LAW 

-1 


TAD* 

INAR 

/BUFFER ADDRESS -t 

DAC* 

(ADCAR) 

/ TO CURRENT ADDRESS REG • 

DZM* 

INFLAPS 

/CLEAR FLAG 

DZM 

INSUB# 

/CLEAR REAL-TIME SUBROUTINE 

ADWI 


/INITIALIZE INTERFACE 

JMP* 

INPAD 

/RETURN 

FOLLOWI NR CODE IS EXECUTED ONLY ONCE 

LAC* 

(.SCOM+55) 

/OFT ENTRY POINT ADDERSS OF 

DAC 

m 


JMS* 
ADSO 
ADI NT 

.-i 

/CALL .SETUP TO CONNECT ADI N, 

DZM* 

<204 


LAC 

(LAC* I NWC 


DAC 

INR 

/MODIFY INSTRUCTION 

JMP 

INR 

/AND JUMP TO IT 


/INTERRUPT SERVICE ROUTINE. FXECUTFP IMMEDIATELY AFTER COMPLETION 
/ OF DATA TRANSFER. DFTERMI NES STATUS OF A-D INTERFACE. SETS 
/ COMPLETION FLAO AND ACTIVATES REAL-TIME SUBROUTINE. 

/ PUNS AT API LEVEL 0. ' 


DBA 


/PA OF ADDRESSI NO MODF 

DAC 

ADSVA 

/SAVE AC 

ADST 


/TIMINO ERROR? 

SKP !CLA ! IAC 

/NO.+ I TO AC 

LAW 

♦!00i 

/YES, ERROR CODE 

DAC* 

INFLAG 

/SET FLAO 

AD CO 


/CLEAR 

ADCT 


/ INTERFACE FLAGS 

LAC 

ADSVA 

/RESTORE AC 

DBR 


/SET TO LEAVE HARDWARE 

JMP* 

ADINT 

/RETURN TO INTERRUPTED 



